[BACK]Return to dpstf2.f CVS log [TXT] [DIR] Up to [local] / rpl / lapack / lapack

Diff for /rpl/lapack/lapack/dpstf2.f between versions 1.9 and 1.15

version 1.9, 2012/12/14 12:30:27 version 1.15, 2017/06/17 11:06:31
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 ===========
 *  *
 * Online html documentation available at   * Online html documentation available at
 *            http://www.netlib.org/lapack/explore-html/   *            http://www.netlib.org/lapack/explore-html/
 *  *
 *> \htmlonly  *> \htmlonly
 *> Download DPSTF2 + dependencies   *> Download DPSTF2 + dependencies
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dpstf2.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dpstf2.f">
 *> [TGZ]</a>   *> [TGZ]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dpstf2.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dpstf2.f">
 *> [ZIP]</a>   *> [ZIP]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dpstf2.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dpstf2.f">
 *> [TXT]</a>  *> [TXT]</a>
 *> \endhtmlonly   *> \endhtmlonly
 *  *
 *  Definition:  *  Definition:
 *  ===========  *  ===========
 *  *
 *       SUBROUTINE DPSTF2( UPLO, N, A, LDA, PIV, RANK, TOL, WORK, INFO )  *       SUBROUTINE DPSTF2( UPLO, N, A, LDA, PIV, RANK, TOL, WORK, INFO )
 *   *
 *       .. Scalar Arguments ..  *       .. Scalar Arguments ..
 *       DOUBLE PRECISION   TOL  *       DOUBLE PRECISION   TOL
 *       INTEGER            INFO, LDA, N, RANK  *       INTEGER            INFO, LDA, N, RANK
Line 29 Line 29
 *       DOUBLE PRECISION   A( LDA, * ), WORK( 2*N )  *       DOUBLE PRECISION   A( LDA, * ), WORK( 2*N )
 *       INTEGER            PIV( N )  *       INTEGER            PIV( N )
 *       ..  *       ..
 *    *
 *  *
 *> \par Purpose:  *> \par Purpose:
 *  =============  *  =============
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:
 *  ========  *  ========
 *  *
 *> \author Univ. of Tennessee   *> \author Univ. of Tennessee
 *> \author Univ. of California Berkeley   *> \author Univ. of California Berkeley
 *> \author Univ. of Colorado Denver   *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.   *> \author NAG Ltd.
 *  *
 *> \date September 2012  *> \date December 2016
 *  *
 *> \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.7.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  *     December 2016
 *  *
 *     .. 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
Removed from v.1.9  
changed lines
  Added in v.1.15

CVSweb interface <joel.bertrand@systella.fr>