[論文レビュー] Multiprecision path tracking

Daniel J. Bates, Andrew J. Sommese+1 more|ArXiv.org|May 3, 2006

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

ひとこと要約

本稿では、数値代数幾何におけるホモトピー続行法のための適応的精度パス追跡アルゴリズムを提案する。この手法は、多項式系の条件数に基づき、計算精度を動的に調整することで、特に悪条件または高次系において高い信頼性と効率性を達成する。チェビシェフ多項式（次数300まで）を用いた実験により、完全な精度が達成された。

ABSTRACT

A path tracking algorithm that adaptively adjusts precision is presented. By adjusting the level of precision in accordance with the numerical conditioning of the path, the algorithm achieves high reliability with less computational cost than would be incurred by raising precision across the board. We develop simple rules for adjusting precision and show how to integrate these into an algorithm that also adaptively adjusts the step size. The behavior of the method is illustrated on several examples arising as homotopies for solving systems of polynomial equations.

研究の動機と目的

• 固定精度パス追跡の非効率性と信頼性の欠如、特に悪条件系や高次系を含む問題に対して、これを是正すること。
• すべてのパスに均一に高精度算術を適用するのを避けることで、計算コストを削減し、数値的条件数に基づき各パスごとに精度を適応的に調整すること。
• 特異点付近や高全次数系において、標準的な倍精度では失敗するが、信頼性のある解パスの追跡を可能にすること。
• 精度適応と適応的ステップサイズ制御を統合し、強固で効率的なパス追跡を実現すること。
• 高次チェビシェフ多項式のような困難な問題、特に厳密な精度要件を満たす問題に対して、本手法の有効性を示すこと。

提案手法

• アルゴリズムは、パスに沿ったホモトピー関数の条件数に基づき、浮動小数点算術で使用する桁数を動的に調整する。
• 精度適応と適応的ステップサイズ制御を統合し、修正子の収束行動に基づきステップサイズを調整する。
• 予測子・修正子フレームワークを用い、Euler法または高次の予測子とニュートン法ベースの修正子を採用し、正しい解パスを分離する制約を課す。
• ヤコビアンと残差ノルムの推定値を用いて、条件数を評価し、精度の増加を誘導する。
• 失敗を回避するために、余剰な桁数を追加する安全マージンを適用し、パラメータは経験的テストにより調整する。
• 精度が不十分であると検出された場合、収束が達成されるまで段階的に精度を向上させる。

実験結果

リサーチクエスチョン

• RQ1 パス追跡中に精度をどのように適応的に調整すれば、均一な高精度算術のコストを負わずに信頼性を向上させられるか？
• RQ2 標準的な倍精度がパス追跡で失敗する条件は何か？ そして、その場合に効率的に精度を向上させる方法は？
• RQ3 適応的精度制御を適応的ステップサイズ制御と効果的に組み合わせることで、信頼性と効率性の両方を向上させられるか？
• RQ4 本手法は、根が密集しているか悪条件の高次多項式系に対して、どのように性能を発揮するか？
• RQ5 適応的精度は、特異な終点に対してエンドゲーム手法を置き換えたり補完したりできるか、その範囲はどの程度か？

主な発見

• 適応的精度アルゴリズムは、次数とパス依存の精度レベルを用いて、次数300のチェビシェフ多項式をすべての根が少なくとも10桁の精度で正しく解いた。
• 次数300の場合、x=1.0付近で最も悪条件の強いパスに対して、最大で約3.4×10^24ビットの精度を要した。これは、根の密集に極めて敏感であることを示している。
• 最も高い精度を要したパスは、常にx=1.0付近に終点をもつパスであり、その領域では根の間隔が最小であったためである。
• 修正子ループの失敗を避けるために、初期ステップサイズを小さくする必要があった。ステップサイズを大きくすると、精度要請が急増し、収束しなくなった。
• 全テスト次数（10から300）において、すべてのパスで完全な精度が達成され、同じ次数の異なるパス間で精度レベルに顕著な差が見られた。
• 良好に条件付けられたセグメントでは不要な高精度計算を回避し、悪条件のセグメントでは信頼性を保証するため、固定精度アプローチを上回る性能を発揮した。