Augmented implicitly restarted Lanczos bidiagonalization methods

被引:205
作者
Baglama, J [1 ]
Reichel, L
机构
[1] Univ Rhode Isl, Dept Math, Kingston, RI 02881 USA
[2] Kent State Univ, Dept Math Sci, Kent, OH 44242 USA
关键词
singular value computation; partial singular value decomposition; iterative method; large-scale computation;
D O I
10.1137/04060593X
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
New restarted Lanczos bidiagonalization methods for the computation of a few of the largest or smallest singular values of a large matrix are presented. Restarting is carried out by augmentation of Krylov subspaces that arise naturally in the standard Lanczos bidiagonalization method. The augmenting vectors are associated with certain Ritz or harmonic Ritz vectors. Computed examples show the new methods to be competitive with available schemes.
引用
收藏
页码:19 / 42
页数:24
相关论文
共 33 条
[1]  
[Anonymous], 1998, ARPACK USERS GUIDE S
[2]  
[Anonymous], 1961, NUMER MATH, DOI DOI 10.1007/BF01386022
[3]   Iterative methods for the computation of a few eigenvalues of a large symmetric matrix [J].
Baglama, J ;
Calvetti, D ;
Reichel, L .
BIT NUMERICAL MATHEMATICS, 1996, 36 (03) :400-421
[4]   Algorithm 827: irbleigs: a MATLAB program for computing a few eigenpairs of a large sparse Hermitian matrix [J].
Baglama, J ;
Calvetti, D ;
Reichel, L .
ACM TRANSACTIONS ON MATHEMATICAL SOFTWARE, 2003, 29 (03) :337-348
[5]   IRBL: An implicitly restarted block-lanczos method for large-scale Hermitian eigenproblems [J].
Baglama, J ;
Calvetti, D ;
Reichel, L .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2003, 24 (05) :1650-1677
[6]  
BJORCK A, 1994, BIT, V34, P510, DOI 10.1007/BF01934265
[7]  
BJORCK A, 1996, NUMERICAL METHODS LE
[8]   UPDATING SINGULAR VALUE DECOMPOSITION [J].
BUNCH, JR ;
NIELSEN, CP .
NUMERISCHE MATHEMATIK, 1978, 31 (02) :111-129
[9]  
Calvetti D., 1994, ELECTRON T NUMER ANA, V2, P1
[10]   TRACKING A FEW EXTREME SINGULAR-VALUES AND VECTORS IN SIGNAL-PROCESSING [J].
COMON, P ;
GOLUB, GH .
PROCEEDINGS OF THE IEEE, 1990, 78 (08) :1327-1343