Matlab, Algorithms, Eigenvalue, Singular Values, Computation.
In all three cases, the QR iteration itself is preceded by a reduction to a compact form. All transformations are orthogonal similarities using Givens and Householder transformations. These are numerically stable, preserve eigenvalues, and preserve any symmetry. The starting matrix for all three