[論文レビュー] New Polynomial Preconditioned GMRES.

Jennifer Loe, Ronald B. Morgan|arXiv (Cornell University)|Nov 16, 2019

Matrix Theory and Algorithms参考文献 23被引用数 4

ひとこと要約

本稿では、最小残差多項式から導かれる新しい多項式予約子を、GMRESに導入し、大規模な線形系の効率的かつ安定的な解法を可能にする。収束を著しく加速し、ドット積の回数を削減する。BiCGStab や FGMRES より優れた性能を示し、特に困難な問題において顕著な効果を発揮する。根の追加と専用の安定性チェックにより、安定性が向上する。

ABSTRACT

A new polynomial preconditioner is given for solving large systems of linear equations. The polynomial is derived from the minimum residual polynomial and is straightforward to compute and implement. It this paper, we study the polynomial preconditioner applied to GMRES, however it could be used with any Krylov solver. This preconditioning algorithm can dramatically improve convergence for difficult problems and can reduce dot products by an even greater margin. Stability control using added roots allows for high degree polynomials, and we give an additional check for stability. This approach is compared to existing methods such as BiCGStab and FGMRES.

研究の動機と目的

• 大規模で困難な線形系をGMRESで解く際に、収束を改善する新しい多項式予約子の開発。
• 計算と実装が簡単でありながら、数値的安定性を維持する予約子の設計。
• Krylov解法におけるドット積の回数を削減し、計算効率を向上。
• 追加された根を用いた安定性制御により、高次多項式の使用を可能に。
• 実用的で信頼性の高い性能を確保するための堅牢な安定性チェックの提供。

提案手法

• 予約子は最小残差多項式から導かれるため、各ステップで残差を最適に減少させる。
• この方法は、Krylov部分空間構造を変更せずに、GMRESフレームワーク内に多項式予約子を適用する。
• 多項式に追加の根を追加することで安定性が向上し、数値的崩壊を回避しながら高次多項式の使用が可能になる。
• 計算中に潜在的な不安定性を監視・制御するため、専用の安定性チェックを実装。
• この手法は一般化可能であり、GMRESに限らず、任意のKrylov部分空間法に適用可能。

実験結果

リサーチクエスチョン

• RQ1 困難な線形系において、GMRESの収束を著しく改善する多項式予約子はどのように構築できるか？
• RQ2 この新しい予約戦略により、GMRESにおけるドット積の回数はどの程度削減できるか？
• RQ3 安定性を損なわずに、高次多項式を効果的に予約に使用できるか？
• RQ4 提案された予約子は、BiCGStab や FGMRES といった既存手法と比較して、性能面でどのように差をつけるか？
• RQ5 追加された根と安定性チェックは、堅牢で高次多項式の予約を可能にする上で、果たす役割は何か？

主な発見

• 新しい多項式予約子は、標準的なGMRESと比較して、困難な線形系の収束を著しく加速する。
• 残差の減少割合よりも、ドット積の回数の削減が顕著に顕著であり、計算効率が向上する。
• 追加された根による安定性制御により、数値的問題を伴わずに高次多項式の使用が可能になる。
• 安定性チェックにより、多様な問題タイプにわたり信頼性の高い性能が確保され、感受性の高いケースでの崩壊を防止する。
• 困難な問題において、収束速度と堅牢性の両面で、BiCGStab や FGMRES を上回る性能を発揮する。