The first step can be done using Householder reflections for a cost of 4mn2 − 4n3/3 flops, assuming that only the singular values are needed and not the singular vectors.
The second step can be done by a variant of the QR algorithm for the computation of eigenvalues, which was first described by .