Diff for /rpl/lapack/lapack/zpotrf.f between versions 1.13 and 1.14

version 1.13, 2014/01/27 09:28:41 version 1.14, 2015/11/26 11:44:25
Line 100 Line 100
 *> \author Univ. of Colorado Denver   *> \author Univ. of Colorado Denver 
 *> \author NAG Ltd.   *> \author NAG Ltd. 
 *  *
 *> \date November 2011  *> \date November 2015
 *  *
 *> \ingroup complex16POcomputational  *> \ingroup complex16POcomputational
 *  *
 *  =====================================================================  *  =====================================================================
       SUBROUTINE ZPOTRF( UPLO, N, A, LDA, INFO )        SUBROUTINE ZPOTRF( UPLO, N, A, LDA, INFO )
 *  *
 *  -- LAPACK computational routine (version 3.4.0) --  *  -- 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..--
 *     November 2011  *     November 2015
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          UPLO        CHARACTER          UPLO
Line 137 Line 137
       EXTERNAL           LSAME, ILAENV        EXTERNAL           LSAME, ILAENV
 *     ..  *     ..
 *     .. External Subroutines ..  *     .. External Subroutines ..
       EXTERNAL           XERBLA, ZGEMM, ZHERK, ZPOTF2, ZTRSM        EXTERNAL           XERBLA, ZGEMM, ZHERK, ZPOTRF2, ZTRSM
 *     ..  *     ..
 *     .. Intrinsic Functions ..  *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN        INTRINSIC          MAX, MIN
Line 172 Line 172
 *  *
 *        Use unblocked code.  *        Use unblocked code.
 *  *
          CALL ZPOTF2( UPLO, N, A, LDA, INFO )           CALL ZPOTRF2( UPLO, N, A, LDA, INFO )
       ELSE        ELSE
 *  *
 *        Use blocked code.  *        Use blocked code.
Line 189 Line 189
                JB = MIN( NB, N-J+1 )                 JB = MIN( NB, N-J+1 )
                CALL ZHERK( 'Upper', 'Conjugate transpose', JB, J-1,                 CALL ZHERK( 'Upper', 'Conjugate transpose', JB, J-1,
      $                     -ONE, A( 1, J ), LDA, ONE, A( J, J ), LDA )       $                     -ONE, A( 1, J ), LDA, ONE, A( J, J ), LDA )
                CALL ZPOTF2( 'Upper', JB, A( J, J ), LDA, INFO )                 CALL ZPOTRF2( 'Upper', JB, A( J, J ), LDA, INFO )
                IF( INFO.NE.0 )                 IF( INFO.NE.0 )
      $            GO TO 30       $            GO TO 30
                IF( J+JB.LE.N ) THEN                 IF( J+JB.LE.N ) THEN
Line 218 Line 218
                JB = MIN( NB, N-J+1 )                 JB = MIN( NB, N-J+1 )
                CALL ZHERK( 'Lower', 'No transpose', JB, J-1, -ONE,                 CALL ZHERK( 'Lower', 'No transpose', JB, J-1, -ONE,
      $                     A( J, 1 ), LDA, ONE, A( J, J ), LDA )       $                     A( J, 1 ), LDA, ONE, A( J, J ), LDA )
                CALL ZPOTF2( 'Lower', JB, A( J, J ), LDA, INFO )                 CALL ZPOTRF2( 'Lower', JB, A( J, J ), LDA, INFO )
                IF( INFO.NE.0 )                 IF( INFO.NE.0 )
      $            GO TO 30       $            GO TO 30
                IF( J+JB.LE.N ) THEN                 IF( J+JB.LE.N ) THEN

Removed from v.1.13  
changed lines
  Added in v.1.14


CVSweb interface <joel.bertrand@systella.fr>