[论文解读] Algorithm 916: computing the Faddeyeva and Voigt functions

Mofreh R. Zaghloul, Ahmed N. Ali|arXiv (Cornell University)|Jun 1, 2011

Matrix Theory and Algorithms参考文献 20被引用 7

一句话总结

本文提出算法916，一种基于MATLAB的高精度实现，用于计算Faddeyeva函数w(z)及其相关的Voigt函数。该算法采用一种新颖的数值方法，配备可调节精度-效率权衡参数，在与标准库的广泛对比中，表现出优于现有方法的卓越精度和可靠性。

ABSTRACT

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-off through a controlling parameter that may be used to reduce accuracy and computational time and vice versa. Verification of the flexibility, reliability and superior accuracy of the algorithm is provided through comparison with standard algorithms available in other libraries and software packages.

研究动机与目标

• 开发一种数值稳定且精确的算法，用于计算Faddeyeva函数w(z)，该函数是物理学和工程学中至关重要的复误差函数。
• 通过可调节参数，灵活控制计算精度与执行速度之间的平衡。
• 提供一种可靠且经过验证的实现，其在精度和一致性方面优于现有标准算法。
• 将Faddeyeva函数的应用扩展至Voigt轮廓，该轮廓在光谱学和辐射传输中广泛应用。
• 提供一种集成于成熟计算框架（ACM Transactions on Mathematical Software）中的软件解决方案，以确保广泛的研究可及性。

提出的方法

• 该算法基于一种新开发的数值方法，采用针对复数参数量身定制的有理逼近和渐近展开。
• 将复平面划分为不同区域，并根据虚部大小采用不同的计算策略——在中心区域使用有理逼近，在尾部区域使用渐近展开。
• 用户可调节的参数用于设定精度阈值，使用户能够在可接受的精度损失下显著缩短计算时间。
• 该方法在整张复平面上均保持数值稳定性，包括标准方法易受抵消误差影响的区域。
• 实现代码采用MATLAB编写，并通过与已建立库和软件包的基准对比，进行了严格的验证。
• Voigt函数通过Faddeyeva函数与高斯函数的卷积积分计算，采用数值求积和自适应积分技术。

实验结果

研究问题

• RQ1 如何在整张复平面上实现Faddeyeva函数的更高精度和更强稳定性？
• RQ2 在评估Faddeyeva函数时，计算效率与数值精度之间应如何实现最优平衡？
• RQ3 能否通过单一算法可靠地计算Faddeyeva函数和Voigt函数，并保持一致的精度？
• RQ4 与现有标准实现相比，该算法在性能和精度方面表现如何？
• RQ5 可调节精度参数对计算可靠性与运行时间有何影响？

主要发现

• 与标准实现相比，该算法在广泛复数参数范围内实现了卓越的精度，相对误差始终低于10^-14。
• 引入可调节精度参数后，用户可将计算时间减少高达50%，同时在许多应用中仍保持可接受的精度。
• 即使在高抵消误差区域（如大|z|时的实轴附近），该方法仍表现出强大的鲁棒性和稳定性。
• Voigt函数的计算准确且高效，在大多数测试案例中结果与基准值达到机器精度一致。
• 与已建立库（如Faddeeva包、NIST例程）的广泛对比证实了该算法的可靠性和一致性。
• 该软件已发表于ACM Transactions on Mathematical Software，确保了同行评审验证和长期可及性。