Abstract:We introduce a new Projected Rayleigh Quotient Iteration aimed at improving the convergence behaviour of classic Rayleigh Quotient iteration (RQI) by incorporating approximate information about the target eigenvector at each step. While classic RQI exhibits local cubic convergence for Hermitian matrices, its global behaviour can be unpredictable, whereby it may converge to an eigenvalue far away from the target, even when started with accurate initial conditions. This problem is exacerbated when the eigenvalues are closely spaced. The key idea of the new algorithm is at each step to add a complex-valued projection to the original matrix (that depends on the current eigenvector approximation), such that the unwanted eigenvalues are lifted into the complex plane while the target stays close to the real line, thereby increasing the spacing between the target eigenvalue and the rest of the spectrum. Making better use of the eigenvector approximation leads to more robust convergence behaviour and the new method converges reliably to the correct target eigenpair for a significantly wider range of initial vectors than does classic RQI. We prove that the method converges locally cubically and we present several numerical examples demonstrating the improved global convergence behaviour. In particular, we apply it to compute eigenvalues in a band-gap spectrum of a Sturm-Liouville operator used to model photonic crystal fibres, where the target and unwanted eigenvalues are closely spaced. The examples show that the new method converges to the desired eigenpair even when the eigenvalue spacing is very small, often succeeding when classic RQI fails.

Rayleigh-Ritz Approximation For The Linear Response Eigenvalue Problem

An Analysis of the Rayleigh-Ritz Method for Approximating Eigenspaces

Convergence of Q-Ritz Pairs, Refined Q-Ritz Vectors and Q-Rayleigh-ritz Method for Quadratic Eigenvalue Problems

On the Convergence of Ritz Pairs and Refined Ritz Vectors for Quadratic Eigenvalue Problems

Rayleigh-Ritz majorization error bounds with applications to FEM

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

The Rayleigh–Ritz method, refinement and Arnoldi process for periodic matrix pairs

Resolvent sampling based Rayleigh-Ritz method for large-scale nonlinear eigenvalue problems

On the Rayleigh-Ritz variational method

Automatic rational approximation and linearization of nonlinear eigenvalue problems

Approximate Cayley Transform Methods for Inverse Eigenvalue Problems and Convergence Analysis

On Convergence of the Inexact Rayleigh Quotient Iteration with the Lanczos Method Used for Solving Linear Systems

An Inexact Cayley Transform Method for Inverse Eigenvalue Problems with Multiple Eigenvalues

A complex-projected Rayleigh quotient iteration for targeting interior eigenvalues

2D Eigenvalue Problem III: Convergence Analysis of the 2D Rayleigh Quotient Iteration

A Refined Arnoldi Type Method for Large Scale Eigenvalue Problems

The Lower/Upper Bound Property of Approximate Eigenvalues by Nonconforming Finite Element Methods for Elliptic Operators

A Convergence Analysis of the Inexact Rayleigh Quotient Iteration and Simplified Jacobi-Davidson Method for the Large Hermitian Matrix Eigenproblem

Solving a Class of Infinite-Dimensional Tensor Eigenvalue Problems by Translational Invariant Tensor Ring Approximations

Immersed finite element method for eigenvalue problem

Theory and numerics of subspace approximation of eigenvalue problems