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

Diff for /rpl/lapack/lapack/zpotrf2.f between versions 1.2 and 1.6

version 1.2, 2016/08/27 15:35:04 version 1.6, 2020/05/21 21:46:10
Line 2 Line 2
 *  *
 *  =========== 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/
 *  *
 *  Definition:  *  Definition:
 *  ===========  *  ===========
 *  *
 *       RECURSIVE SUBROUTINE ZPOTRF2( UPLO, N, A, LDA, INFO )  *       RECURSIVE SUBROUTINE ZPOTRF2( UPLO, N, A, LDA, INFO )
 *   *
 *       .. Scalar Arguments ..  *       .. Scalar Arguments ..
 *       CHARACTER          UPLO  *       CHARACTER          UPLO
 *       INTEGER            INFO, LDA, N  *       INTEGER            INFO, LDA, N
Line 17 Line 17
 *       .. Array Arguments ..  *       .. Array Arguments ..
 *       COMPLEX*16         A( LDA, * )  *       COMPLEX*16         A( LDA, * )
 *       ..  *       ..
 *    *
 *  *
 *> \par Purpose:  *> \par Purpose:
 *  =============  *  =============
 *>  *>
 *> \verbatim  *> \verbatim
 *>  *>
 *> ZPOTRF2 computes the Cholesky factorization of a real symmetric  *> ZPOTRF2 computes the Cholesky factorization of a Hermitian
 *> positive definite matrix A using the recursive algorithm.  *> positive definite matrix A using the recursive algorithm.
 *>  *>
 *> The factorization has the form  *> The factorization has the form
Line 41 Line 41
 *>  *>
 *> The subroutine calls itself to factor A11. Update and scale A21  *> The subroutine calls itself to factor A11. Update and scale A21
 *> or A12, update A22 then call itself to factor A22.  *> or A12, update A22 then call itself to factor A22.
 *>   *>
 *> \endverbatim  *> \endverbatim
 *  *
 *  Arguments:  *  Arguments:
Line 63 Line 63
 *> \param[in,out] A  *> \param[in,out] A
 *> \verbatim  *> \verbatim
 *>          A is COMPLEX*16 array, dimension (LDA,N)  *>          A is COMPLEX*16 array, dimension (LDA,N)
 *>          On entry, the symmetric matrix A.  If UPLO = 'U', the leading  *>          On entry, the Hermitian matrix A.  If UPLO = 'U', the leading
 *>          N-by-N upper triangular part of A contains the upper  *>          N-by-N upper triangular part of A contains the upper
 *>          triangular part of the matrix A, and the strictly lower  *>          triangular part of the matrix A, and the strictly lower
 *>          triangular part of A is not referenced.  If UPLO = 'L', the  *>          triangular part of A is not referenced.  If UPLO = 'L', the
Line 94 Line 94
 *  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 November 2015  *> \date December 2016
 *  *
 *> \ingroup complex16POcomputational  *> \ingroup complex16POcomputational
 *  *
 *  =====================================================================  *  =====================================================================
       RECURSIVE SUBROUTINE ZPOTRF2( UPLO, N, A, LDA, INFO )        RECURSIVE SUBROUTINE ZPOTRF2( UPLO, N, A, LDA, INFO )
 *  *
 *  -- LAPACK computational routine (version 3.6.0) --  *  -- 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..--
 *     November 2015  *     December 2016
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          UPLO        CHARACTER          UPLO
Line 128 Line 128
       PARAMETER          ( CONE = (1.0D+0, 0.0D+0) )        PARAMETER          ( CONE = (1.0D+0, 0.0D+0) )
 *     ..  *     ..
 *     .. Local Scalars ..  *     .. Local Scalars ..
       LOGICAL            UPPER                    LOGICAL            UPPER
       INTEGER            N1, N2, IINFO        INTEGER            N1, N2, IINFO
       DOUBLE PRECISION   AJJ        DOUBLE PRECISION   AJJ
 *     ..  *     ..
Line 193 Line 193
          IF ( IINFO.NE.0 ) THEN           IF ( IINFO.NE.0 ) THEN
             INFO = IINFO              INFO = IINFO
             RETURN              RETURN
          END IF               END IF
 *  *
 *        Compute the Cholesky factorization A = U**H*U  *        Compute the Cholesky factorization A = U**H*U
 *  *
Line 205 Line 205
      $                  A( 1, 1 ), LDA, A( 1, N1+1 ), LDA )       $                  A( 1, 1 ), LDA, A( 1, N1+1 ), LDA )
 *  *
 *           Update and factor A22  *           Update and factor A22
 *            *
             CALL ZHERK( UPLO, 'C', N2, N1, -ONE, A( 1, N1+1 ), LDA,              CALL ZHERK( UPLO, 'C', N2, N1, -ONE, A( 1, N1+1 ), LDA,
      $                  ONE, A( N1+1, N1+1 ), LDA )       $                  ONE, A( N1+1, N1+1 ), LDA )
             CALL ZPOTRF2( UPLO, N2, A( N1+1, N1+1 ), LDA, IINFO )              CALL ZPOTRF2( UPLO, N2, A( N1+1, N1+1 ), LDA, IINFO )
Removed from v.1.2  
changed lines
  Added in v.1.6

CVSweb interface <joel.bertrand@systella.fr>