version 1.7, 2010/12/21 13:53:33
|
version 1.8, 2011/07/22 07:38:08
|
Line 22
|
Line 22
|
* (upper or lower) bidiagonal matrix with diagonal D and offdiagonal |
* (upper or lower) bidiagonal matrix with diagonal D and offdiagonal |
* E, accumulating the transformations if desired. Letting B denote |
* E, accumulating the transformations if desired. Letting B denote |
* the input bidiagonal matrix, the algorithm computes orthogonal |
* the input bidiagonal matrix, the algorithm computes orthogonal |
* matrices Q and P such that B = Q * S * P' (P' denotes the transpose |
* matrices Q and P such that B = Q * S * P**T (P**T denotes the transpose |
* of P). The singular values S are overwritten on D. |
* of P). The singular values S are overwritten on D. |
* |
* |
* The input matrix U is changed to U * Q if desired. |
* The input matrix U is changed to U * Q if desired. |
* The input matrix VT is changed to P' * VT if desired. |
* The input matrix VT is changed to P**T * VT if desired. |
* The input matrix C is changed to Q' * C if desired. |
* The input matrix C is changed to Q**T * C if desired. |
* |
* |
* See "Computing Small Singular Values of Bidiagonal Matrices With |
* See "Computing Small Singular Values of Bidiagonal Matrices With |
* Guaranteed High Relative Accuracy," by J. Demmel and W. Kahan, |
* Guaranteed High Relative Accuracy," by J. Demmel and W. Kahan, |
Line 84
|
Line 84
|
* |
* |
* VT (input/output) DOUBLE PRECISION array, dimension (LDVT, NCVT) |
* VT (input/output) DOUBLE PRECISION array, dimension (LDVT, NCVT) |
* On entry, contains a matrix which on exit has been |
* On entry, contains a matrix which on exit has been |
* premultiplied by P', dimension N-by-NCVT if SQRE = 0 |
* premultiplied by P**T, dimension N-by-NCVT if SQRE = 0 |
* and (N+1)-by-NCVT if SQRE = 1 (not referenced if NCVT=0). |
* and (N+1)-by-NCVT if SQRE = 1 (not referenced if NCVT=0). |
* |
* |
* LDVT (input) INTEGER |
* LDVT (input) INTEGER |
Line 104
|
Line 104
|
* |
* |
* C (input/output) DOUBLE PRECISION array, dimension (LDC, NCC) |
* C (input/output) DOUBLE PRECISION array, dimension (LDC, NCC) |
* On entry, contains an N-by-NCC matrix which on exit |
* On entry, contains an N-by-NCC matrix which on exit |
* has been premultiplied by Q' dimension N-by-NCC if SQRE = 0 |
* has been premultiplied by Q**T dimension N-by-NCC if SQRE = 0 |
* and (N+1)-by-NCC if SQRE = 1 (not referenced if NCC=0). |
* and (N+1)-by-NCC if SQRE = 1 (not referenced if NCC=0). |
* |
* |
* LDC (input) INTEGER |
* LDC (input) INTEGER |