[論文レビュー] Truncated Normal Forms for Solving Polynomial Systems: Generalized and Efficient Algorithms

Bernard Mourrain, Simon Telen+1 more|arXiv (Cornell University)|Mar 21, 2018

Polynomial and algebraic computation参考文献 34被引用数 12

ひとこと要約

本稿では、特にチェビシェフ多項式のような直交多項式を含む非単項式基底を導入することで、零次元多項式系を解くための切り捨て正規形（TNF）法を一般化・高速化する。これにより、数値的安定性と精度が向上し、特に解が有界領域に集中する系において顕著な性能向上が達成される。SVDに基づく基底選択を用いて直交TNFを構築し、高次においても、単項式基底に比べて最大10⁷倍の残差精度向上を実現した実験結果が得られている。

ABSTRACT

We consider the problem of finding the isolated common roots of a set of polynomial functions defining a zero-dimensional ideal I in a ring R of polynomials over C. Normal form algorithms provide an algebraic approach to solve this problem. The framework presented in Telen et al. (2018) uses truncated normal forms (TNFs) to compute the algebra structure of R/I and the solutions of I. This framework allows for the use of much more general bases than the standard monomials for R/I. This is exploited in this paper to introduce the use of two special (nonmonomial) types of basis functions with nice properties. This allows, for instance, to adapt the basis functions to the expected location of the roots of I. We also propose algorithms for efficient computation of TNFs and a generalization of the construction of TNFs in the case of non-generic zero-dimensional systems. The potential of the TNF method and usefulness of the new results are exposed by many experiments.

研究の動機と目的

• 標準の単項式基底に制限されない、非一般的な零次元多項式系へも適用可能な切り捨て正規形（TNF）法の拡張。
• 単項式基底の代わりにチェビシェフ多項式のような直交多項式基底を用いることで、多項式系の解法における数値的安定性と精度の向上。
• 効率的なカーネル計算アルゴリズムを用いてTNF計算の計算量を低減。
• SVDによる基底選択が、特に実数解や凝集した解を持つ系において、根の位置をより適切に反映する直交射影をもたらすことを示すこと。
• 高次系において、非単項式基底が、標準の単項式TNFに比べて著しく改善された後向き誤差をもたらすことを示すこと。

提案手法

• 単項式に制限されない任意の基底を許容する一般化されたTNFフレームワークを提案。これにより、商代数 R/I の基底として任意の基底を用いることが可能となる。
• QRに代えてSVDを用いた基底選択を導入し、直交TNFを構築。これにより、R/I への直交射影が得られる。
• 計算コストを低減するため、修正されたマカウレイ型行列構築法を用いてカーネル写像を効率的に計算。
• 根が [−1, 1]² のような有界領域に集中すると予想される系に対して、積のチェビシェフ基底を用い、数値的条件を改善。
• 大規模な構造的行列に対して列ピボット付きQRおよびSVDを適用し、R/I の数値的に安定した基底を選択。
• 結果写像の余核写像を高速に計算するアルゴリズムを開発し、TNF計算パイプラインの高速化を実現。

実験結果

リサーチクエスチョン

• RQ1 チェビシェフ多項式のような非単項式基底（特に直交多項式）は、多項式系のTNFに基づくソルバーの精度と安定性を向上させ得るか？
• RQ2 TNF構築において、SVDに基づく基底選択はQRに比べて数値的品質と残差誤差の面で優れているか？
• RQ3 高次系において、チェビシェフ多項式を基底関数として用いることで、後向き誤差と解の精度にどのような影響を与えるか？
• RQ4 単項式順序付けに依存せずに、TNF法を非一般的な零次元系へ効果的に一般化できるか？
• RQ5 実際の応用において、基底の選択がTNFソルバーの計算量とパフォーマンスに与える影響は何か？

主な発見

• チェビシェフ多項式を基底関数として用いることで、次数 d = 25 において最大後向き誤差が 10⁻¹ から 10⁻⁶ に低下し、精度が10⁵倍向上した。
• SVDに基づくTNFは残差 5.16 × 10⁻¹² を達成したのに対し、QRでは 2.84 × 10⁻¹¹ にとどまり、数値的安定性に優れたことが示された。
• d = 20 の場合、チェビシェフTNFは単項式TNFに比べて残差が1桁分改善され、古典的なマカウレイ結果式法に比べて10⁷倍の改善が得られた。
• n = 4, d = 4 の場合、DBDアルゴリズムによるカーネル計算はPHCpackと同等の速度を達成し、Bertiniに比べて1.6倍速かった。
• n = 7, d = 2 の場合、ホモトピー法ソルバーがTNFソルバーを上回ったが、TNF法は小規模な系では競争力を持ち、構造的問題に対しては有望な結果を示した。
• SVDを用いた直交TNFは、根の分布をより適切に捉えており、解の位置付近で大きさが顕著に高くなる基底関数の挙動が、等高線図で可視化された。