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.

2D Eigenvalue Problem II: Rayleigh Quotient Iteration and Applications

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

On Rayleigh Quotient Iteration for Dual Quaternion Hermitian Eigenvalue Problem

Rayleigh Quotient Iteration, cubic convergence, and second covariant derivative

RAYLEIGH QUOTIENT ITERATION FOR COMPUTING TENSOR EIGENPAIRS

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

A Note on the Approximation of Two-Sided Rayleigh Quotient to Eigenvalues

A complex-projected Rayleigh quotient iteration for targeting interior eigenvalues

2D Eigenvalue Problems I: Existence and Number of Solutions

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

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

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

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

A Newton‐type Method for Two‐dimensional Eigenvalue Problems

A Parallelized Contour Integral Rayleigh–Ritz Method for Computing Critical Eigenvalues of Large-Scale Power Systems

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

On Convergence of the Inexact Rayleigh Quotient Iteration with MINRES.

Eigenvalue backward errors of Rosenbrock systems and optimization of sums of Rayleigh quotient

The Rayleigh Quotient Iteration Method For Computing Eigenvalue Bounds Of Structures With Bounded Uncertain Parameters

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

An Improved Two‐Step Method for Inverse Eigenvalue Problems