[论文解读] Polynomial root clustering and explicit deflation

Rémi Imbach, Victor Y. Pan|arXiv (Cornell University)|Jun 11, 2019

Advanced Numerical Analysis Techniques参考文献 7被引用 5

一句话总结

本文提出了一种显式降次技术，用于复平面上多项式根的聚类，通过在隔离根簇后降低多项式次数，提升了基于划分的算法效率。该方法在高次多项式（例如，512次）上实现了最高达35%的计算加速，通过早期降次根，实测在伯努利、曼德布罗特和螺旋多项式上速度提升了1.3–1.7倍。

ABSTRACT

We seek complex roots of a univariate polynomial $P$ with real or complex coefficients. We address this problem based on recent algorithms that use subdivision and have a nearly optimal complexity. They are particularly efficient when only roots in a given Region Of Interest (ROI) are sought. We propose two improvements for root finders. The first one is applied to polynomials having only real coefficients; their roots are either real or appear in complex conjugate pairs. We show how to adapt the subdivision scheme to focus the computational effort on the imaginary positive part of the ROI. In our second improvement we deflate $P$ to decrease its degree and the arithmetic cost of the subdivision.

研究动机与目标

• 解决当仅需关注感兴趣区域（ROI）中根的子集时，全局根求解方法效率低下的问题。
• 提升基于划分的局部根聚类算法在高次多项式且ROI中根较少情况下的实际效率。
• 开发并评估一种降次策略，通过移除已聚类的根来降低多项式次数和后续迭代的计算成本。
• 证明显式降次在伯努利、曼德布罗特和螺旋多项式等实际多项式族中可带来可测量的性能提升。

提出的方法

• 使用基于划分的根聚类算法，在给定ROI中隔离根的自然簇，返回带有重数信息的不相交隔离圆盘。
• 在识别出具有重数的根簇后，利用根的幂和与多项式重构技术，显式提取对应的多项式因子。
• 迭代执行降次过程：计算已隔离根的幂和，重构降次后的多项式，并在低次因子上递归处理。
• 将降次机制集成到Ccluster算法的修改版本（命名为CclusterD）中，使用基于预言机的接口处理随位精度逐步提高的多项式系数。
• 在划分树中采用深度优先搜索策略，逐步细化簇，并在每一步计算高精度的幂和近似值。
• 使用PSToCoeffs和CoeffsToPS例程在幂和与多项式系数之间转换，控制位精度，确保数值稳定性。

实验结果

研究问题

• RQ1 在根聚类后显式降次多项式因子，是否能降低复根求解中后续划分步骤的计算成本？
• RQ2 在高次多项式且根集中于局部区域的情况下，增强降次的划分算法与未降次的对应算法相比，性能如何？
• RQ3 在根聚类处理流程中，降次的最佳频率和时机为何，才能在不引入数值不稳定性的情况下最大化加速效果？
• RQ4 对于具有聚类或对称根分布的伯努利、曼德布罗特和螺旋多项式等实际多项式族，降次在多大程度上改善了运行时间？

主要发现

• 在256次伯努利多项式上，增强降次的CclusterD算法相比Ccluster实现了1.35倍的中位数加速，运行时间减少25%。
• 在512次伯努利多项式上，CclusterD耗时144秒，而Ccluster为192秒，性能提升25%。
• 在511次曼德布罗特多项式上，CclusterD耗时149秒，Ccluster为212秒，加速30%。
• 在512次螺旋多项式上，CclusterD实现23%的加速（423秒 vs. 560秒），其中降次开销最高，达85.1秒（主要来自幂和计算）。
• 在OracleForQ例程（降次计算）中累计花费的时间相对较低，占总时间的0.07–0.20，表明降次计算开销轻量。
• 隔离圆盘所需的最大位精度达到3392位（512次螺旋多项式），反映出对聚类根的高精度要求。