[論文レビュー] PHCpack in Macaulay2

Elizabeth Gross, Sonja Petrović+1 more|May 24, 2011

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

ひとこと要約

本稿では、数値ホモトピー続行法を用いて多項式系ソルバであるPHCpackを統合したMacaulay2パッケージ、PHCpackを提示する。このパッケージにより、大規模な多項式系に対する効率的な数値解法が可能になる。ブラックボックスソルブ、ミックスドボリューム計算、解のフィルタリング、パス追跡、数値的既約分解を提供し、21変数・67解までの系において性能が実証されている。

ABSTRACT

The Macaulay2 package PHCpack.m2 provides an interface to PHCpack, a general-purpose polynomial system solver that uses homotopy continuation. The main method is a numerical blackbox solver which is implemented for all Laurent systems. The package also provides a fast mixed volume computation, the ability to filter solutions, homotopy path tracking, and a numerical irreducible decomposition method. As the size of many problems in applied algebraic geometry often surpasses the capabilities of symbolic software, this package will be of interest to those working on problems involving large polynomial systems.

研究の動機と目的

• 大規模な多項式系の数値的解法を目的とした、Macaulay2とPHCpackの間の高性能インターフェースを提供すること。
• 応用代数幾何学における大規模問題を扱う際の記号的ソフトウェアの限界を是正すること。
• 解のフィルタリング、パス追跡、既約分解などの高度な数値的代数幾何学的タスクを支援すること。
• トーラス内における孤立解の数に対する上界としてのミックスドボリュームの効率的計算を可能にすること。
• 完全な診断フィードバックを備えた高レベルのMacaulay2インターフェースを通じて、低レベルのPHCpack機能にアクセス可能にすること。

提案手法

• 既知の解を持つスタート系の解から、ターゲット系の解に至るパスを数値的ホモトピー続行法で追跡する。
• ユーザーが複雑さを抽象化できるブラックボックスソルバを採用し、自動的にスタート系とホモトピーを構築する。
• ポリトープ的メソッドを用いた高速なミックスドボリューム計算を実装し、標準的および安定的ミックスドボリュームの両方をサポートする。
• ゼロフィルターやノンゼロフィルタ一を用いて、座標の大きさに基づき解を分離する解のフィルタリングを提供する。
• 係数を高精度で再計算することで、解の高精度化を可能にする refineSolutions をサポートする。
• モノドロミーに基づく分解を用いて数値的既約分解を実現し、各既約成分に対してWitnessSetsを返す。

実験結果

リサーチクエスチョン

• RQ1 数値的ホモトピー続行法を、大規模な多項式系を解くためのMacaulay2環境に効果的に統合する方法は何か？
• RQ2 ミックスドボリューム計算は、トーラス内における孤立解の数に対する高速かつ信頼性の高い上界として、どの程度有効に機能するか？
• RQ3 PHCpackによる数値的既約分解は、多項式系によって定義される代数的多様体の次元と次数を正確に回復できるか？
• RQ4 有理数係数をもつ系および高精度の解の精錬処理において、このインターフェースはどの程度効率的に動作するか？
• RQ5 大規模なゼロ次元系において、数値的手法が記号的計算に比べてどの程度の性能向上をもたらすか？

主な発見

• solveSystem関数は21変数の系について0.056秒で67個の孤立解を正しく計算し、系の次数が67であることを確認した。
• 同じ系におけるミックスドボリューム計算は0.011秒で完了し、上界として75を算出し、実際の孤立解数を上回った。
• 同じ系についての記号的次数計算には767.432秒を要し、数値的手法の性能優位性が顕著に示された。
• refineSolutionsを用いた解の精錬により、解の12番目の座標が2.23e-20から1.43e-67にまで低下し、数値的安定性と精度が確認された。
• 数値的既約分解により、次元5で次数が4, 2, 2の3つの既約成分が正しく同定され、記号的プライマリ分解と一致した。
• インターフェースは一時ファイルを通じてPHCpackの全機能にアクセス可能であり、ユーザーが生の計算ログや診断情報を確認できる。