An inverse free preconditioned Krylov subspace method for symmetric generalized eigenvalue problems

被引:89
作者
Golub, GH [1 ]
Ye, Q
机构
[1] Stanford Univ, Sci Comp & Computat Math Program, Dept Comp Sci, Stanford, CA 94305 USA
[2] Univ Kentucky, Dept Math, Lexington, KY 40506 USA
关键词
Krylov subspace; preconditioning; eigenvalue problems;
D O I
10.1137/S1064827500382579
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
In this paper, we present an inverse free Krylov subspace method for finding some extreme eigenvalues of the symmetric definite generalized eigenvalue problem Ax = lambdaBx. The basic method takes a form of inner-outer iterations and involves no inversion of B or any shift-and-invert matrix A - lambda(0)B. A convergence analysis is presented that leads to a preconditioning scheme for accelerating convergence through some equivalent transformations of the eigenvalue problem. Numerical examples are given to illustrate the convergence properties and to demonstrate the competitiveness of the method.
引用
收藏
页码:312 / 334
页数:23
相关论文
共 42 条
[1]  
[Anonymous], ELECT T NUMER ANAL
[2]  
Bai Zhaojun, 2000, TEMPLATES SOLUTION A
[3]  
Bergamaschi L, 2000, NUMER LINEAR ALGEBR, V7, P99, DOI 10.1002/(SICI)1099-1506(200004/05)7:3<99::AID-NLA188>3.0.CO
[4]  
2-5
[5]  
BOURAS A, 2000, TRPA0016 CERFACS
[6]   A subspace preconditioning algorithm for eigenvector/eigenvalue computation [J].
Bramble, JH ;
Pasciak, JE ;
Knyazev, AV .
ADVANCES IN COMPUTATIONAL MATHEMATICS, 1996, 6 (02) :159-189
[7]   THE DAVIDSON METHOD [J].
CROUZEIX, M ;
PHILIPPE, B ;
SADKANE, M .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 1994, 15 (01) :62-76
[8]  
Demmel J.W., 1997, APPL NUMERICAL LINEA
[9]   MINIMIZATION OF THE COMPUTATIONAL LABOR IN DETERMINING THE 1ST EIGENVALUES OF DIFFERENTIAL-OPERATORS [J].
DYAKONOV, EG ;
OREKHOV, MY .
MATHEMATICAL NOTES, 1980, 27 (5-6) :382-391
[10]   Jacobi-Davidson style QR and QZ algorithms for the reduction of matrix pencils [J].
Fokkema, DR ;
Sleijpen, GLG ;
Van der Vorst, HA .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 1998, 20 (01) :94-125