[论文解读] Exact Symbolic-Numeric Computation of Planar Algebraic Curves

Eric Berberich, Pavel Emeliyanenko+2 more|arXiv (Cornell University)|Jan 7, 2012

Polynomial and algebraic computation参考文献 5被引用 5

一句话总结

本文提出了一种经过认证、完整且高效的算法，用于通过结合符号与数值技术计算平面代数曲线的排列。该算法引入了两个核心组件——Bisolve 用于隔离双变量代数系统的实解，以及 GeoTop 用于分析单条曲线——仅使用结式和最大公因式作为符号运算，通过 GPU 计算加速，并避免坐标变换，从而在复杂实例上的性能优于现有方法。

ABSTRACT

We present a novel certified and complete algorithm to compute arrangements of real planar algebraic curves. It 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. From a high-level perspective, the overall method splits into two main subroutines, namely an algorithm denoted Bisolve to isolate the real solutions of a zero-dimensional bivariate system, and an algorithm denoted GeoTop to analyze a single algebraic curve. Compared to existing approaches based on elimination techniques, we considerably improve the corresponding lifting steps in both subroutines. As a result, generic position of the input system is never assumed, and thus our algorithm never demands for any change of coordinates. In addition, we significantly limit the types of involved exact operations, that is, we only use resultant and gcd computations as purely symbolic operations. The latter results are achieved by combining techniques from different fields such as (modular) symbolic computation, numerical analysis and algebraic geometry. We have implemented our algorithms as prototypical contributions to the C++-project CGAL. They exploit graphics hardware to expedite the symbolic computations. We have also compared our implementation with the current reference implementations, that is, LGP and Maple's Isolate for polynomial system solving, and CGAL's bivariate algebraic kernel for analyses and arrangement computations of algebraic curves. For various series of challenging instances, our exhaustive experiments show that the new implementations outperform the existing ones.

研究动机与目标

• 开发一种无需坐标变换即可计算实平面代数曲线排列的完整且经过认证的算法。
• 通过仅将符号计算限制在结式和最大公因式运算上，最小化符号计算量，减少对昂贵代数运算的依赖。
• 通过利用 GPU 加速算术运算，提高对几何复杂曲线（包括奇点曲线和稀疏曲线）的效率与鲁棒性。
• 即使在非通用输入下，也能实现在原始坐标系中对曲线排列进行精确且完整的几何拓扑分析。
• 提供一个实际的高性能实现，并集成到 CGAL 中，其性能优于当前的参考实现。

提出的方法

• 该算法将排列计算分解为两个主要子程序：Bisolve 用于隔离零维双变量代数系统的实解，以及 GeoTop 用于分析单条代数曲线的拓扑结构。
• Bisolve 使用结式计算将解沿 x 和 y 方向投影，然后应用一种新颖的包含判定条件，结合结式数据与同伦论证，以在不进行坐标变换的情况下认证隔离框。
• 仅使用两种符号运算：双变量多项式结式计算和单变量多项式最大公因式计算，从而最小化符号计算开销。
• 该方法利用 GPU 加速结式和最大公因式的计算，与基于 CPU 的计算相比，速度提升超过 100 倍。
• 该方法避免了任何坐标变换，使得即使在具有共享坐标的非通用输入下，也能直接在原始坐标系中进行计算。
• 通过使用过滤技术与经过认证的近似算术，减少昂贵的符号运算次数，将计算聚焦于关键步骤。

实验结果

研究问题

• RQ1 混合符号-数值方法是否能够以完全认证和完整的方式计算平面代数曲线的排列，同时最小化符号运算？
• RQ2 GPU 加速的结式与最大公因式计算在多大程度上可以提升代数曲线符号-数值算法的性能？
• RQ3 消除法在多大程度上可以避免坐标变换而不损失鲁棒性或效率？
• RQ4 所提出的方法是否能够比现有方法更高效地处理几何复杂的曲线（如奇点曲线、稀疏曲线或相交曲线）？
• RQ5 与 Lgp、Maple 的 Isolate 以及 CGAL 的双变量内核等参考实现相比，新算法在挑战性实例上的性能提升如何？

主要发现

• 所提出的算法 GeoTopAK_2 在困难实例（如投影交线和特殊曲线组合）上优于现有实现，部分基准测试中速度提升高达 10 倍。
• 对于 10 条特殊曲线组合类型的曲线，GeoTopAK_2 耗时 1201.34 秒，而 Ak_2 在 4000 秒后仍超时，显示出显著的可扩展性优势。
• 当曲线数量达到 8 条或以上时，Ak_2 超时，而 GeoTopAK_2 在 1200 秒内完成，表明其对复杂度具有显著的鲁棒性。
• GPU 加速的结式与最大公因式计算消除了符号计算的瓶颈，使其不再成为算法的性能限制因素。
• 该方法在不进行坐标变换的情况下实现了完全认证与完整性，即使在具有共享坐标的非通用输入下也成立。
• 新算法集成到网络演示中后已获得广泛使用，表明其具有实际相关性与可用性。