We present a new simple algorithm for efficient, and relatively accurate computation of the Faddeyeva function w(z). The algorithm carefully exploits previous approximations by Hui et al [1978] and Humlicek [1982] along with asymptotic expressions from Laplace continued fractions. Over a wide and fine grid of the complex argument, z=x+iy, numerical results from the present approximation show a maximum relative error less than 4.0x10-5 for both real and imaginary parts of w while running in a relatively shorter execution time than other competitive techniques. In addition to the calculation of the Faddeyeva function, w, partial derivatives of the real and imaginary parts of the function can easily be calculated and returned as optional output.

This remark describes efficiency improvements to Algorithm 916 [Zaghloul and Ali 2011]. It is shown that the execution time required by the algorithm, when run at its highest accuracy, may be improved by more than a factor of two. A better accuracy vs efficiency trade off scheme is also implemented; this requires the user to supply the number of significant figures desired in the computed values as an extra input argument to the function. Using this trade-off, it is shown that the efficiency of the algorithm may be further improved significantly while maintaining reasonably accurate and safe results that are free of the pitfalls and complete loss of accuracy seen in other competitive techniques. The current version of the code is provided in Matlab and Scilab in addition to a Fortran translation prepared to meet the needs of real-world problems where very large numbers of function evaluations would require the use of a compiled language. To fulfill this last requirement, a recently proposed reformed version of Humlicek's w4 routine, shown to maintain the claimed accuracy of the algorithm over a wide and fine grid is implemented in the present Fortran translation for the case of 4 significant figures. This latter modification assures the reliability of the code to be employed in the solution of practical problems requiring numerous evaluation of the function for applications tolerating low accuracy computations (<10-4).

We present a simple reform for treating the reported problem of loss-of-accuracy near the real axis of Humlicek's w4 algorithm, widely used for the calculation of the Faddeyeva or complex probability function. The reformed routine maintains the claimed accuracy of the algorithm over a wide and fine grid that covers all the domain of the real part, x, of the complex input variable, z=x+iy, and values for the imaginary part in the range y=[10-30, 10+30]