[论文解读] An implementation of CAD in Maple utilising McCallum projection

Matthew England|arXiv (Cornell University)|Feb 26, 2013

Advanced Numerical Analysis Techniques参考文献 23被引用 6

一句话总结

本文在 Maple 中完整实现了 McCallum 的 CADW 算法，采用其改进的投影算子与界定多项式，生成顺序不变的柱状代数分解。这是唯一完整的 CADW 实现，能够解决其他工具（如 Qepcad）因零化问题而失败的问题，并支持特定应用所需的最小化 CAD。

ABSTRACT

Cylindrical algebraic decomposition (CAD) is an important tool for the investigation of semi-algebraic sets. Originally introduced by Collins in the 1970s for use in quantifier elimination it has since found numerous applications within algebraic geometry and beyond. Following from his original work in 1988, McCallum presented an improved algorithm, CADW, which offered a huge increase in the practical utility of CAD. In 2009 a team based at the University of Western Ontario presented a new and quite separate algorithm for CAD, which was implemented and included in the computer algebra system Maple. As part of a wider project at Bath investigating CAD and its applications, Collins and McCallum's CAD algorithms have been implemented in Maple. This report details these implementations and compares them to Qepcad and the Ontario algorithm. The implementations were originally undertaken to facilitate research into the connections between the algorithms. However, the ability of the code to guarantee order-invariant output has led to its use in new research on CADs which are minimal for certain problems. In addition, the implementation described here is of interest as the only full implementation of CADW, (since Qepcad does not currently make use of McCallum's delineating polynomials), and hence can solve problems not admissible to other CAD implementations.

研究动机与目标

• 在 Maple 中实现 McCallum 的 CADW 算法，提供其改进投影算子与界定多项式的完整、可运行版本。
• 支持对不同 CAD 计算方法（特别是投影-提升法与正规链方法）之间结构与算法联系的研究。
• 支持需要顺序不变 CAD 的新应用，例如特定逻辑问题的最小化分解。
• 解决现有工具（如 Qepcad）的局限性，后者未实现 McCallum 的界定多项式，且在某些输入下可能失败。
• 在 Maple 中提供一个可靠、可扩展的 CAD 框架，适用于研究与实际应用，支持自定义变量排序与输出控制。

提出的方法

• 将 McCallum 1988 年的 CADW 算法适配为 Maple 实现，使用基于可界定性理论与最小界定多项式的相同投影算子。
• 集成 Brown 的最小界定多项式，以确保最终 CAD 输出的顺序不变性，这是高级应用的关键特性。
• 在 Maple 中实现完整的投影-提升流程，支持变量排序、单元构建以及每个单元上的堆栈形成。
• 允许用户控制输出类型，包括在输入多项式未严格要求时仍强制实现顺序不变性。
• 在一组基准问题上，将本实现（PCAD）与 Qepcad（QCAD）以及基于 RegularChains 的 CAD（TCAD）进行比较，以评估性能与正确性。
• 在提升过程中处理零化情况，允许计算继续进行，避免在安全时过早失败。

实验结果

研究问题

• RQ1 McCallum 的 CADW 算法是否能在 Maple 中完整且正确地实现，包括使用界定多项式以确保顺序不变性？
• RQ2 是否存在 Qepcad 因零化问题而失败，但新的 Maple 实现 CADW 仍能成功求解的问题？
• RQ3 该实现的顺序不变输出在多大程度上支持新应用，例如特定逻辑问题的最小化 CAD？
• RQ4 与 Qepcad 和基于 RegularChains 的 CAD 相比，新的基于 Maple 的 CADW 实现的性能（以单元数量与计算时间衡量）如何？
• RQ5 该实现是否可扩展以集成进一步改进（如 Brown-McCallum 算子），以减少单元数量并提高效率？

主要发现

• 该实现是 McCallum CADW 算法唯一完整且可运行的实现，因为 Qepcad 未实现其界定多项式。
• 在 Whitney 雨伞与 x 轴椭圆问题中，该实现成功计算了 CAD，而 Qepcad 因零化问题失败，证明了其鲁棒性。
• 在 x 轴椭圆问题中，Maple 实现生成了 20,225 个单元，与 Qepcad 输出一致，但 Qepcad 因零化报告失败。
• 当强制实现顺序不变性时，该实现正确计算出 $f = zy - x^2$ 的 23 个单元，而符号不变模式下为 21 个单元。
• 平均而言，Qepcad 在速度与单元数量上优于 Maple 实现，但 Maple 版本在 Qepcad 失败的问题中不可或缺。
• 由于保证了顺序不变性并正确处理了界定多项式，该实现支持为特定问题生成最小化 CAD，符合近期研究需求。