Efficient algorithmic implementation of the Voigt/complex error function based on exponential series approximation

Loading...
Thumbnail Image

Date

2011-11-01

Authors

Abrarov, S. M.
Quine, B. M.

Journal Title

Journal ISSN

Volume Title

Publisher

Elsevier, Applied Mathematics and Computation

Abstract

We show that a Fourier expansion of the exponential multiplier yields an exponential series that can compute high-accuracy values of the complex error function in a rapid algorithm. Numerical error analysis and computational test reveal that with essentially higher accuracy it is as fast as FFT-based Weideman’s algorithm at a regular size of the input array and considerably faster at an extended size of the input array. As this exponential series approximation is based only on elementary functions, the algorithm can be implemented utilizing freely available functions from the standard libraries of most programming languages. Due to its simplicity, rapidness, high-accuracy and coverage of the entire complex plane, the algorithm is efficient and practically convenient in numerical methods related to the spectral line broadening and other applications requiring error-function evaluation over extended input arrays.

Description

Keywords

Complex error function, Voigt function, Faddeeva function, Weideman’s algorithm, Complex probability function, Plasma dispersion function, Spectral line broadening

Citation

S. M. Abrarov and B. M. Quine, “Efficient algorithmic implementation of the Voigt/complex error function based on exponential series approximation,” Applied Mathematics and Computation, vol. 218, no. 5, pp. 1894–1902, Nov. 2011.