[论文解读] Multiprecision path tracking

Daniel J. Bates, Andrew J. Sommese+1 more|ArXiv.org|May 3, 2006

Polynomial and algebraic computation参考文献 15被引用 4

一句话总结

本文提出了一种用于数值代数几何中同伦延续的自适应精度路径跟踪算法，根据多项式系统的条件数动态调整计算精度。通过将精度控制与自适应步长相结合，该方法在病态或高次多项式系统中表现出高可靠性和高效性，尤其在度数高达300的切比雪夫多项式上实现了完整精度。

ABSTRACT

A path tracking algorithm that adaptively adjusts precision is presented. By adjusting the level of precision in accordance with the numerical conditioning of the path, the algorithm achieves high reliability with less computational cost than would be incurred by raising precision across the board. We develop simple rules for adjusting precision and show how to integrate these into an algorithm that also adaptively adjusts the step size. The behavior of the method is illustrated on several examples arising as homotopies for solving systems of polynomial equations.

研究动机与目标

• 为解决病态或高次多项式系统中固定精度路径跟踪的低效与不可靠问题。
• 通过根据数值条件性按路径自适应调整精度，而非对所有路径统一使用高精度算术，从而降低计算成本。
• 在奇点附近或总次数较高的系统中实现可靠解路径跟踪，此类情况标准双精度会失效。
• 将精度自适应与自适应步长控制相结合，实现鲁棒且高效的路径跟踪。
• 在具有严格精度要求的高难度问题（如高次切比雪夫多项式）上验证该方法的有效性。

提出的方法

• 该算法根据同伦函数沿路径的条件数，自适应调整浮点运算中使用的十进制位数。
• 将精度自适应与自适应步长控制相结合，其中步长根据校正器的收敛行为进行调整。
• 采用预测-校正框架，使用欧拉法或更高阶预测器及基于牛顿的校正器，并施加约束以隔离正确解路径。
• 通过雅可比矩阵和残差范数估计条件数，以指导精度的提升。
• 通过增加额外数位设置安全余量，以防止校正过程中失败，参数通过经验测试进行调优。
• 当检测到精度不足时，算法会逐步增加精度，直至实现收敛。

实验结果

研究问题

• RQ1 在路径跟踪过程中，如何自适应调整精度，以在不引入统一高精度算术成本的前提下提升可靠性？
• RQ2 标准双精度在路径跟踪中失效的条件是什么？如何高效地响应以提升精度？
• RQ3 自适应精度控制能否与自适应步长有效结合，以同时增强鲁棒性与效率？
• RQ4 该方法在高次多项式系统（尤其是根聚集或病态系统）上的表现如何？
• RQ5 自适应精度在多大程度上可以替代或补充奇异终点的末端方法？

主要发现

• 自适应精度算法成功求解了度数为300的切比雪夫多项式，所有根均达到至少10位有效数字的精度，且精度水平依度数和路径而异。
• 对于度数300的系统，最病态的路径在x=1.0附近需要高达约3.4×10^24比特的精度，显示出对根聚集的极端敏感性。
• 所需精度最高的路径始终是终点靠近x=1.0的那条，因为该区域根间距最小。
• 为避免校正循环失败，必须采用较小的初始步长；若步长过大，则会导致精度需求过高并引发非收敛。
• 该方法在所有测试度数（10至300）上均实现了完整精度，且同一度数下不同路径的精度水平存在显著差异。
• 与固定精度方法相比，该算法通过避免在条件良好段落进行不必要的高精度计算，同时确保病态段落的可靠性，从而表现出更优性能。