opinion

Remark on “Algorithm 916: Computing the Faddeyeva and Voigt Functions”: Efficiency Improvements and Fortran Translation

Author:

Mofreh R. ZaghloulAuthors Info & Claims

ACM Transactions on Mathematical Software (TOMS), Volume 42, Issue 3

Article No.: 26, Pages 1 - 9

https://doi.org/10.1145/2806884

Published: 23 May 2016 Publication History

Get Access

Abstract

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 2. A better accuracy vs efficiency tradeoff 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 tradeoff, 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 Humlíček'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 four significant figures. This latter modification assures the reliability of the code in the solution of practical problems requiring numerous evaluation of the function for applications requiring low-accuracy computations (<10⁻⁴).

References

[1]

M. Abramowitz and I. A. Stegun. 1965. Handbook of Mathematical Functions. Dover, New York, NY.

Digital Library

Google Scholar

[2]

W. J. Cody. 1969. Rational chebyshev approximations for the error function. Math. Comp. 23, 107, 631--638.

Google Scholar

[3]

S. Eddins. 2010. Timeit Benchmarking Function, MATLAB Central File Exchange, http://www.mathworks.com/matlabcentral/fileexchange/.

Google Scholar

[4]

V. N. Faddeyeva and N. M. Terent’ev. 1961. Tables of Values of the Function w(z) = e^&mins;z²(1 + 2i/&sqrt;π∫&nolimits;₀^z e^t² dt) for Complex Argument, Gosud. Izdat. The.-Teor. Lit., Moscow, 195 English translation, Pergamon Press, New York, NY

Google Scholar

[5]

W. Gautschi. 1970. Efficient computation of the complex error function. SIAM J. Numer. Anal. 7, 1, 187--198

Crossref

Google Scholar

[6]

J. Humlíček. 1982. Optimized computation of the voigt and complex probability functions. J. Quant. Spectrosc. and Radiat. Transfer 27, 4, 437--444.

Crossref

Google Scholar

[7]

ISO/IEC. 2011. Information Technology--Programming Languages--Part 1, Base Language—Fortran (ISO/IEC 1539:2010 (E)). ISO/IEC Copyright Office, Geneva, Switzerland, 2011.

Google Scholar

[8]

K. L. Letchworth and D. C. Benner. 2007. Rapid and accurate calculation of the voigt function. J. Quant. Spectrosc. and Radiat. Transf. 107, 173--192.

Crossref

Google Scholar

[9]

G. P. M. Poppe and C. M. J. Wijers. 1990. Algorithm 680, evaluation of the complex error function. ACM Trans. Math. Software 16, 1, 47.

Digital Library

Google Scholar

[10]

J. A. C. Weideman. 1994. Computation of the complex error function. SIAM J. Numer. Anal. 31, 5, 1497--1518

Digital Library

Google Scholar

[11]

Wolfram Research, Inc. 2008. Mathematica, Version 7.0, Champaign, IL.

Google Scholar

[12]

M. R. Zaghloul and A. N. Ali. 2011. Algorithm 916: Computing the faddeyeva and voigt functions. ACM Trans. Math. Software 38, 2, 15, 1--22.

Digital Library

Google Scholar

[13]

M. R. Zaghloul. 2015. A simple reform for treating the loss of accuracy of humliceck's W4 algorithm near the real axis. arXiv:1505.05596v1 [astro-ph.IM].

Google Scholar

Cited By

View all

de Campos MRomero Sda Silva Lde Castro J(2024)Shape Anisotropy and Magnetic Texture Determination in Anisotropic and Isotropic Alnico MagnetsJOM10.1007/s11837-024-06586-3Online publication date: 14-May-2024
https://doi.org/10.1007/s11837-024-06586-3
Zaghloul M(2024)Efficient numerical algorithms for multi-precision and multi-accuracy calculation of the error functions and Dawson integral with complex argumentsNumerical Algorithms10.1007/s11075-023-01727-297:2(869-887)Online publication date: 1-Oct-2024
https://dl.acm.org/doi/10.1007/s11075-023-01727-2
Zaghloul M(2024)Efficient multiple-precision computation of the scaled complementary error function and the Dawson integralNumerical Algorithms10.1007/s11075-023-01608-895:3(1291-1308)Online publication date: 1-Mar-2024
https://dl.acm.org/doi/10.1007/s11075-023-01608-8
Show More Cited By

Index Terms

Remark on “Algorithm 916: Computing the Faddeyeva and Voigt Functions”: Efficiency Improvements and Fortran Translation
1. Mathematics of computing
  1. Mathematical analysis
    1. Functional analysis
      1. Approximation
    2. Numerical analysis
      1. Arbitrary-precision arithmetic
  2. Mathematical software
2. Theory of computation
  1. Design and analysis of algorithms
    1. Approximation algorithms analysis

