[論文レビュー] Error Bounds for the Krylov Subspace Methods for Computations of Matrix Exponentials

Hao Wang, Qiang Ye|arXiv (Cornell University)|Mar 23, 2016

Matrix Theory and Algorithms参考文献 22被引用数 3

ひとこと要約

本稿では、行列指数関数 $e^{-\tau A}v$ の計算におけるKrylov部分空間法について、$A$ のエルミート部および反エルミート部の固有値を用いて長方形領域による場の値の包囲を定義することで、新たな事前および事後誤差境界を提示する。これらの境界は、超線形収束および初期停滞を説明し、Saad よりも、Hochbruch-Lubich よりも顕著に改善された古典的境界を凌駕する。

ABSTRACT

In this paper, we present new a posteriori and a priori error bounds for the Krylov subspace methods for computing $e^{-τA}v$ for a given $τ>0$ and $v \in C^n$, where $A$ is a large sparse non-Hermitian matrix. The {\em a priori} error bounds relate the convergence to $λ_{\min}\left(\frac{A+A^*}{2} ight)$, $λ_{\max}\left(\frac{A+A^*}{2} ight)$ (the smallest and the largest eigenvalue of the Hermitian part of $A$) and $|λ_{\max}\left(\frac{A-A^*}{2} ight)|$ (the largest eigenvalue in absolute value of the skew-Hermitian part of $A$), which define a rectangular region enclosing the field of values of $A$. In particular, our bounds explain an observed superlinear convergence behavior where the error may first stagnate for certain iterations before it starts to converge. The special case that $A$ is skew-Hermitian is also considered. Numerical examples are given to demonstrate the theoretical bounds.

研究の動機と目的

• 大規模なスパースな非エルミート行列に対して $e^{-\tau A}v$ を計算するKrylov部分空間法の、より鋭い事前および事後誤差境界の開発。
• 行列のエルミート部および反エルミート部のスペクトル特性を用いて、観察された超線形収束挙動（初期停滞を含む）を説明すること。
• 円形またはノルムに基づく推定ではなく、場の値の長方形包囲を用いることで、既存の境界（特にSaad よりも Hochbruch-Lubich よりも）を改善すること。
• バンド行列関数の減衰特性に基づく、計算可能で鋭い事後誤差推定を提供すること。

提案手法

• バンド行列関数の減衰特性を用いて、事後誤差境界を導出することで、鋭く計算可能な誤差推定を可能にする。
• 行列 $A$ の場の値を包囲する長方形領域上で指数関数のファーベル多項式近似を構築する。長方形は $a = \lambda_{\min}\left(\frac{A+A^*}{2}\right)$, $b = \lambda_{\max}\left(\frac{A+A^*}{2}\right)$, および $c = \left|\lambda_{\max}\left(\frac{A-A^*}{2}\right)\right|$ で定義される。
• 長方形領域へのコンformal写像をJacobiの楕円関数を用いて構築し、これを簡略化することで取り扱いやすい誤差境界を導出する。
• 非エルミート正定値および反エルミート行列に対して、簡略化された事前境界を確立し、収束速度が包囲長方形のサイズおよび形状にどのように関連するかを明らかにする。
• さまざまな $\tau$、行列タイプ（非エルミート正定値、反エルミート）およびKrylov反復回数を用いた数値的検証により、境界の有効性を検証する。
• 本稿の境界を古典的境界（Saad、Hochbruch-Lubich）と比較し、タイトネスおよび正確性の顕著な向上を示す。

実験結果

リサーチクエスチョン

• RQ1 場の値の長方形包囲を用いることで、Krylov部分空間法の行列指数関数計算における事前誤差境界をどのように改善できるか？
• RQ2 Krylov近似の誤差が、収束に至るまでにしばしば初期停滞を示すのはなぜか？これは行列 $A$ のスペクトル特性で説明可能か？
• RQ3 バンド行列の減衰特性に基づいて、計算可能で鋭い事後誤差境界を構築できるか？
• RQ4 本稿の境界は、Saad よりも Hochbruch-Lubich の境界と比較して、タイトネスおよび実用的有用性の観点でどのように差をつけるか？
• RQ5 反エルミート行列（$A = iH$）の場合は、解がノルムを保存するため、境界は収束挙動をどの程度説明できるか？

主な発見

• 事前誤差境界は、$A$ のエルミート部および反エルミート部の極端な固有値を用いて導出され、場の値長方形による幾何的解釈が可能となる。
• 境界は観察された超線形収束および初期停滞を説明する：収束は通常 $k \approx 2\tau\rho$ で開始され、ここで $\rho$ はエルミート部のスペクトルギャップである。
• ${\tau = 50}$ の場合、事前境界はSaadのノルムに基づく境界よりも著しくタイトであり、$\|\tau A\|$ が大きいために過剰に悲観的になる古典的境界を凌駆する。
• 事後誤差推定は、すべてのテストケースで鋭く、非エルミート正定値および反エルミート行列の両方の例で真の誤差を一貫して追跡する。
• HochbruchとLubichの円形包囲境界よりも、特に大きな $\tau$ に対して、本稿の境界は場の値のよりタイトな長方形包囲により顕著に改善される。
• 反エルミート行列の場合（$A = iH$）、境界は観察された収束曲線をよく再現し、初期停滞に続いて急激な減衰が生じるが、収束開始の予測時刻 $k \approx 2\tau\rho$ が数値結果と一致する。