R. Byers et al., THE MATRIX SIGN FUNCTION-METHOD AND THE COMPUTATION OF INVARIANT SUBSPACES, SIAM journal on matrix analysis and applications, 18(3), 1997, pp. 615-632
A perturbation analysis shows that if a numerically stable procedure i
s used to compute the matrix sign function, then it is competitive wit
h conventional methods for computing invariant subspaces. Stability an
alysis of the Newton iteration improves an earlier result of Byers and
confirms that ill-conditioned iterates may cause numerical instabilit
y. Numerical examples demonstrate the theoretical results.