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.

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

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

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

An Analysis of the Rayleigh-Ritz Method for Approximating Eigenspaces

Composite orthogonal projection methods for large matrix eigenproblems

Some Theoretical Comparisons of Refined Ritz Vectors and Ritz Vectors

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

A Theoretical Comparison of Two Classes of Projection Methods

An Ulm-like Cayley Transform Method for Inverse Eigenvalue Problems with Multiple Eigenvalues

Residuals of refined projection methods for large matrix eigenproblems

A Refined Subspace Iteration Algorithm For Large Sparse Eigenproblems

Accelerating Convergence by Augmented Rayleigh-Ritz Projections for Large-Scale Eigenpair Computation

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

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

On Rayleigh Quotient Iteration for Dual Quaternion Hermitian Eigenvalue Problem

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

Rayleigh-Ritz Approximation For The Linear Response Eigenvalue Problem

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

Convergence Rate of the Relaxed CQ Algorithm under Hölderian Type Error Bound Property

A complex-projected Rayleigh quotient iteration for targeting interior eigenvalues

2D Eigenvalue Problem II: Rayleigh Quotient Iteration and Applications