[论文解读] Computations modulo regular chains

Xin Li, Marc Moreno Maza+1 more|ArXiv.org|Mar 21, 2009

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

一句话总结

本文提出了一种基于模方法和快速多项式算术的新算法，用于计算正则链模下的正则结式最大公因式（GCD）与正则性测试，显著加速了多项式系统的三角分解。该方法利用结式理论与基于FFT的求值/插值技术，相较于现有的Maple和Magma实现，性能提升达数个数量级。

ABSTRACT

The computation of triangular decompositions are based on two fundamental operations: polynomial GCDs modulo regular chains and regularity test modulo saturated ideals. We propose new algorithms for these core operations relying on modular methods and fast polynomial arithmetic. Our strategies take also advantage of the context in which these operations are performed. We report on extensive experimentation, comparing our code to pre-existing Maple implementations, as well as more optimized Magma functions. In most cases, our new code outperforms the other packages by several orders of magnitude.

研究动机与目标

• 为解决多项式系统三角分解中的性能瓶颈，特别是正则链模下正则GCD与正则性测试的计算问题。
• 实现在非整环上无需因式分解的正则GCD高效计算，尤其适用于正维正则链。
• 通过在基域上而非模饱和理想上计算，减少表达式膨胀并避免重复中间计算。
• 将快速多项式算术与模方法集成到三角分解的核心运算中。
• 提供可扩展、高性能的实现，适用于大规模密集多项式系统。

提出的方法

• 该方法以结式链为基础计算正则GCD，并建立了子结式作为正则GCD的理论条件。
• 采用基于FFT的求值/插值技术，高效地在模算术中计算结式链。
• 仅计算一次结式链，重复用于正则GCD与结式计算，避免冗余。
• 使用模方法在有限域上执行算术运算，并通过提升至基域控制表达式膨胀。
• 通过在基域上的多项式环中工作，避免直接计算模饱和理想，降低开销。
• 采用自底向上的策略替代传统的自顶向下PRS方法，更好地集成快速算术与模技术。

实验结果

研究问题

• RQ1 当饱和理想非根时，子结式多项式在何种条件下是正则链模下的正则GCD？
• RQ2 如何利用快速多项式算术与模方法优化正则GCD计算？
• RQ3 子结式链能否被重用于同时计算正则GCD与结式，以减少冗余计算？
• RQ4 所提方法在大规模密集多项式系统上与Maple和Magma现有实现相比性能如何？
• RQ5 中间分裂对性能有何影响？新算法如何处理此类情况？

主要发现

• 新实现的正则性测试与正则GCD计算性能相比Maple的RegularChains库提升了数个数量级。
• 对于2变量随机密集系统，快速正则性测试速度比Maple中原始的Regularize命令快达3000倍。
• 在含多个分裂的3变量密集系统中，新算法仍显著快于通用Maple实现，尽管由于接口开销略慢于Magma的优化代码。
• 快速正则性测试在输入次数达30时可实现亚秒级性能，而原始Maple命令对类似输入耗时超过1000秒。
• 基于FFT的求值/插值技术降低了计算成本与表达式膨胀，使大规模系统得以高效处理。
• 该方法具有良好的可扩展性：对于3变量系统，次数分别为12、22、33时，新算法在300秒内完成，而旧版Maple版本则超时。