[論文レビュー] Divide-And-Conquer Computation of Cylindrical Algebraic Decomposition

Adam Strzeboński|arXiv (Cornell University)|Feb 4, 2014

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

ひとこと要約

本論文は、複雑な多項式系をより小さな部分式に分割し、それぞれでCADを計算した後、ブール結合アルゴリズムを用いて統合する、分割統治型の円柱代数的分解（CAD）の手法を提示する。特に、グラフに基づくヒューリスティクスとパラメータ $ p = 0.75 $ を用いることで、直接CADと比較して平均で最大3.96倍の高速化を達成し、項間の多項式被覆が少ない論理和系に対して特に有効である。

ABSTRACT

We present a divide-and-conquer version of the Cylindrical Algebraic Decomposition (CAD) algorithm. The algorithm represents the input as a Boolean combination of subformulas, computes cylindrical algebraic decompositions of solution sets of the subformulas, and combines the results. We propose a graph-based heuristic to find a suitable partitioning of the input and present empirical comparison with direct CAD computation.

研究の動機と目的

• 大規模または複雑な実多項式系に対する円柱代数的分解（CAD）の効率を向上させること。
• 多数の論理和項を含む系において、直接CADの高い計算コストを軽減すること。
• 共有多項式を最小限に抑えた部分式への分解を通じて、スケーラブルでモジュラーなCADアプローチを構築すること。
• グラフベースの分割とパrameter化されたヒューリスティクスを用いた分割統治型CADの性能を評価すること。

提案手法

• アルゴリズムは、実多項式系 $ S $ をブール結合 $ \bigvee_{i=1}^m P_i $ に分解する。ここで各 $ P_i $ は部分式を表す。
• グラフベースのヒューリスティクスは、論理和項を頂点とし、共有多項式に基づいてエッジを張るグラフを構築する。重み閾値パラメータ $ p $ を用いる。
• グラフ内の連結成分が部分式 $ P_1, \dots, P_m $ を定義し、それらの間の多項式被覆を最小化する。
• 各 $ P_i $ は標準CADを用いて処理され、円柱代数的式（CAF）が得られる。
• 結果はCAFCombineアルゴリズムを用いて統合され、元の $ S $ と等価なCAFが再構築される。
• パラメータ $ p $ は分割の感度を制御する：$ p $ が大きいほど、粗い、より不連続な分割が得られる。

実験結果

リサーチクエスチョン

• RQ1 分割統治戦略は、大規模な論理和系におけるCADの計算コストを低減できるか？
• RQ2 分割ヒューリスティクスの選択が、CAD計算の性能にどのように影響するか？
• RQ3 グラフベースの分割アルゴリズムにおけるパラメータ $ p $ のどの値が多様な入力に対して最良の性能を発揮するか？
• RQ4 どのような状況で分割統治型CADが直接CADを上回り、その差はどの程度か？

主な発見

• 12個のベンチマーク例（[5]より）において、$ p = 0.75 $ の分割統治型CADは、直接CADと比較して平均2.25倍の高速化を達成した。
• 14個のランダム例において、$ p = 0.75 $ の分割統治型CADは、直接CADと比較して最低でも3.96倍速く、2例では直接CADが1200秒以内に完了しなかった。
• [5]の12例中9例において、$ p = 0.75 $ の分割統治型CADは直接CADより速かった。残りの3例では遅かった。
• 入力が論理和であり、部分式同士が少ない多項式を共有する場合に、この手法は特に有効であり、グラフベースの分割ヒューリスティクスによって裏付けられた。
• ベンチマーク例およびランダム例の両方において、$ p = 0.75 $ は $ p = 0.25, 0.5, 1 $ と比較して平均性能で常に優れていた。
• 28例中10例では、$ p = 1 $ ですら直接CADより遅くなった。これは、あまりに細かすぎる分割が有害である可能性を示している。