[論文レビュー] alphaCertified: certifying solutions to polynomial systems

Jonathan D. Hauenstein, Frank Sottile|arXiv (Cornell University)|Nov 4, 2010

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

ひとこと要約

この論文は、スメールのα理論を用いて多項式系の数値解を厳密に認証する計算ツールである alphaCertified を紹介する。このツールは、正確な有理数または任意精度浮動小数点演算を用いて、ニュートン法の2次収束を検証し、異なる解を区別し、実解を特定する。正方系および過剰決定系への応用があり、特にモーメント曲線への接平面を含むシューベルト問題において126個の実解を検証した。

ABSTRACT

Smale's alpha-theory uses estimates related to the convergence of Newton's method to give criteria implying that Newton iterations will converge quadratically to solutions to a square polynomial system. The program alphaCertified implements algorithms based on alpha-theory to certify solutions to polynomial systems using both exact rational arithmetic and arbitrary precision floating point arithmetic. It also implements an algorithm to certify whether a given point corresponds to a real solution to a real polynomial system, as well as algorithms to heuristically validate solutions to overdetermined systems. Examples are presented to demonstrate the algorithms.

研究の動機と目的

• 多項式系の数値解をヒューリスティックに検証する数値ホモトピー・ソルバーにおいて、厳密な認証の欠如に対処すること。
• 正確な有理数または任意精度浮動小数点算術を用いて、数学的に厳密な数値解の認証を提供するα理論のアルゴリズムを実装すること。
• 近似解が実解に対応するかどうかを認証できるようにすること、特に実多項式系に対して特に。
• 固定点が実際に解でない可能性のある過剰決定系に対して、ランダムな正方部分系を用いたヒューリスティックな検証により、認証を拡張すること。
• 数値代数幾何学の研究者にとって、数学的確実性をもって解集合を検証できる実用的で使いやすいツールを提供すること。

提案手法

• alphaCertified は、スメールのα理論から得られる α(f,x) 定数を計算し、多項式系 f の点 x においてニュートン法が2次収束するかどうかを判定する。
• 正方多項式系の解に対して完全に厳密な認証を行うために、正確な有理数算術を用いる。
• 高精度の認証には、誤差の制御された境界を有する任意精度浮動小数点算術を用いて α(f,x) を評価する。
• 解が実数であるかどうかを、系の実構造と複素共役に関するニュートン反復の挙動を分析することで特定する。
• 過剰決定系（N > n）に対しては、複数のランダムな正方部分系を生成し、それぞれにα理論を適用して結果を比較し、近似解が与えられたδ内に近接していることをヒューリスティックに検証する。
• 複雑な代数幾何学的問題の入力ファイル作成を支援するために、Maple と統合されたソフトウェアを提供する。

実験結果

リサーチクエスチョン

• RQ1 α理論は、多項式系の数値解に対する厳密な認証を効果的に実装できるか？
• RQ2 alphaCertified における正確な有理数算術は、正方系の解認証において数学的正しさを保証できるか？
• RQ3 α理論は、多項式系において異なる解をどれほど正確に区別でき、実解を特定できるか？
• RQ4 固定点が実際に解でない可能性のある過剰決定系において、α理論をヒューリスティックに解の検証に適応できるか？
• RQ5 モーメント曲線への接平面を含むシューベルト問題における実解の数の下界は何か？

主な発見

• alphaCertified は、すべての接平面が実数である1000件のシャピロ予想の例において、正確に126個の実解を持つことを認証し、この設定における予想の妥当性を確認した。
• k=1 の場合、実数と複素共役の混合な接平面を含む42,000件の例では、最小100個の実解が検出され、非自明な下界の仮説を支持した。
• k=2 の場合、24,000件の例では最小100個の実解が検出され、1,000件の例では126個の実解が得られ、下界仮説の強い一貫性を示した。
• k=3 の23,500件の例では最小100個の実解が検出され、k=4 の22,500件の例でも最小100個の実解が検出され、100解の下界の存在を強化した。
• 256ビット精度の下で、シャピロ予想の1,000件のテストケースにおいて、すべての126個の解が正しく実数として特定された。これは高い信頼性を示している。
• k=0 から 4 までの100,000件の例における実解数の頻度表は、一貫した下界を示し、このようなシューベルト問題における最小実解数の予想を支持した。