The Convergence of Harmonic Ritz Values, Harmonic Ritz Vectors, and Refined Harmonic Ritz Vectors

ZX Jia
DOI: https://doi.org/10.1090/s0025-5718-04-01684-9
IF: 2.118
2005-01-01
Mathematics of Computation
Abstract:This paper concerns the harmonic projection methods for computing an approximation to an eigenpair (Λ, of a large matrix . Given a target point τ and a subspace W that contains an approximation to i, a harmonic projection method returns an approximation to Three convergence results are established as the deviation c of from W approaches zero. First, the harmonic Ritz value converges to Λ if certain Rayleigh quotient is uniformly nonsingular. Second, the harmonic Ritz vector converges to if the Rayleigh quotient is uniformly nonsingular and remains well separated from the other harmonic Ritz values. Third, better error bounds for the convergence of are derived when converges. However, we show that the harmonic projection method can fail to find the desired eigenvalue Λ or in other words they can miss Λ if it is very close to r. We then propose to compute the Rayleigh quotient of with respect to and take it as a new approximate eigenvalue, which is shown to converge to Λ no matter how r is close to Λ whenever tends to Finally, we deduce from a known result that under the assumption that the Rayleigh quotient is uniformly nonsingular certain refined harmonic Ritz vector or more generally refined eigenvector approximation introduced by the author converge. We construct examples to illustrate some of our theory.
What problem does this paper attempt to address?