TY - JOUR
T1 - Efficient multiple-precision computation of the scaled complementary error function and the Dawson integral
AU - Zaghloul, Mofreh R.
N1 - Publisher Copyright:
© The Author(s) 2023.
PY - 2024/3
Y1 - 2024/3
N2 - We present algorithms to approximate the scaled complementary error function, expx2erfc(x), and the Dawson integral, e-x2∫x0et2dt, to the best accuracy in the standard single, double, and quadruple precision arithmetic. The algorithms are based on expansion in Chebyshev subinterval polynomial approximations together with expansion in terms of Taylor series and/or Laplace continued fraction. The present algorithms, implemented as Fortran elemental modules, have been benchmarked versus competitive algorithms available in the literature and versus functions built-in in modern Fortran compilers, in addition to comprehensive tables generated with variable precision computations using the Matlab™ symbolic toolbox. The present algorithm for calculating the scaled complementary error function showed an overall significant efficiency improvement (factors between 1.3 and 20 depending on the compiler and tested dataset) compared to the built-in function “Erfc_Scaled” in modern Fortran compilers, whereas the algorithm for calculating the Dawson integral is exceptional in calculating the function to 32 significant digits (compared to 19 significant digits reported in the literature) while being more efficient than competitive algorithms as well.
AB - We present algorithms to approximate the scaled complementary error function, expx2erfc(x), and the Dawson integral, e-x2∫x0et2dt, to the best accuracy in the standard single, double, and quadruple precision arithmetic. The algorithms are based on expansion in Chebyshev subinterval polynomial approximations together with expansion in terms of Taylor series and/or Laplace continued fraction. The present algorithms, implemented as Fortran elemental modules, have been benchmarked versus competitive algorithms available in the literature and versus functions built-in in modern Fortran compilers, in addition to comprehensive tables generated with variable precision computations using the Matlab™ symbolic toolbox. The present algorithm for calculating the scaled complementary error function showed an overall significant efficiency improvement (factors between 1.3 and 20 depending on the compiler and tested dataset) compared to the built-in function “Erfc_Scaled” in modern Fortran compilers, whereas the algorithm for calculating the Dawson integral is exceptional in calculating the function to 32 significant digits (compared to 19 significant digits reported in the literature) while being more efficient than competitive algorithms as well.
KW - Dawson integral
KW - Fortran
KW - Scaled complementary error function
KW - Special functions
UR - http://www.scopus.com/inward/record.url?scp=85167361905&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85167361905&partnerID=8YFLogxK
U2 - 10.1007/s11075-023-01608-8
DO - 10.1007/s11075-023-01608-8
M3 - Article
AN - SCOPUS:85167361905
SN - 1017-1398
VL - 95
SP - 1291
EP - 1308
JO - Numerical Algorithms
JF - Numerical Algorithms
IS - 3
ER -