Power method with a dominant eigenvalue cluster (source code)

= Power method with a dominant eigenvalue cluster

When several eigenvalues share the largest modulus, normalized power iterates need not converge individually, but every limit lies in their joint invariant subspace.