Otherwise, it is the storage work area used by this subroutine. Its size is specified by naux . The work area is deallocated before control is returned to the calling program. You have the option of having the minimum required value for naux dynamically returned to your program. I A is an m by n real general matrix.
V is a real general orthogonal matrix of order n .
On output, V overwrites the first n rows and n columns of A . U T B is an n by nb real general matrix. If m or n is equal to 0, no computation is performed. Reduce the real general matrix A to bidiagonal form using Householder transformations. Iteratively reduce the bidiagonal form to diagonal form using a variant of the QR algorithm.