Deflation techniques for an implicitly restarted Arnoldi iteration

被引:439
作者
Lehoucq, RB [1 ]
Sorensen, DC [1 ]
机构
[1] RICE UNIV, DEPT COMPUTAT & APPL MATH, HOUSTON, TX 77251 USA
关键词
Arnoldi method; Lanczos method; eigenvalues; deflation; implicit restarting;
D O I
10.1137/S0895479895281484
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
A deflation procedure is introduced that is designed to improve the convergence of an implicitly restarted Arnoldi iteration for computing a few eigenvalues of a large matrix. As the iteration progresses, the Ritz value approximations of the eigenvalues converge at different rates. A numerically stable scheme is introduced that implicitly deflates the converged approximations from the iteration. We present two forms of implicit deflation. The first, a locking operation, decouples converged Ritz values and associated vectors from the active part of the iteration. The second, a purging operation, removes unwanted but converged Ritz pairs. Convergence of the iteration is improved and a reduction in computational effort is also achieved. The deflation strategies make it possible to compute multiple or clustered eigenvalues with a single vector restart method. A block method is not required. These schemes are analyzed with respect to numerical stability, and computational results are presented.
引用
收藏
页码:789 / 821
页数:33
相关论文
共 42 条
[1]  
[Anonymous], 237 U MANCH
[3]  
BAI ZJ, 1993, LINEAR ALGEBRA APPL, V186, P73
[4]   ALGORITHM - SOLUTION OF MATRIX EQUATION AX+XB = C [J].
BARTELS, RH ;
STEWART, GW .
COMMUNICATIONS OF THE ACM, 1972, 15 (09) :820-&
[5]  
BENNANJI M, 1993, TRPA9325 CERFACS
[6]  
BOJANCZYK A., 1993, NATO ASI SERIES E, P351
[7]  
Braconnier T., 1993, TRPA9325 CERFACS
[8]  
Calvetti D., 1994, ELECTRON T NUMER ANA, V2, P1
[9]  
Chatelin F., 1993, EIGENVALUES MATRICES
[10]  
CHATELIN F, 1993, QUALITATIVE COMPUTIN