The Voigt and complex error function: Humlíček's rational approximation generalized. Schreier, F. Monthly Notices of the Royal Astronomical Society, June, 2018.
The Voigt and complex error function: Humlíček's rational approximation generalized [link]Paper  doi  abstract   bibtex   
Accurate yet efficient computation of the Voigt and complex error function is a challenge since decades in astrophysics and other areas of physics. Rational approximations have attracted considerable attention and are used in many codes, often in combination with other techniques. The 12-term code “cpf12” of Humlíček (1979) achieves an accuracy of five to six significant digits throughout the entire complex plane. Here we generalize this algorithm to a larger (even) number of terms. The n = 16 approximation has a relative accuracy better than 10−5 for almost the entire complex plane except for very small imaginary values of the argument even without the correction term required for the cpf12 algorithm. With 20 terms the accuracy is better than 10−6. In addition to the accuracy assessment we discuss methods for optimization and propose a combination of the 16-term approximation with the asymptotic approximation of Humlíček (1982) for high efficiency.
@article{schreier_voigt_2018,
	title = {The {Voigt} and complex error function: {Humlíček}'s rational approximation generalized},
	issn = {0035-8711, 1365-2966},
	shorttitle = {The {Voigt} and complex error function},
	url = {https://academic.oup.com/mnras/advance-article/doi/10.1093/mnras/sty1680/5045263},
	doi = {10.1093/mnras/sty1680},
	abstract = {Accurate yet efficient computation of the Voigt and complex error function is a challenge since decades in astrophysics and other areas of physics. Rational approximations have attracted considerable attention and are used in many codes, often in combination with other techniques. The 12-term code “cpf12” of Humlíček (1979) achieves an accuracy of five to six significant digits throughout the entire complex plane. Here we generalize this algorithm to a larger (even) number of terms. The n = 16 approximation has a relative accuracy better than 10−5 for almost the entire complex plane except for very small imaginary values of the argument even without the correction term required for the cpf12 algorithm. With 20 terms the accuracy is better than 10−6. In addition to the accuracy assessment we discuss methods for optimization and propose a combination of the 16-term approximation with the asymptotic approximation of Humlíček (1982) for high efficiency.},
	language = {en},
	urldate = {2018-07-03},
	journal = {Monthly Notices of the Royal Astronomical Society},
	author = {Schreier, Franz},
	month = jun,
	year = {2018},
	keywords = {mentions sympy, spectroscopy},
}

Downloads: 0