version 1.8, 2010/12/21 13:53:24
|
version 1.9, 2011/07/22 07:38:04
|
Line 1
|
Line 1
|
SUBROUTINE DBDSDC( UPLO, COMPQ, N, D, E, U, LDU, VT, LDVT, Q, IQ, |
SUBROUTINE DBDSDC( UPLO, COMPQ, N, D, E, U, LDU, VT, LDVT, Q, IQ, |
$ WORK, IWORK, INFO ) |
$ WORK, IWORK, INFO ) |
* |
* |
* -- LAPACK routine (version 3.2.2) -- |
* -- LAPACK routine (version 3.3.1) -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* June 2010 |
* -- April 2011 -- |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER COMPQ, UPLO |
CHARACTER COMPQ, UPLO |
Line 77
|
Line 77
|
* |
* |
* VT (output) DOUBLE PRECISION array, dimension (LDVT,N) |
* VT (output) DOUBLE PRECISION array, dimension (LDVT,N) |
* If COMPQ = 'I', then: |
* If COMPQ = 'I', then: |
* On exit, if INFO = 0, VT' contains the right singular |
* On exit, if INFO = 0, VT**T contains the right singular |
* vectors of the bidiagonal matrix. |
* vectors of the bidiagonal matrix. |
* For other values of COMPQ, VT is not referenced. |
* For other values of COMPQ, VT is not referenced. |
* |
* |
Line 287
|
Line 287
|
CALL DLASCL( 'G', 0, 0, ORGNRM, ONE, N, 1, D, N, IERR ) |
CALL DLASCL( 'G', 0, 0, ORGNRM, ONE, N, 1, D, N, IERR ) |
CALL DLASCL( 'G', 0, 0, ORGNRM, ONE, NM1, 1, E, NM1, IERR ) |
CALL DLASCL( 'G', 0, 0, ORGNRM, ONE, NM1, 1, E, NM1, IERR ) |
* |
* |
EPS = DLAMCH( 'Epsilon' ) |
EPS = (0.9D+0)*DLAMCH( 'Epsilon' ) |
* |
* |
MLVL = INT( LOG( DBLE( N ) / DBLE( SMLSIZ+1 ) ) / LOG( TWO ) ) + 1 |
MLVL = INT( LOG( DBLE( N ) / DBLE( SMLSIZ+1 ) ) / LOG( TWO ) ) + 1 |
SMLSZP = SMLSIZ + 1 |
SMLSZP = SMLSIZ + 1 |