version 1.7, 2010/12/21 13:53:53
|
version 1.8, 2011/07/22 07:38:19
|
Line 2
|
Line 2
|
$ S, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, |
$ S, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, |
$ RWORK, INFO ) |
$ 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 116
|
Line 116
|
* AF (input or output) COMPLEX*16 array, dimension (LDAF,N) |
* AF (input or output) COMPLEX*16 array, dimension (LDAF,N) |
* If FACT = 'F', then AF is an input argument and on entry |
* If FACT = 'F', then AF 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, in the same storage |
* factorization A = U**H *U or A = L*L**H, in the same storage |
* format as A. If EQUED .ne. 'N', then AF is the factored form |
* format as A. If EQUED .ne. 'N', then AF is the factored form |
* of the equilibrated matrix diag(S)*A*diag(S). |
* of the equilibrated matrix diag(S)*A*diag(S). |
* |
* |
* If FACT = 'N', then AF is an output argument and on exit |
* If FACT = 'N', then AF 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 original |
* factorization A = U**H *U or A = L*L**H of the original |
* matrix A. |
* matrix A. |
* |
* |
* If FACT = 'E', then AF is an output argument and on exit |
* If FACT = 'E', then AF 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 319
|
Line 319
|
* |
* |
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. |
* |
* |
CALL ZLACPY( UPLO, N, N, A, LDA, AF, LDAF ) |
CALL ZLACPY( UPLO, N, N, A, LDA, AF, LDAF ) |
CALL ZPOTRF( UPLO, N, AF, LDAF, INFO ) |
CALL ZPOTRF( UPLO, N, AF, LDAF, INFO ) |