[論文レビュー] A universal matrix-free split preconditioner for the fixed-point iterative solution of non-symmetric linear systems

Tom Vettenburg, Ivo M. Vellekoop|arXiv (Cornell University)|Jul 28, 2022

Matrix Theory and Algorithms被引用数 4

ひとこと要約

本稿では、任意の分解 $A = L + V$ に対して $\|V\| < 1$ を満たす場合に、非対称で発散的（accretive）な線形系 $A\bm{x} = \bm{y}$ に対して固定点反復の単調収束を保証する、普遍的な行列非依存の分割プリコンディショナーを提案する。特筆すべきは、唯一のプリコンディショナー形式 $(L + I)(I - V)^{-1}$ が、内部反復を伴わず、行列形式で評価可能であり、前方問題を排除することで反復コストを半減させつつ、多様な偏微分方程式（PDE）および遅延方程式において収束を保証できることである。

ABSTRACT

We present an efficient preconditioner for linear problems $A x=y$. It guarantees monotonic convergence of the memory-efficient fixed-point iteration for all accretive systems of the form $A = L + V$, where $L$ is an approximation of $A$, and the system is scaled so that the discrepancy is bounded with $\lVert V Vert&lt;1$. In contrast to common splitting preconditioners, our approach is not restricted to any particular splitting. Therefore, the approximate problem can be chosen so that an analytic solution is available to efficiently evaluate the preconditioner. We prove that the only preconditioner with this property has the form $(L+I)(I - V)^{-1}$. This unique form moreover permits the elimination of the forward problem from the preconditioned system, often halving the time required per iteration. We demonstrate and evaluate our approach for wave problems, diffusion problems, and pantograph delay differential equations. With the latter we show how the method extends to general, not necessarily accretive, linear systems.

研究の動機と目的

• 非対称線形系における固定点反復の不安定性と高メモリコストを解消し、単調収束を保証すること。
• 大規模な作用素の反復的逆行列計算を必要とする伝統的プリコンディショナーの限界を克服すること。
• 任意の分解 $A = L + V$ に対して $\|V\| < 1$ を満たす、行列非依存のプリコンディショナーを構築し、閉形式での評価が可能となるようにすること。
• プリコンディショニング系における前方問題の計算コストを排除し、反復あたりの時間コストを著しく削減すること。
• 非発散的（non-accretive）な一般系に対しても、同等の発散的形に変換することでこの手法を拡張すること。

提案手法

• 任意の $A = L + V$ および $\|V\| < 1$ を満たす前提で、$\Gamma = (L + I)(I - V)^{-1}$ の形を有する普遍的プリコンディショナーを提案し、$\|I - \Gamma^{-1}A\| < 1$ を満たすことで固定点反復の単調収束を保証する。
• 逆行列が閉形式で得られる $L$ を選べることを活用し、反復ソルバーを用いずにプリコンディショナーを評価可能とする。
• プリコンディショニング系から前方問題を排除する変換を導出し、反復コストを最大50％まで削減する。
• 提案されたプリコンディショナーが、発散的系の任意の分解に対して単調収束を保証する唯一の形式であることを証明する。
• 波動方程式、拡散方程式、遅延微分方程式にこの手法を適用し、多様なPDEにわたり高い頑健性と効率性を示す。
• 作用素ノルムの上限とスペクトル解析を用いて収束速度の推定値を導出し、一般系では $\|M\| \leq 1 - \frac{1}{S}$、ヘルミート系ではよりタイトな境界を導出する。

実験結果

リサーチクエスチョン

• RQ1 任意の分解 $A = L + V$ に対して $\|V\| < 1$ を満たす行列非依存のプリコンディショナーを構築でき、固定点反復の単調収束を保証できるか？
• RQ2 任意の分解選択に関わらず、すべての発散的系に対して収束を保証する唯一のプリコンディショナー形式は存在するか？
• RQ3 精度を損なわず、プリコンディショニング反復から前方問題の計算コストを排除できるか？
• RQ4 普遍的分割プリコンディショナーの収束速度は、シフト分割法や他の標準的プリコンディショナーと比較して、さまざまなPDEにおいてどのように異なるか？
• RQ5 非発散的系に対しても、同等の発散的形に変換することでこの手法を拡張可能か？

主な発見

• 固定点反復の単調収束をすべての発散的系で保証する唯一の形式は、$(L + I)(I - V)^{-1}$ である。
• 本手法により、内部反復を伴わず、行列非依存での評価が可能となり、高い効率性と低メモリ使用量が実現される。
• 実際の応用では、プリコンディショニング系から前方問題を排除でき、反復あたりのコストを約50％まで削減できる。
• 収束速度の境界が導出された：一般系では $\|M\| \leq 1 - \frac{1}{S}$、ヘルミート系では $\|M\| \leq 1 - \frac{1}{S}$、ここで $S = \|A^{-1}\|\|V\|$ である。
• 数値実験により、GMRES、BiCGSTAB、固定点反復の各手法が、ヘルムホルツ方程式、拡散方程式、パンタグラフ遅延方程式において有効に収束することが確認された。
• 非発散的系に対しても、同等の発散的形に変換することで本手法を拡張可能であり、適用範囲が広がった。