Bai Z, Stewart G.W.  1997.  Algorithm 776: SRRIT: a Fortran subroutine to calculate the dominant invariant subspace of a nonsymmetric matrix. ACM Trans. Math. Softw.. 23(4):494-513.