[論文レビュー] An Incremental Algorithm for Computing Cylindrical Algebraic Decompositions

Changbo Chen, Marc Moreno Maza|arXiv (Cornell University)|Oct 19, 2012

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

ひとこと要約

本稿では、まず複素円筒木（CCD）を構築し、その後それを実CADに精錬する、增量的アルゴリズムを提示する。この手法は、GCDおよび平方自由計算を用いて、既存の木の枝を精錬する増分的「Intersect」演算を用いる。等式制約を活用することで、重複計算を回避し、多数のベンチマーク例において従来のCAD手法に比べ顕著な性能向上を達成する。

ABSTRACT

In this paper, we propose an incremental algorithm for computing cylindrical algebraic decompositions. The algorithm consists of two parts: computing a complex cylindrical tree and refining this complex tree into a cylindrical tree in real space. The incrementality comes from the first part of the algorithm, where a complex cylindrical tree is constructed by refining a previous complex cylindrical tree with a polynomial constraint. We have implemented our algorithm in Maple. The experimentation shows that the proposed algorithm outperforms existing ones for many examples taken from the literature.

研究の動機と目的

• 従来の投影・リフトスキームにおける重複計算を回避する増分的アルゴリズムを、円筒代数的分解（CAD）の計算に開発すること。
• 係数が正次元的成分上で同時に消える場合に失敗するMcCallumおよびBrownの投影作用素の限界を克服すること。
• 増分的精錬プロセスを通じて等式制約を体系的に統合し、効率を向上させること。
• 特に複雑または冗長な構造を有する系に対して、従来のCAD構築手法のより効率的な代替手段を提供すること。
• 三角分解を介して、直接に非交差的かつ円筒的に配置された分解を計算可能とし、実代数幾何や量化子除去への応用を容易にすること。

提案手法

• アルゴリズムは、新たな多項式制約を既存のCCDに適用する「Intersect」操作を用いて、複素円筒分解（CCD）を増分的に構築する。
• コアとなる操作「Intersect」は、現在の変数レベルに関して多項式の平方自由部分とGCDを計算することで、木のパスを精錬する。
• 段階的投影とパス抽出を用いて符号情報を伝搬させ、必要に応じてのみ頂点を精錬する。
• 多項式が消えるか否かに基づいて細胞を分割するために、GCDと余因子を計算することで等式制約を活用する。
• 各パスが単純系に対応する木データ構造を維持し、円筒的配置と非交差性を保証する。
• 各段階で符号不変性のチェックを統合し、精錬された細胞全体にわたって多項式の符号を更新する。

実験結果

リサーチクエスチョン

• RQ1 従来の投影・リフト手法と比較して、CAD構築の増分的アプローチは重複計算を削減できるか？
• RQ2 CAD精錬の過程で等式制約を体系的に活用することで、効率を向上させることは可能か？
• RQ3 三角分解と複素円筒木の使用は、従来の投影ベース手法と比較して、より効率的なCAD計算パイプラインを実現できるか？
• RQ4 実行時間および細胞数の観点から、増分的精錬戦略はバッチ構築に比べてどの程度優れているか？
• RQ5 増分的CCD構築を用いて、等投影的またはThomas分解のような正規表現を拡張して計算可能か？

主な発見

• 提示された増分的アルゴリズムは、文脈文献の多数のベンチマーク例において、従来のCAD手法を上回る性能を示す。特に、複雑または冗長な構造を有する系において顕著である。
• 係数が正次元的成分上で同時に消える場合を体系的に処理することで、McCallumおよびBrownの投影作用素の欠陥を回避する。
• GCDおよび平方自由計算に基づく増分的「Intersect」操作を用いることで、不要な分岐と計算を削減する。
• アルゴリズムは自然に非交差的かつ円筒的に配置された分解を計算可能であり、等投影的またはThomas分解のような正規表現の導出に利用可能である。
• 有限個の複素解を有する系に対して、アルゴリズムは臨界ペアを含まない三角分解に特化し、単に「Split & Merge」アルゴリズムの「Merge」ステップのみを必要とする。
• Mapleにおける実装は、非増分的かつブラックボックスベースの手法と比較して、増分的CCD構築が計算時間を顕著に短縮することを示している。