AN ANALYSIS OF THE RAYLEIGH--RITZ AND REFINED RAYLEIGH--RITZ METHODS FOR REGULAR NONLINEAR EIGENVALUE PROBLEMS

We establish a general convergence theory of the Rayleigh--Ritz method and the refined Rayleigh--Ritz method for computing some simple eigenpair (\lambda \ast , x\ast ) of a given analytic regular nonlinear eigenvalue problem (NEP). In terms of the deviation \varepsilon of x\ast from a given subspace W, we establish a priori convergence results on the Ritz value, the Ritz vector, and the refined Ritz vector. The results show that, as \varepsilon ! 0, there exists a Ritz value that unconditionally converges to \lambda \ast , as does the corresponding refined Ritz vector, but the Ritz vector converges conditionally and may fail to converge and even may not be unique. We also present an error bound for the approximate eigenvector in terms of the computable residual norm of a given approximate eigenpair and give lower and upper bounds for the error of the refined Ritz vector and the Ritz vector as well as for that of the corresponding residual norms. These results nontrivially extend some convergence results on these two methods for the linear eigenvalue problem to the NEP. Examples are constructed to illustrate the main results.