[論文レビュー] Certifying solutions to overdetermined and singular polynomial systems over Q

Tulay Ayyildiz Akoglu, Jonathan D. Hauenstein+1 more|arXiv (Cornell University)|Aug 12, 2014

Polynomial and algebraic computation参考文献 41被引用数 6

ひとこと要約

この論文は、ℚ 上の過決定および特異な多項式系の近似解を証明するためのハイブリッド記号的・数値的手法を提示する。数値的根から正確な有理数一変数表現（RUR）を計算することで、ニュートン反復と正確算術を用いて、悪条件な系を適切に制約されたRURに変換し、α理論による証明を可能にする。多項式時間の計算量と、RURが得られると早期終了する。

ABSTRACT

This paper is concerned with certifying that a given point is near an exact root of an overdetermined or singular polynomial system with rational coefficients. The difficulty lies in the fact that consistency of overdetermined systems is not a continuous property. Our certification is based on hybrid symbolic-numeric methods to compute the exact "rational univariate representation" (RUR) of a component of the input system from approximate roots. For overdetermined polynomial systems with simple roots, we compute an initial RUR from approximate roots. The accuracy of the RUR is increased via Newton iterations until the exact RUR is found, which we certify using exact arithmetic. Since the RUR is well-constrained, we can use it to certify the given approximate roots using alpha-theory. To certify isolated singular roots, we use a determinantal form of the "isosingular deflation", which adds new polynomials to the original system without introducing new variables. The resulting polynomial system is overdetermined, but the roots are now simple, thereby reducing the problem to the overdetermined case. We prove that our algorithms have complexity that are polynomial in the input plus the output size upon successful convergence, and we use worst case upper bounds for termination when our iteration does not converge to an exact RUR. Examples are included to demonstrate the approach.

研究の動機と目的

• 一貫性が連続的でない、普遍的な下界があまりに楽観的すぎるという問題に対処すること。
• 過剰に保守的な境界を避けるために、近似根から直接正確な有理数一変数表現（RUR）を計算する実用的な証明手法を開発すること。
• 定性的な系の変形と反復的精錬を用いて、特異な根を単純な根に変換することで、α理論の証明を過決定および特異な系に拡張すること。
• 数値的根の探索と正確算術を組み合わせることで、収束時に入力および出力サイズの多項式時間の計算量を達成し、正しさと効率性を保証すること。

提案手法

• 過決定または特異な系に対して、ホモトピー続行法を用いてランダムに選ばれた適切に制約された部分系の近似根を計算する。
• ヒューリスティックなフィルタリングを用いて、孤立したまたは有理数的成分に対応する候補の近似根を選択する。
• 選択された近似根からラグランジュ補間と有理数再構成を用いて初期のRURを構築する。
• 局所的またはグローバルなニュートン反復を用いてRURを精錬し、正確算術により正確なRURに到達するまで反復する。
• 行列式形式による同一型の退化を用いて、変数を追加せずに多項式を追加し、特異な根を単純な根に変換する。
• 最終的なRURをα理論またはヘルミート行列を用いて証明し、元の近似根が正確な解の近傍にあることを確認する。

実験結果

リサーチクエスチョン

• RQ1 ℚ 上の過決定多項式系の近似根を、楽観的な普遍的下界ではなく、正確なRURを用いて高い信頼性で証明することは可能か？
• RQ2 多項式系における特異な根を、標準的な証明手法を可能にするために、どのように単純な根に変換できるか？
• RQ3 近似根から正確なRURを計算する計算量は何か？ また、実世界の系において実用的に行えるか？
• RQ4 ハイブリッド記号的・数値的手法は、適切に制約されたRURにおいて、入力および出力サイズの多項式時間収束を達成できるか？

主な発見

• 提案手法は、正しい根候補が選択された場合に、ニュートン反復を用いて正確なRURを正確に計算でき、収束が保証される。
• キャプラス系では、各次数2の有理数的成分が4つ特定され、分解体 ℚ[√3] が確認され、8つの重複根の構造が裏付けられた。
• サイクリック-4系では、有理数係数の12次RURが計算され、8つの埋め込まれた同一型の特異点と2つの非可約曲線が正しく特定された。
• 162個の解を有するサイクリック-9系では、54個の正則点および54個の二重点に対して、高次のRUR（最大12次）が計算され、イデアルの構造が確認された。
• ワーニングサイズが小さい場合に早期終了が達成され、普遍的下界よりも実用的であることが示された。
• 収束に成功した場合、計算量は入力および出力サイズの多項式時間であり、非収束ケースに対しても最悪ケースの境界が提供された。