Two-sided Grassmann-Rayleigh quotient iteration

From MaRDI portal
Publication:2267769



Abstract: The two-sided Rayleigh quotient iteration proposed by Ostrowski computes a pair of corresponding left-right eigenvectors of a matrix C. We propose a Grassmannian version of this iteration, i.e., its iterates are pairs of p-dimensional subspaces instead of one-dimensional subspaces in the classical case. The new iteration generically converges locally cubically to the pairs of left-right p-dimensional invariant subspaces of C. Moreover, Grassmannian versions of the Rayleigh quotient iteration are given for the generalized Hermitian eigenproblem, the Hamiltonian eigenproblem and the skew-Hamiltonian eigenproblem.


The authors merge Ostrowski's two-sided Rayleigh quotient iteration (RQI) for a single-vector pair with the block (or Grassmann) RQI to obtain a two-sided Grassmann iteration (2sGRQI) that operates on pairs of \(p\)-dimensional subspaces. It is shown that this iteration converges locally cubically to the pair of left-right \(p\)-dimensional eigenspaces of a nonnormal matrix \(C\). Several lemmas are proved on their own right. The practical implementation of 2sGRQI is a non trivial issue because of the need to solve Sylvester equations that become increasingly ill-conditioned as the iteration progresses. This inherent ill-condition confers to 2sGRQI its superlinear convergence. It is shown that the 2sGRQI algorithm induces particular one-sided formulations for some structured eigenproblems. It is observed numerically that certain ways of solving the Sylvester equations are preferable to others.



Cites work



Describes a project that uses

Uses Software






This page was built for publication: Two-sided Grassmann-Rayleigh quotient iteration

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2267769)