数值计算与计算机应用
數值計算與計算機應用
수치계산여계산궤응용
JOURNAL ON NUMERICAL METHODS AND COMPUTER APPLICATIONS
2003年
2期
101-110
,共10页
Ritz pairs%Refined Ritz pairs%the Arnoldi process%the Refined Arnoldi method%Rayleigh quotient%Approximate eigenvalue
The refined Arnoldi method proposed by Jia is used for computing some eigen-pairs of large matrices. In contrast to the Arnoldi method, the fundamental dif-ference is that the refined method seeks certain refined Ritz vectors, which aredifferent from the Ritz vectors obtained by the Arnoldi method, from a projection space with minimal residuals to approximate the desired eigenvectors. In com-parison with the Ritz vectors, the refined Ritz vectors are guaranteed to converge theoretically and can converge much faster numerically. In this paper we propose to replace the Ritz values, obtained by the Arnoldi method with respect to a Krylovsubspace, by the ones obtained with respect to the subspace spanned by the refined Ritz vectors. We discuss how to compute these new approximations cheaply and reliably. Theoretical error bounds between the original Ritz values and the new Ritz values are established. Finally, we present a variant of the refined Arnoldi al-gorithm for an augmented Krylov subspace and discuss restarting issue. Numerical results confirm efficiency of the new algorithm.