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

version 1.2, 2016/08/27 15:35:04 version 1.3, 2017/06/17 10:54:25
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:
 *  =============  *  =============
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 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.3


CVSweb interface <joel.bertrand@systella.fr>