[論文レビュー] Khovanskii-Rolle continuation for real solutions

D. Bates, Frank Sottile|ArXiv.org|Aug 31, 2009

Polynomial and algebraic computation被引用数 5

ひとこと要約

本稿では、多項式系の非退化な実解をすべて同定するための新しい数値的手法、Khovanskii-Rolle継続法を紹介する。この手法は、Gale双対表現における実曲線を追跡することで、複素解を完全に回避する。ホモトピー法とは異なり、この手法は少数項多項式（fewnomial）の境界に依存し、複素解の数ではなく実解の数に依存して、少数項多項式系において著しく効率的である。

ABSTRACT

We present a new continuation algorithm to find all nondegenerate real solutions to a system of polynomial equations. Unlike homotopy methods, it is not based on a deformation of the system; instead, it traces real curves connecting the solutions of one system of equations to those of another, eventually leading to the desired real solutions. It also differs from homotopy methods in that it follows only real paths and computes no complex solutions of the original equations. The number of curves traced is bounded by the fewnomial bound for real solutions, and the method takes advantage of any slack in that bound.

研究の動機と目的

• 多項式系の実解のみを求める数値的継続法を開発し、複素解の計算を回避すること。
• 少数項理論とKhovanskiiによるロルの定理の一般化を活用し、実解の数を上限づけ、経路追跡を誘導すること。
• 実解の数が複素解の数に比べてはるかに小さい少数項多項式系において、効率的な手法を設計すること。
• Gale双対系における実曲線を追跡するプロトタイプアルゴリズムを実装し、複素数演算を一切用いずにすべての実解を発見することを保証すること。
• 次元の呪いに見られる他の実解法とは異なり、少数項境界の余剰（slack）を自然に活用するフレームワークを提供すること。

提案手法

• 元の $ n $ 変数の $ n $ 個の多項式、$ l+n+1 $ 個の単項式を持つ系を、$ \mathbb{R}^l $ 上の $ l $ 個のマスターファンクション系、$ l+n+1 $ 個の超平面の補集合上で定義されたGale双対系に変換する。
• 元の系の実解とGale双対系の実解との間の全単射を用い、非退化性と正の性質を保存する。
• Gale双対空間における実経路に沿って曲線追跡を実行し、単体 $ \Delta $ の頂点から出発する。予測子・修正子ステップを用い、ステップ長を動的に制御し、ニュートン修正を適用する。
• 曲線が特異となる頂点では、標準的手法が失敗するため、$ y_2 = \alpha \cdot y_1^\beta $ の単項式近似を用いて、単項式曲線上の近傍点から曲線追跡を開始する。
• 各実解が、追跡された弧の両端点として2回発見されることを保証し、出力の基本的検証チェックを可能にする。
• 事前計算はBertiniを用いるが、核心的な経路追跡はMapleで実装され、すべて実数演算で処理される。

実験結果

リサーチクエスチョン

• RQ1 複素解を計算せずに、多項式系の実解のみを求める数値的継続法を設計できるか？
• RQ2 実解の少数項境界を活用することで、数値的継続における経路追跡数を削減できるか？
• RQ3 Gale双対性は、実曲線追跡に適した同値な系に多項式系を変換する役割を果たすか？
• RQ4 単体 $ \Delta $ の頂点のような特異点において、標準的手法が失敗する状況で、曲線追跡をどのように安定化できるか？
• RQ5 実数のみを扱う継続法は、少数項多項式系において、標準的なホモトピー継続法を上回る効率性を示せるか？

主な発見

• Khovanskii-Rolle継続法は、複素解を一切計算せず、実解のみを追跡するため、実解探索においてホモトピー法に比べて根本的に効率的である。
• 追跡する経路数は、複素解の数ではなく、実解の少数項境界によって上限づけられており、少数項多項式系では顕著な効率的向上が達成される。
• $ l=2 $ の場合、Mapleで実装されたプロトタイプ実装は、Bertini や PHCpack に比べ、テスト例において優れた性能を示した。
• 単項式近似を用いることで、単体の頂点における特異点を効果的に処理し、標準的な予測子・修正子追跡に滑らかに移行できる。
• 各解 $ S_i $ が弧の端点として2回発見されることを検証する組み込み検証チェックにより、曲線追跡の欠落を検出可能である。
• 理論的境界よりも実際の実解数が小さい場合、少数項境界の余剰を自然に活用でき、追跡する経路数が減少する。