version 1.3, 2010/08/13 21:03:56
|
version 1.5, 2011/07/22 07:38:10
|
Line 21
|
Line 21
|
* pivoting of a real symmetric positive semidefinite matrix A. |
* pivoting of a real symmetric positive semidefinite matrix A. |
* |
* |
* The factorization has the form |
* The factorization has the form |
* P' * A * P = U' * U , if UPLO = 'U', |
* P**T * A * P = U**T * U , if UPLO = 'U', |
* P' * A * P = L * L', if UPLO = 'L', |
* P**T * A * P = L * L**T, if UPLO = 'L', |
* where U is an upper triangular matrix and L is lower triangular, and |
* where U is an upper triangular matrix and L is lower triangular, and |
* P is stored as vector PIV. |
* P is stored as vector PIV. |
* |
* |
Line 170
|
Line 170
|
* |
* |
IF( UPPER ) THEN |
IF( UPPER ) THEN |
* |
* |
* Compute the Cholesky factorization P' * A * P = U' * U |
* Compute the Cholesky factorization P**T * A * P = U**T * U |
* |
* |
DO 140 K = 1, N, NB |
DO 140 K = 1, N, NB |
* |
* |
Line 257
|
Line 257
|
* |
* |
ELSE |
ELSE |
* |
* |
* Compute the Cholesky factorization P' * A * P = L * L' |
* Compute the Cholesky factorization P**T * A * P = L * L**T |
* |
* |
DO 180 K = 1, N, NB |
DO 180 K = 1, N, NB |
* |
* |