[论文解读] Arrangement Computation for Planar Algebraic Curves

Eric Berberich, Pavel Emeliyanenko+2 more|arXiv (Cornell University)|Mar 24, 2011

Advanced Numerical Analysis Techniques参考文献 26被引用 4

一句话总结

本文提出了一种认证且完备的算法，用于仅使用结式和最大公因式作为符号运算来计算平面代数曲线排列，并结合通过认证复根求解器实现的混合符号-数值提升步骤。该方法通过利用GPU加速算术运算并避免坐标剪切，相较于先前的最先进方法（CGAL的Ak_2），在困难实例上实现了高达120倍的性能提升，尤其适用于高次或奇异曲线。

ABSTRACT

We present a new certified and complete algorithm to compute arrangements of real planar algebraic curves. Our algorithm provides a geometric-topological analysis of the decomposition of the plane induced by a finite number of algebraic curves in terms of a cylindrical algebraic decomposition of the plane. Compared to previous approaches, we improve in two main aspects: Firstly, we significantly reduce the amount of exact operations, that is, our algorithms only uses resultant and gcd as purely symbolic operations. Secondly, we introduce a new hybrid method in the lifting step of our algorithm which combines the usage of a certified numerical complex root solver and information derived from the resultant computation. Additionally, we never consider any coordinate transformation and the output is also given with respect to the initial coordinate system. We implemented our algorithm as a prototypical package of the C++-library CGAL. Our implementation exploits graphics hardware to expedite the resultant and gcd computation. We also compared our implementation with the current reference implementation, that is, CGAL's curve analysis and arrangement for algebraic curves. For various series of challenging instances, our experiments show that the new implementation outperforms the existing one.

研究动机与目标

• 开发一种无需坐标变换即可计算实平面代数曲线排列的认证且完备的算法。
• 减少在曲线拓扑计算中对昂贵符号运算（如子结式和Gröbner基）的依赖。
• 提升高次曲线和奇异曲线的计算效率，这些曲线在先前方法中常因计算瓶颈而受阻。
• 通过GPU加速结式和最大公因式计算，实现实际性能提升。
• 在最小化符号计算并避免剪切变换的同时，保持几何精度和完备性。

提出的方法

• 该算法使用平面的柱状代数分解来计算代数曲线的排列，通过几何-拓扑分解分析单个和成对曲线。
• 仅使用结式和最大公因式作为符号运算，并将这些运算卸载到图形硬件以实现加速。
• 混合提升步骤结合了认证的数值复根求解器与Teissier关于复根数量的理论边界，以高效隔离实根。
• 使用一种称为FastLift的符号-数值过滤器作为主要根隔离方法，仅在极少数情况下失败。
• 当FastLift失败时，采用基于[4]的回退方法以确保完备性。
• 整个计算在原始坐标系中进行，保留了输入的稀疏性，并避免了剪切带来的系数膨胀。

实验结果

研究问题

• RQ1 是否能显著降低代数曲线排列计算中符号运算的计算成本，同时保持正确性和完备性？
• RQ2 混合符号-数值方法在曲线拓扑计算的提升步骤中的根隔离效果如何？
• RQ3 GPU加速在结式和最大公因式等符号计算任务中对代数曲线性能的提升程度如何？
• RQ4 在奇异或非一般位置的曲线排列中，避免坐标剪切对性能和数值稳定性的影响如何？
• RQ5 在高难度基准实例上，尤其是高次曲线或投影交线时，认证且完备的算法能否超越当前参考实现（CGAL的Ak_2）？

主要发现

• 对于非奇异且处于一般位置的曲线等简单实例，新实现FastKernel的平均速度是GPU加速版Ak_2的两倍。
• 对于曲面的投影交线或具有奇异性的特殊曲线，性能提升在5倍至120倍之间，且随着曲线次数增加而提升更显著。
• 在最困难的测试用例（11条X型曲线）中，FastKernel耗时2696.15秒完成，而Ak_2在4000秒内超时。
• 符号-数值过滤器FastLift在大多数情况下成功处理了所有纤维，仅在极少数情况下触发了回退的Lift方法。
• 性能提升在高次曲线中最为显著，因为Ak_2的子结式计算成为主要瓶颈。
• 该方法避免了剪切，保留了系数的稀疏性，防止了先前方法中因系数位数膨胀而导致的性能下降。