title: Which algorithm does DGESVD in LAPACK implement?
body: <p>I'm curious about the DGESVD function used to calculate SVD in MATLAB. As far as I can tell from "Matrix Computations" by Gene H. Golub and Charles F. Van Loan, there are two possible bidiagonalization schemes used - the Householder Bidiagonalization and R-Bidiagonalization. However, I wasn't able to determine from the LAPACK documentation which is being used. Anyone have any ideas?</p>