version 1.7, 2010/12/21 13:53:53
|
version 1.8, 2011/07/22 07:38:18
|
Line 2
|
Line 2
|
$ EQUED, S, B, LDB, X, LDX, RCOND, FERR, BERR, |
$ EQUED, S, B, LDB, X, LDX, RCOND, FERR, BERR, |
$ WORK, RWORK, INFO ) |
$ WORK, RWORK, INFO ) |
* |
* |
* -- LAPACK driver routine (version 3.2) -- |
* -- LAPACK driver 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..-- |
* November 2006 |
* -- April 2011 -- |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER EQUED, FACT, UPLO |
CHARACTER EQUED, FACT, UPLO |
Line 118
|
Line 118
|
* AFB (input or output) COMPLEX*16 array, dimension (LDAFB,N) |
* AFB (input or output) COMPLEX*16 array, dimension (LDAFB,N) |
* If FACT = 'F', then AFB is an input argument and on entry |
* If FACT = 'F', then AFB is an input argument and on entry |
* contains the triangular factor U or L from the Cholesky |
* contains the triangular factor U or L from the Cholesky |
* factorization A = U**H*U or A = L*L**H of the band matrix |
* factorization A = U**H *U or A = L*L**H of the band matrix |
* A, in the same storage format as A (see AB). If EQUED = 'Y', |
* A, in the same storage format as A (see AB). If EQUED = 'Y', |
* then AFB is the factored form of the equilibrated matrix A. |
* then AFB is the factored form of the equilibrated matrix A. |
* |
* |
* If FACT = 'N', then AFB is an output argument and on exit |
* If FACT = 'N', then AFB is an output argument and on exit |
* returns the triangular factor U or L from the Cholesky |
* returns the triangular factor U or L from the Cholesky |
* factorization A = U**H*U or A = L*L**H. |
* factorization A = U**H *U or A = L*L**H. |
* |
* |
* If FACT = 'E', then AFB is an output argument and on exit |
* If FACT = 'E', then AFB is an output argument and on exit |
* returns the triangular factor U or L from the Cholesky |
* returns the triangular factor U or L from the Cholesky |
* factorization A = U**H*U or A = L*L**H of the equilibrated |
* factorization A = U**H *U or A = L*L**H of the equilibrated |
* matrix A (see the description of A for the form of the |
* matrix A (see the description of A for the form of the |
* equilibrated matrix). |
* equilibrated matrix). |
* |
* |
Line 351
|
Line 351
|
* |
* |
IF( NOFACT .OR. EQUIL ) THEN |
IF( NOFACT .OR. EQUIL ) THEN |
* |
* |
* Compute the Cholesky factorization A = U'*U or A = L*L'. |
* Compute the Cholesky factorization A = U**H *U or A = L*L**H. |
* |
* |
IF( UPPER ) THEN |
IF( UPPER ) THEN |
DO 40 J = 1, N |
DO 40 J = 1, N |