A sampling-based approximation of the complex error function and its implementation without poles

Abstract

Recently we developed a new sampling methodology based on incomplete cosine expansion of the sinc function and applied it in numerical integration in order to obtain a rational approximation for the complex error function $w\left(z\right)=e^{-z^2}\left(1+\frac{2i}{\sqrt{\pi }}{\int }_0^z{e}^{t^2}dt\right),$ where $z=x+iy$. As a further development, in this work we show how this sampling-based rational approximation can be transformed into alternative form for efficient computation of the complex error function $w\left(z\right)$ at smaller values of the imaginary argument $y=\mathrm{Im}\left[z\right]$. Such an approach enables us to avoid poles in implementation and to cover the entire complex plain with high accuracy in a rapid algorithm. An optimized Matlab code utilizing only three rapid approximations is presented.

arXiv comment20 pages

References in corpus (7)

• Decay time integrals in neutral meson mixing and their efficient evaluation
• A rational approximation of the Dawson's integral for efficient computation of the complex error function
• A rational approximation for efficient computation of the Voigt function in quantitative spectroscopy
• Sampling by incomplete cosine expansion of the sinc function: application to the Voigt/complex error function
• Accurate approximations for the complex error function with small imaginary argument
• Analytical evaluation and asymptotic evaluation of Dawson's integral and related functions in mathematical physics
• The Fourier expansion approximation for high-accuracy computation of the Voigt/complex error function at small imaginary argument

Cited by in corpus (7)

• Efficient application of the Radiance Enhancement method for detection of the forest fires due to combustion-originated reflectance
• A single-domain implementation of the Voigt/complex error function by vectorized interpolation
• A two-domain MATLAB implementation for efficient computation of the Voigt/complex error function
• Effect of the instrument slit function on upwelling radiance from a wavelength dependent surface reflectance
• Analytical valuation of some non-elementary integrals involving some exponential, hyperbolic and trigonometric elementary functions and derivation of new probability measures generalizing the gamma-type and normal distributions
• Rapid computation of the total band radiance by using the Spectrally Integrated Voigt Function
• Efficient application of the Voigt functions in the Fourier transform