クリロフ部分空間

提供: testwiki
2025年2月17日 (月) 19:24時点における133.86.227.82 (トーク)による版 (参考文献)
(差分) ← 古い版 | 最新版 (差分) | 新しい版 → (差分)
ナビゲーションに移動 検索に移動

クリロフ部分空間(クリロフぶぶんくうかん、テンプレート:Lang-en線型代数において、テンプレート:Mvar正方行列テンプレート:Mvarテンプレート:Mvarベクトルテンプレート:Mvarによって生成されるテンプレート:Mvar次クリロフ部分空間は、テンプレート:Mvarテンプレート:Mvarのべき乗の像が張る線型部分空間である。

𝒦r(𝑨,𝒃)=span{𝒃,𝑨𝒃,𝑨2𝒃,,𝑨r1𝒃}.

ロシアの応用数学者で海軍技術者であったテンプレート:Illにちなんで名づけられた。

大規模疎行列の1個または少数の固有値の計算や、大規模な連立一次方程式の求解に用いられる現代的な反復法では、行列を消去法などで順次変型すると疎行列の構造が崩れてしまい次第に密化するので演算量や記憶を保持する量が共に増大してしまい,ついには扱いきれなくなりがちである。そこでクリロフ系の解法では,元の疎行列を変型せずに,ベクトルに対する線型の作用素としてだけ利用する。つまり与えられたベクトルに対して行列を乗じるという計算を,行列の疎性を活かして(行列が対称であれば対称性も)行うのである。 𝒃を初期ベクトルとすると、𝑨を順に掛けて𝑨𝒃𝑨2𝒃を得るといった方法を取る。このようなアルゴリズムを総称して、クリロフ部分空間法と呼ぶ。これは数値線形代数において最も成功した手法の一つである。

ベクトル列は急速に線型従属に近づくため、クリロフ部分空間を用いる方法には、エルミート行列に対してはランチョス法、非エルミート行列に対してはアーノルディ法などの直交化手法が含まれることが多い。

主なクリロフ部分空間法として、アーノルディ法ランチョス法GMRES法(generalized minimum residual)、 BiCGSTAB法 (stabilized biconjugate gradient, 共役勾配法の一つ)、QMR法 (quasi minimal residual)[1][2]、 TFQMR法 (transpose-free QMR)[3]、MINRES法 (minimal residual)[4] などが知られている。

出典

テンプレート:Reflist

参考文献

  • テンプレート:Cite book
  • Charles George Broyden and Maria Teresa Vespucci(2004): Krylov Solvers for Linear Algebraic Systems, Elsevier(Studies in Computational Mathematics 11), ISBN 0-444-51474-0.
  • David S. Watkins (2008). The Matrix Eigenvalue Problem: GR and Krylov Subspace Methods, SIAM.
  • Liesen, J. and Strakos, Z. (2012). Krylov subspace methods: principles and analysis. OUP Oxford.
  • Gerard Meurant and Jurjen Duintjer Tebbens (2020). Krylov methods for nonsymmetric linear systems - From theory to computations, Springer Series in Computational Mathematics, vol.57. テンプレート:ISBN2, テンプレート:DOI.
  • Iman Farahbakhsh: Krylov Subspace Methods with Application in Incompressible Fluid Flow Solvers, Wiley, ISBN 978-1-11961868-3 (2020).
  • 藤野清次、阿部邦美、杉原正顯、中嶋徳正:「線形方程式の反復解法」、丸善、 ISBN 978-4-62108741-1(2013年9月27日)。

外部リンク

テンプレート:Linear algebra テンプレート:Math-stub

  1. Freund, R. and Nachtigal, N. "QMR: A Quasi-Minimal Residual Method for Non-Hermitian Linear Systems." Numer. Math. 60, 315-339, 1991.
  2. Freund, R. and Nachtigal, N. "An Implementation of the QMR Method Based on Coupled Two-Term Recurrences." SIAM J. Sci. Statist. Comput. 15, 313-337, 1994.
  3. Roland W. Freund, A transpose-free quasi-minimal residual algorithm for non-Hermitian linear systems, en:SIAM Journal on Scientific Computing 1993; 14(2):470–482.
  4. Christopher C. Paige and Michael A. Saunders, Solution of sparse indefinite systems of linear equations, en:SIAM Journal on Numerical Analysis 1975; 12(4):617–629.