[論文レビュー] Arrangement Computation for Planar Algebraic Curves

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

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

ひとこと要約

本稿では、結果式とGCDのみを記号的操作として用い、認証済みの複素根ソルバを介したハイブリッド記号的数値リフトステップを組み合わせることで、平面上の代数的曲線の配置を認証済みかつ完全に計算するアルゴリズムを提示する。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.

研究の動機と目的

• 座標変換を用いずに、実平面上の代数的曲線の配置を認証済みかつ完全に計算するアルゴリズムを開発すること。
• 曲線トポロジー計算における高価な記号的操作（例えば、部分結果式やグレブナー基底）への依存度を低減すること。
• 従来の手法が計算のボトルネックに陥る高次または特異な曲線に対して、効率を向上させること。
• 結果式およびGCDの計算をGPUで加速することで、実用的な性能向上を実現すること。
• 記号的計算を最小限に抑え、歪み変換を回避しながら、幾何的正確性と完全性を維持すること。

提案手法

• アルゴリズムは、平面の円筒的代数的分解を用い、単一および対ごとの曲線を幾何的トポロジカル分解により分析することで、代数的曲線の配置を計算する。
• 記号的操作として結果式とGCDのみを用い、これらをグラフィックスハードウェアにオフロードして加速する。
• 認証済みの数値的複素根ソルバと、Teissierによる複素根数の理論的上限を組み合わせたハイブリッドリフトステップにより、実根を効率的に分離する。
• 主な根分離法として、記号的数値フィルタであるFastLiftを用い、稀な場合を除き失敗しない。
• FastLiftが失敗する場合には、[4]に基づくフォールバック法を用いて完全性を保証する。
• すべての計算を元の座標系で実行することで、入力のスパarsityを保持し、歪みによる係数の爆発を回避する。

実験結果

リサーチクエスチョン

• RQ1 代数的曲線の配置計算における記号的操作の計算コストを、正しさと完全性を維持したまま顕著に低減できるか？
• RQ2 リフトステップにおける根分離に、ハイブリッド記号的数値的手法がどの程度有効であるか？
• RQ3 GPUアクセラレーションは、代数的曲線の結果式やGCDのような記号的計算タスクにおいて、性能向上にどの程度寄与するか？
• RQ4 座標の歪みが特異的または非一般な曲線の配置において、性能と数値的安定性にどのような影響を与えるか？
• RQ5 認証済みかつ完全なアルゴリズムは、特に高次または射影された交差曲線において、現在の基準実装（CGALのAk_2）をハードなベンチマークインスタンスで上回れるか？

主な発見

• 容易なインスタンス（一般位置に配置された非特異曲線など）において、FastKernelはGPUアクセラレーテッド版Ak_2の平均2倍の速度を達成した。
• 表面の交差曲線や特異点を持つ特殊な曲線においては、5倍から120倍の高速化が得られ、曲線の次数が高くなるほど改善が顕著になった。
• 最も挑戦的なテストケース（X型の11本の曲線）では、FastKernelは2696.15秒で完了したが、Ak_2は4000秒後にタイムアウトした。
• ほとんどの場合、記号的数値フィルタであるFastLiftはすべてのファイバーで成功し、失敗は稀なケースに限られ、その際はフォールバックのLift手法が起動された。
• 性能向上は特に高次曲線において顕著であり、Ak_2の部分結果式計算が主なボトルネックとなるためである。
• 歪みを回避することで、係数のスパarsityが保持され、従来のアプローチで性能を劣化させる係数ビットサイズの爆発を防いだ。