[论文解读] Software for Computing the Spheroidal Wave Functions Using Arbitrary Precision Arithmetic

Ross Adelman, Nail A. Gumerov+1 more|arXiv (Cornell University)|Aug 1, 2014

Numerical Methods and Algorithms参考文献 15被引用 20

一句话总结

本论文提出了一款名为'spheroidal'的高精度软件库，用于在C++和MATLAB中使用任意精度算术计算旋转椭球波函数（prolate and oblate spheroidal wave functions）。该库采用自适应系数选择、基于朗斯基安（Wronskian）的径向函数方法切换，以及优化的递推关系，实现了卓越的数值稳定性和精度，尤其在高频和高模态数情况下表现优异。

ABSTRACT

The spheroidal wave functions, which are the solutions to the Helmholtz equation in spheroidal coordinates, are notoriously difficult to compute. Because of this, practically no programming language comes equipped with the means to compute them. This makes problems that require their use hard to tackle. We have developed computational software for calculating these special functions. Our software is called spheroidal and includes several novel features, such as: using arbitrary precision arithmetic; adaptively choosing the number of expansion coefficients to compute and use; and using the Wronskian to choose from several different methods for computing the spheroidal radial functions to improve their accuracy. There are two types of spheroidal wave functions: the prolate kind when prolate spheroidal coordinates are used; and the oblate kind when oblate spheroidal coordinate are used. In this paper, we describe both, methods for computing them, and our software. We have made our software freely available on our webpage.

研究动机与目标

• 解决标准编程语言中缺乏可靠、高精度椭球波函数实现的问题。
• 克服双精度计算中的数值不稳定性和精度损失，特别是在高频和高模态数问题中。
• 提供一个灵活、可扩展的软件库，支持旋转椭球波函数与扁椭球波函数，并具备可定制的精度和展开系数控制能力。
• 通过使用任意精度算术和基于朗斯基安的方法选择，提升径向函数计算的精度。
• 通过持久存储中间展开系数和优化递推计算，实现高效且可复现的计算。

提出的方法

• 利用GNU MPFR库执行任意精度算术，允许用户在每个计算步骤中自定义精度级别。
• 采用两种展开系数选择模式：固定用户指定数量，或自适应计算直至系数低于用户定义的阈值（例如10⁻²⁰⁰）。
• 对增长序列（如球贝塞尔诺伊曼函数）使用前向递推，对衰减序列使用连分数法，以保持数值稳定性。
• 使用三对角矩阵法求解非齐次递推关系，特别适用于初始增长后衰减的系数（如$ B_{2r}^{mn}(-ic) $）。
• 通过递归分解优化阶乘和上升Pochhammer符号的计算，利用增量更新将每项的复杂度从O(n²)降低至O(1)。
• 通过最小化朗斯基安中的误差，动态选择最准确的算法路径，实现椭球径向函数的最优计算方法选择。

实验结果

研究问题

• RQ1 当标准双精度算术失效时，如何在高频和高波数模态下实现椭球波函数的高精度计算？
• RQ2 在递推关系先增长后衰减的情况下，应采用何种计算策略以确保数值稳定性？
• RQ3 如何有效将任意精度算术集成到特殊函数软件库中，同时将性能开销降至最低？
• RQ4 应采用何种标准自动选择计算椭球径向函数的最准确方法？
• RQ5 阶乘和Pochhammer符号项的因式分解在级数展开中能在多大程度上提升计算效率？

主要发现

• 该软件在高频和高模态数问题中显著优于双精度实现，精度大幅提升。
• 自适应系数选择确保仅计算到所需精度，减少不必要的计算，同时保持高精度。
• 通过用户自定义精度级别的任意精度算术，可在广泛物理参数范围内实现可靠计算。
• 基于朗斯基安的方法选择通过动态选择最准确的算法路径，显著降低了径向函数计算中的误差。
• 阶乘和Pochhammer符号项的递归优化将每项的计算复杂度从O(n²)降低至O(1)，极大提升了级数展开的性能。
• 该软件免费提供，并支持计算系数的持久存储，实现跨多次计算的复现与重用。