線形連立方程式を解くためのクリロフ部分空間法の収束の背後にある原理は何ですか?


24

私が理解しているように、線形連立方程式を解くための反復法には2つの主要なカテゴリがあります。

  1. 定常法(Jacobi、Gauss-Seidel、SOR、Multigrid)
  2. クリロフ部分空間法(共役勾配法、GMRESなど)

ほとんどの定常手法は、エラーのフーリエモードを繰り返し緩和(平滑化)することで機能することを理解しています。私が理解しているように、共役勾配法(クリロフ部分空間法)は、番目の残差に適用される行列のべきから最適な探索方向のセットを「ステップ実行」することによって機能します。この原理は、すべてのクリロフ部分空間法に共通ですか?そうでない場合、クリロフ部分空間法の収束の背後にある原則を一般的にどのように特徴付けますか?n


2
定常法の分析は、単純なモデルの問題によって偏っています。これらはフーリエモードの観点から分析できるためです。また、交互方向の暗黙的(ADI)および他の多くの方法を無視します。ほとんどの「定常法」のポイントは、多くの単純な「近似部分」ソルバーを1つの 反復ソルバーに結合することです。Krylovメソッドのポイントは、特定の定常線形反復の収束を加速(または強制)することです。
トーマスクリンペル

4
あなたの質問に答えるために書かれたと思う論文は、イプセンとマイヤー、クリロフ法の背後にあるアイデア、アメールです。数学。月刊105(1998)pp。889-899。それは、ここで入手できる、素晴らしくよく書かれた明確な論文です
アンドリューT.バーカー

@ AndrewT.Barker:素晴らしい!アンドリューありがとう!:)
ポール

回答:


21

一般に、すべてのクリロフ法は、行列のスペクトルで評価したときに小さい多項式を本質的に探します。特に、クリロフ法の番目の残差(初期推定がゼロ)は、次の形式で記述できます。n

rn=Pn(A)b

ここで、は次数モニック多項式です。 nPnn

場合有する、対角化可能である、我々はA = V Λ V - 1AA=VΛV1

rnVPn(Λ)V1b=κ(V)Pn(Λ)b.

た場合に、(例えば、対称またはユニタリ)正常で我々は知っている CGが異なる内積を使用して多項式を構築しながら(参照、アーノルディ反復介してGMRES構築ような多項式この回答を詳細)。同様に、BiCGは非対称Lanczosプロセスを介して多項式を構築しますが、チェビシェフ反復はスペクトルに関する事前情報(通常、対称定行列の最大および最小の固有値の推定)を使用します。κ V = 1。Aκ(V)=1.

クールな例(Trefethen + Bauが動機)として、スペクトルが次のようなマトリックスを考えます。

行列のスペクトル

MATLABでは、次を使用してこれを構築しました。

A = rand(200,200);
[Q R] = qr(A);
A = (1/2)*Q + eye(200,200);

次数すべてのモニック多項式にわたって実際残差を最小化する多項式を構築するGMRESを考慮すると、候補多項式を調べることで簡単に残差履歴を予測できますn

Pn(z)=(1z)n

私たちの場合、それは与える

|Pn(z)|=12n

以下のためののスペクトルにおける。AzA

ランダムなRHSでGMRESを実行し、残差履歴をこの多項式と比較すると、それらは非常によく似ているはずです(ため、多項式の候補値はGMRES残差よりも小さくなります)。b2>1

残存履歴


「マトリックスのスペクトル上で小さい」とはどういう意味かを明確にしていただけますか?
ポール

2
複素多項式としてと、多項式はのスペクトルを含む複素平面の領域に小さな係数を持ちます。固有値の散布図に重ねられた等高線図を想像してください。どれくらい小さいですか?問題、が正常かどうか、および右側に依存しただし基本的な考え方は、多項式のシーケンスがスペクトル上で次第に小さくなり、答えの残差推定値がなる傾向があることです。 A A b P n0PnAAb.(Pn)0
-Reid.Atcheson

@ Reid.Atcheson:非常にうまくいっています。書くことをお勧めします として、それが正規行列用であることを言及していますか?κ V VV1κ(V)
ジャックポールソン

最適なSORで前処理されたラプラシアンは、この行列例と非常に似たスペクトルを持っています。詳細はこちら:scicomp.stackexchange.com/a/852/119
ジェドブラウン

厳密に言えば、CGNEは特異値のみに依存するため、スペクトルから独立しています。
ジェドブラウン

17

規範について

Reid.Atchesonの回答の補遺として、規範に関するいくつかの問題を明確にしたいと思います。で反復、GMRES多項式見つけ最小にする残差のノルムを P n 2nthPn2

rn=Axnb=(Pn(A)1)bb=Pn(A)b.

仮定ので、SPDでありノルムを誘導し、そう。それからAAA1

rnA1=rnTA1rn=(Aen)TA1Aen=enTAen=enA

エラーを使用した場所

en=xnx=xnA1b=A1rn

こうして誤差のノルムは、と等価である残差のノルム。共役勾配は、誤差のノルムを最小化するため、低エネルギーモードの解決において比較的正確になります。GMRESが最小化する残差のノルムは、誤差のノルムに似ており、したがって、低エネルギーモードがあまりよく解決されないという意味で弱くなります。低エネルギーモードではさらに弱いため、残差のノルムは本質的に価値がないことに注意してください。A 1 A 2 A T A AAA1A2ATAA

収束境界の鮮明さ

最後に、さまざまなクリロフ法とGMRES収束の微妙さ、特に非正規演算子に関する興味深い文献があります。


あなたはOlavi Nevanlinnaことにより、優れた本をオフに左:books.google.com/...
マットKnepley

11

一言で言えば反復法:

  1. 静止法は本質的に固定点反復です:を解くには、可逆行列を選択し、固定点を見つけます。 これは、Banachの不動点定理によって収束します。さまざまなメソッドは、特定の選択に対応します(たとえば、Jacobi反復の場合、、ここでは対角要素を含む対角行列です)。、C 、X = X + C B - C A X I - C A < 1 C CAx=bC

    x=x+CbCAx
    ICA<1C D AC=D1DA
  2. クリロフ法の部分空間法は本質的に投影法にあります:部分空間を選択し、残差が直交するようにを探します。クリロフ法の場合、もちろんは、初期残差に適用される累乗にまたがる空間です。さまざまな方法は、特定の選択に対応します(たとえば、CGの場合は、GMRESの場合は)。XU B - A X V U A V V = U V = A UU,VCnx~UbAx~VUAVV=UV=AU

    Vx~UU

    UAx~

    これは、反復法に関する Youcef Saadの本でうまく説明されています

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.