低ランクの変更はKrylovメソッドの収束にどのように影響しますか?


14

線形システムがあり、すべてのbに対して適切なクリロフ法(CGやGMRESなど)を使用して迅速に収束するとします。場合Bが低いランクを有する行列であり、システム上の同じクリロフ法であろうまた、(理想的にはおおよそにのみ依存する反復の余分な数で迅速に収束)?Ax=bbBA + B x = b rr(A+B)x=br

そのようなシステムの例は、十分に前処理された膜の弾性と曲げに加えて、高密度の外部製品構造を備えた前処理されていない気圧条件です。

あるため、質問は前提条件の有無にかかわらず同じであることに注意してください。P(A+B)Q=PAQ+PBQP A Qのランク修正で。rPAQ

回答:


7

Krylov部分空間がべき乗に基づいている場合、収束は多くても補正のランクの反復回数だけ遅延します。A T Aの累乗に基づいている場合、この数値の最大2倍です。AATA

私は文学でそのような声明を見たことがありません。しかし、最初の場合に妥当性を確認するために、あることを示すのに十分である行列のクリロフ空間番目A + U S V T U Vを有するr個の列が低ランクの修正なしではなく、Aに対応する空間内に収容されていますそれに対応して、より高いインデックスk + r。これは簡単に確認できます。kA+USVTU,Vrk+r


累乗に基づいて」とはどういう意味ですか?クリロフソルバは情報を与えられているA + Bのみではなく、A直接的。AA+BA
ジェフリーアーヴィング

気にしないでください:おそらくあなたは問題の行列のべき乗を意味するので、この場合はです。A+B
ジェフリーアーヴィング

はい。メソッドにはパラメーターとして行列があり、この行列は通常示されます。A
アーノルドノイマイアー

たぶん、さらに関心をあなたが上で、いくつかの要件にあなたの方程式(または溶液)を書き換える可能性が、X = E + Σ のk = 1A - 1つの B のk A - 1件の Bたのかもしれない助けている場合Bは冪零ですかA - 1 Bの小さなノルム。また、問題の解決策への依存を認識しています。Bバツ=E+k=1A1BkA1bBA1Bundisturbed
バスティアンエベリング
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.