version 1.7, 2010/12/21 13:53:36
|
version 1.8, 2011/07/22 07:38:09
|
Line 1
|
Line 1
|
SUBROUTINE DPOTF2( UPLO, N, A, LDA, INFO ) |
SUBROUTINE DPOTF2( UPLO, N, A, LDA, INFO ) |
* |
* |
* -- LAPACK routine (version 3.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..-- |
* November 2006 |
* -- April 2011 -- |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER UPLO |
CHARACTER UPLO |
Line 20
|
Line 20
|
* positive definite matrix A. |
* positive definite matrix A. |
* |
* |
* The factorization has the form |
* The factorization has the form |
* A = U' * U , if UPLO = 'U', or |
* A = U**T * U , if UPLO = 'U', or |
* A = L * L', if UPLO = 'L', |
* A = L * L**T, if UPLO = 'L', |
* where U is an upper triangular matrix and L is lower triangular. |
* where U is an upper triangular matrix and L is lower triangular. |
* |
* |
* This is the unblocked version of the algorithm, calling Level 2 BLAS. |
* This is the unblocked version of the algorithm, calling Level 2 BLAS. |
Line 48
|
Line 48
|
* triangular part of A is not referenced. |
* triangular part of A is not referenced. |
* |
* |
* On exit, if INFO = 0, the factor U or L from the Cholesky |
* On exit, if INFO = 0, the factor U or L from the Cholesky |
* factorization A = U'*U or A = L*L'. |
* factorization A = U**T *U or A = L*L**T. |
* |
* |
* LDA (input) INTEGER |
* LDA (input) INTEGER |
* The leading dimension of the array A. LDA >= max(1,N). |
* The leading dimension of the array A. LDA >= max(1,N). |
Line 107
|
Line 107
|
* |
* |
IF( UPPER ) THEN |
IF( UPPER ) THEN |
* |
* |
* Compute the Cholesky factorization A = U'*U. |
* Compute the Cholesky factorization A = U**T *U. |
* |
* |
DO 10 J = 1, N |
DO 10 J = 1, N |
* |
* |
Line 131
|
Line 131
|
10 CONTINUE |
10 CONTINUE |
ELSE |
ELSE |
* |
* |
* Compute the Cholesky factorization A = L*L'. |
* Compute the Cholesky factorization A = L*L**T. |
* |
* |
DO 20 J = 1, N |
DO 20 J = 1, N |
* |
* |