Power iteration eigenpair diagnostics (${\rho }_k=x_k^{\ast }A{x}_k,r_k=A{x}_k-{\rho }_k{x}_k$)

For each unit vector iterate compute the Rayleigh quotient ${\rho }_k=x_k^{\ast }A{x}_k$ and the eigenpair residual $r_k=A{x}_k-{\rho }_k{x}_k$. Use the same multiplication $A{x}_k$ for both diagnostics and the next normalized iterate. For any complex matrix, this Rayleigh quotient minimizes $\Vert A{x}_k-\mu x_k{\Vert }_2$ over complex scalars $\mu$, because $x_k^{\ast }{r}_k=0$. Under power iteration convergence for a diagonalizable matrix, ${\rho }_k\to {\lambda }_1$ and $\Vert r_k{\Vert }_2\to 0$, even when the vectors alternate signs or phases.^[ref]

A dimensionless stopping quantity iswhen the denominator is positive; $\Vert A{\Vert }_2$ is the operator norm, and a known upper bound may replace it if computing it is costly. A small eigenpair residual certifies a small backward error of an approximate eigenpair, but does not establish that the eigenvalue found is dominant.

For a Hermitian matrix with orthonormal eigenbasis $u_j$, write $x_k={\sum }_j{a}_j{u}_j$. Then the independently useful identityshows that some eigenvalue lies within $\Vert r_k{\Vert }_2$ of ${\rho }_k$. If a specific target eigenvalue ${\lambda }_1$ is separated from all the others by ${\delta }_k={\min }_{j>1}|{\lambda }_j-{\rho }_k|>0$, it also gives ${\epsilon }_k\le \Vert r_k{\Vert }_2/{\delta }_k$. This separation is measured from the computed Rayleigh quotient; it differs from the dominant-modulus spectral gap. For a nonnormal matrix, the same orthogonal expansion is unavailable and a small eigenpair residual alone need not imply a comparably small distance to the spectrum.