A new stable bidiagonal reduction algorithm

The singular values of B are those of X perturbed by no more than f(m, n)εM∥X∥F where f(m, n) is a modestly growing function and εM is the machine unit. Under certain assumptions, relative error bounds on the singular values are possible.