version 1.17, 2018/05/29 07:18:09
|
version 1.18, 2020/05/21 21:46:02
|
Line 39
|
Line 39
|
*> the Bunch-Kaufman diagonal pivoting method. The form of the |
*> the Bunch-Kaufman diagonal pivoting method. The form of the |
*> factorization is |
*> factorization is |
*> |
*> |
*> A = U*D*U**T or A = L*D*L**T |
*> A = U**T*D*U or A = L*D*L**T |
*> |
*> |
*> where U (or L) is a product of permutation and unit upper (lower) |
*> where U (or L) is a product of permutation and unit upper (lower) |
*> triangular matrices, and D is symmetric and block diagonal with |
*> triangular matrices, and D is symmetric and block diagonal with |
Line 144
|
Line 144
|
*> |
*> |
*> \verbatim |
*> \verbatim |
*> |
*> |
*> If UPLO = 'U', then A = U*D*U**T, where |
*> If UPLO = 'U', then A = U**T*D*U, where |
*> U = P(n)*U(n)* ... *P(k)U(k)* ..., |
*> U = P(n)*U(n)* ... *P(k)U(k)* ..., |
*> i.e., U is a product of terms P(k)*U(k), where k decreases from n to |
*> i.e., U is a product of terms P(k)*U(k), where k decreases from n to |
*> 1 in steps of 1 or 2, and D is a block diagonal matrix with 1-by-1 |
*> 1 in steps of 1 or 2, and D is a block diagonal matrix with 1-by-1 |
Line 262
|
Line 262
|
* |
* |
IF( UPPER ) THEN |
IF( UPPER ) THEN |
* |
* |
* Factorize A as U*D*U**T using the upper triangle of A |
* Factorize A as U**T*D*U using the upper triangle of A |
* |
* |
* K is the main loop index, decreasing from N to 1 in steps of |
* K is the main loop index, decreasing from N to 1 in steps of |
* KB, where KB is the number of columns factorized by DLASYF; |
* KB, where KB is the number of columns factorized by DLASYF; |