version 1.16, 2018/05/29 06:55:20
|
version 1.18, 2020/05/21 21:46:01
|
Line 106
|
Line 106
|
*> triangular part of the matrix A, and the strictly upper |
*> triangular part of the matrix A, and the strictly upper |
*> triangular part of A is not referenced. |
*> triangular part of A is not referenced. |
*> On exit, if iterative refinement has been successfully used |
*> On exit, if iterative refinement has been successfully used |
*> (INFO.EQ.0 and ITER.GE.0, see description below), then A is |
*> (INFO = 0 and ITER >= 0, see description below), then A is |
*> unchanged, if double precision factorization has been used |
*> unchanged, if double precision factorization has been used |
*> (INFO.EQ.0 and ITER.LT.0, see description below), then the |
*> (INFO = 0 and ITER < 0, see description below), then the |
*> array A contains the factor U or L from the Cholesky |
*> array A contains the factor U or L from the Cholesky |
*> factorization A = U**T*U or A = L*L**T. |
*> factorization A = U**T*U or A = L*L**T. |
*> \endverbatim |
*> \endverbatim |
Line 413
|
Line 413
|
30 CONTINUE |
30 CONTINUE |
* |
* |
* If we are at this place of the code, this is because we have |
* If we are at this place of the code, this is because we have |
* performed ITER=ITERMAX iterations and never satisified the |
* performed ITER=ITERMAX iterations and never satisfied the |
* stopping criterion, set up the ITER flag accordingly and follow |
* stopping criterion, set up the ITER flag accordingly and follow |
* up on double precision routine. |
* up on double precision routine. |
* |
* |