Numerical Considerations in Computing Invariant Subspaces