[论文解读] General Complex Polynomial Root Solver and Its Further Optimization for Binary Microlenses

J. Skowron, Andrew Gould|arXiv (Cornell University)|Mar 5, 2012

Astronomy and Astrophysical Research参考文献 2被引用 24

一句话总结

本文提出了一种新型复多项式根求解器，根据判别式动态选择拉盖尔法、牛顿法以及一种新颖的中间方法，相较于商业求解器ZROOTS，速度提升1.6–3倍。该求解器专为双星微引力透镜优化，通过问题特异性启发式方法提升根求解效率，并确保数值鲁棒性。

ABSTRACT

We present a new algorithm to solve polynomial equations, and publish its code, which is 1.6-3 times faster than the ZROOTS subroutine that is commercially available from Numerical Recipes, depending on application. The largest improvement, when compared to naive solvers, comes from a fail-safe procedure that permits us to skip the majority of the calculations in the great majority of cases, without risking catastrophic failure in the few cases that these are actually required. Second, we identify a discriminant that enables a rational choice between Laguerre's Method and Newton's Method (or a new intermediate method) on a case-by-case basis. We briefly review the history of root solving and demonstrate that "Newton's Method" was discovered neither by Newton (1671) nor by Raphson (1690), but only by Simpson (1740). Some of the arguments leading to this conclusion were first given by the British historian of science Nick Kollerstrom in 1992, but these do not appear to have penetrated the astronomical community. Finally, we argue that Numerical Recipes should voluntarily surrender its copyright protection for non-profit applications, despite the fact that, in this particular case, such protection was the major stimulant for developing our improved algorithm.

研究动机与目标

• 开发一种更快、更鲁棒的复多项式根求解器，用于天文学应用，特别是双星微引力透镜。
• 降低建模引力透镜光曲线时的计算成本，其中根求解过程消耗高达80%的处理时间。
• 通过智能选择拉盖尔法、牛顿法和一种新中间根求解方法，提升数值稳定性和效率。
• 专门针对双透镜方程中出现的五次复多项式进行优化。
• 通过发布采用宽松许可的开源代码，解决现有商业求解器（如《数值计算手册》）在可访问性和透明度方面的不足，推动开放科学。

提出的方法

• 采用容错机制，在大多数情况下跳过大部分计算，同时在边界情况下仍能保证不失败。
• 引入基于判别式的决策规则，按根为单位选择最优根求解方法（拉盖尔法、中间法或牛顿法）。
• 采用两阶段方法：首先在每个根被定位后通过连续除法进行鲁棒根求解；其次使用完整多项式对根进行精炼。
• 实现一种新型动态方法（CMPLX_LAGUERRE2NEWTON），根据收敛行为在拉盖尔法、一种新型中间方法（SG）和牛顿法之间动态切换。
• 针对双星微引力透镜实施问题特异性优化，例如按根间距排序，并对二次和三次子问题使用解析求解器。
• 采用基于收敛阈值和迭代次数限制的停止准则，以在速度与精度之间取得平衡。

实验结果

研究问题

• RQ1 能否在不牺牲数值鲁棒性的前提下，显著加速通用复多项式根求解器？
• RQ2 何种判别式可可靠地指导拉盖尔法、牛顿法与中间方法之间的选择，以实现最佳性能？
• RQ3 如何利用双透镜方程的结构特征，进一步优化根求解效率，超越通用多项式求解器的性能？
• RQ4 透镜方程求解器的数值精度实际极限是什么？这些极限在天文学界是否已知？
• RQ5 商业数值库（如《数值计算手册》）是否应自愿放弃非营利学术用途的版权，以促进开放科学？

主要发现

• 新求解器在不同应用场景下比商用ZROOTS子程序快1.6–3倍。
• 基于判别式的动态方法选择显著提升了性能，避免了次优的根求解策略。
• 容错机制使大多数情况下可跳过大部分计算，从而在不损害可靠性的情况下实现速度提升。
• 通过利用五次复多项式在双星透镜中的特定结构，该算法在双星微引力透镜中实现了显著的性能提升。
• 作者识别出透镜求解器中存在数值精度极限，这些极限在天文学界可能尚不为人知，但在数值数学领域可能已有认知。
• 开源代码已以宽松许可证发布，鼓励学术使用与透明化，作者也倡导在数值软件中推行类似的版权豁免政策。