version 1.11, 2014/01/27 09:28:26
|
version 1.12, 2015/11/26 11:44:19
|
Line 1
|
Line 1
|
*> \brief \b DPSTF2 computes the Cholesky factorization with complete pivoting of a real symmetric or complex Hermitian positive semi-definite matrix. |
*> \brief \b DPSTF2 computes the Cholesky factorization with complete pivoting of a real symmetric positive semidefinite matrix. |
* |
* |
* =========== DOCUMENTATION =========== |
* =========== DOCUMENTATION =========== |
* |
* |
Line 121
|
Line 121
|
*> < 0: If INFO = -K, the K-th argument had an illegal value, |
*> < 0: If INFO = -K, the K-th argument had an illegal value, |
*> = 0: algorithm completed successfully, and |
*> = 0: algorithm completed successfully, and |
*> > 0: the matrix A is either rank deficient with computed rank |
*> > 0: the matrix A is either rank deficient with computed rank |
*> as returned in RANK, or is indefinite. See Section 7 of |
*> as returned in RANK, or is not positive semidefinite. See |
*> LAPACK Working Note #161 for further information. |
*> Section 7 of LAPACK Working Note #161 for further |
|
*> information. |
*> \endverbatim |
*> \endverbatim |
* |
* |
* Authors: |
* Authors: |
Line 133
|
Line 134
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date September 2012 |
*> \date November 2015 |
* |
* |
*> \ingroup doubleOTHERcomputational |
*> \ingroup doubleOTHERcomputational |
* |
* |
* ===================================================================== |
* ===================================================================== |
SUBROUTINE DPSTF2( UPLO, N, A, LDA, PIV, RANK, TOL, WORK, INFO ) |
SUBROUTINE DPSTF2( UPLO, N, A, LDA, PIV, RANK, TOL, WORK, INFO ) |
* |
* |
* -- LAPACK computational routine (version 3.4.2) -- |
* -- LAPACK computational routine (version 3.6.0) -- |
* -- 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..-- |
* September 2012 |
* November 2015 |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
DOUBLE PRECISION TOL |
DOUBLE PRECISION TOL |
Line 216
|
Line 217
|
AJJ = A( PVT, PVT ) |
AJJ = A( PVT, PVT ) |
END IF |
END IF |
END DO |
END DO |
IF( AJJ.EQ.ZERO.OR.DISNAN( AJJ ) ) THEN |
IF( AJJ.LE.ZERO.OR.DISNAN( AJJ ) ) THEN |
RANK = 0 |
RANK = 0 |
INFO = 1 |
INFO = 1 |
GO TO 170 |
GO TO 170 |