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 Rayleigh Quotient Iteration for Dual Quaternion Hermitian Eigenvalue Problem

Generalized Jacobi Method for Computing Eigenvalues of Dual Quaternion Hermitian Matrices

A Power Method for Computing the Dominant Eigenvalue of a Dual Quaternion Hermitian Matrix

Applications of Dual Complex Adjoint Matrix in Eigenvalue Computation of Dual Quaternion Hermitian Matrix

2D Eigenvalue Problem II: Rayleigh Quotient Iteration and Applications

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

Jacobi method for dual quaternion Hermitian eigenvalue problems and applications

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

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

RAYLEIGH QUOTIENT ITERATION FOR COMPUTING TENSOR EIGENPAIRS

A complex-projected Rayleigh quotient iteration for targeting interior eigenvalues

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

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

Eigenvalues of Dual Hermitian Matrices with Application in Formation Control

Quasi-determinant and right eigenvalues of dual quaternion matrices

Dual Complex Adjoint Matrix: Applications in Dual Quaternion Research

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

Rayleigh Quotient Iteration, cubic convergence, and second covariant derivative

Dual quaternion QR decompositon and its corresponding complex structure-preserving algorithms

Fast parallelizable methods for the Hermitian eigenvalue problem