[论文解读] 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，以计算柱代数分解（CAD）。该方法使用一种增量的“Intersect”操作，通过计算最大公因式（GCD）和平方自由部分，对现有树分支进行细化，利用等式约束避免冗余计算，从而在许多基准示例中显著提升了现有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构造方法更高效的替代方案。
• 通过三角分解实现直接计算不相交、柱状排列的分解，促进实代数几何与量词消除中的应用。

提出的方法

• 该算法通过使用‘Intersect’操作，将新多项式约束与先前的CCD进行细化，增量式地构建复柱分解（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构造显著减少了计算时间。