Recommendations

Algorithm 916: Computing the Faddeyeva and Voigt Functions

We present a MATLAB function for the numerical evaluation of the Faddeyeva function w(z). The function is based on a newly developed accurate algorithm. In addition to its higher accuracy, the software provides a flexible accuracy vs efficiency trade-...
Algorithm 985: Simple, Efficient, and Relatively Accurate Approximation for the Evaluation of the Faddeyeva Function

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 Humlíček (1982) along with asymptotic expressions from ...
Remark on “Algorithm 680: Evaluation of the Complex Error Function”: Cause and Remedy for the Loss of Accuracy Near the Real Axis

In this remark, we identify the cause of the loss of accuracy in the computation of the Faddeyeva function, w(z), near the real axis when using Algorithm 680. We provide a simple correction to this problem that allows us to restore this code as one of ...

Comments

Information & Contributors

Information

Published In

cover image ACM Transactions on Mathematical Software

ACM Transactions on Mathematical Software Volume 42, Issue 3

June 2016

208 pages

ISSN:0098-3500

EISSN:1557-7295

DOI:10.1145/2935754

Editor:
Michael A. Heroux
Sandia National Laboratories, USA

Issue’s Table of Contents

Permission to make digital or hard copies of part or all of this work for personal or classroom use is granted without fee provided that copies are not made or distributed for profit or commercial advantage and that copies bear this notice and the full citation on the first page. Copyrights for third-party components of this work must be honored. For all other uses, contact the Owner/Author.

Publisher

Association for Computing Machinery

New York, NY, United States

Publication History

Published: 23 May 2016

Accepted: 01 July 2015

Revised: 01 June 2015

Received: 01 January 2015

Published in TOMS Volume 42, Issue 3

Check for updates

Author Tags

Qualifiers

Opinion
Research
Refereed

Contributors

Other Metrics

View Article Metrics

Bibliometrics & Citations

Bibliometrics

Article Metrics

7
Total Citations
View Citations
217
Total Downloads

Downloads (Last 12 months)21
Downloads (Last 6 weeks)4

Reflects downloads up to 24 Sep 2024

Other Metrics

View Author Metrics

Citations

Cited By

View all

de Campos MRomero Sda Silva Lde Castro J(2024)Shape Anisotropy and Magnetic Texture Determination in Anisotropic and Isotropic Alnico MagnetsJOM10.1007/s11837-024-06586-3Online publication date: 14-May-2024
https://doi.org/10.1007/s11837-024-06586-3
Zaghloul M(2024)Efficient numerical algorithms for multi-precision and multi-accuracy calculation of the error functions and Dawson integral with complex argumentsNumerical Algorithms10.1007/s11075-023-01727-297:2(869-887)Online publication date: 1-Oct-2024
https://dl.acm.org/doi/10.1007/s11075-023-01727-2
Zaghloul M(2024)Efficient multiple-precision computation of the scaled complementary error function and the Dawson integralNumerical Algorithms10.1007/s11075-023-01608-895:3(1291-1308)Online publication date: 1-Mar-2024
https://dl.acm.org/doi/10.1007/s11075-023-01608-8
Zaghloul MAlrawas L(2023)Calculation of Fresnel integrals of real and complex arguments up to 28 significant digitsNumerical Algorithms10.1007/s11075-023-01654-296:2(489-506)Online publication date: 14-Sep-2023
https://doi.org/10.1007/s11075-023-01654-2
Mahler AWilliams JSu NYang W(2022)Localized orbital scaling correction for periodic systemsPhysical Review B10.1103/PhysRevB.106.035147106:3Online publication date: 28-Jul-2022
https://doi.org/10.1103/PhysRevB.106.035147
Yan XCai HZhang WHao LYao PLi XZou SZhu SHe X(2019)Anomalous mix induced by a collisionless shock wave in an inertial confinement fusion hohlraumNuclear Fusion10.1088/1741-4326/ab32cf59:10(106016)Online publication date: 21-Aug-2019
https://doi.org/10.1088/1741-4326/ab32cf
Schippers S(2018)Analytical expression for the convolution of a Fano line profile with a gaussianJournal of Quantitative Spectroscopy and Radiative Transfer10.1016/j.jqsrt.2018.08.003219(33-36)Online publication date: Nov-2018
https://doi.org/10.1016/j.jqsrt.2018.08.003

View Options

Get Access

Login options

Check if you have access through your login credentials or your institution to get full access on this article.

Cited By

Index Terms

Recommendations

Algorithm 916: Computing the Faddeyeva and Voigt Functions

Algorithm 985: Simple, Efficient, and Relatively Accurate Approximation for the Evaluation of the Faddeyeva Function

Remark on “Algorithm 680: Evaluation of the Complex Error Function”: Cause and Remedy for the Loss of Accuracy Near the Real Axis