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

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

version 1.13, 2014/01/27 09:28:26 version 1.14, 2015/11/26 11:44:19
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 doublePOcomputational  *> \ingroup doublePOcomputational
 *  *
 *  =====================================================================  *  =====================================================================
       SUBROUTINE DPOTRF( UPLO, N, A, LDA, INFO )        SUBROUTINE DPOTRF( 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 136 Line 136
       EXTERNAL           LSAME, ILAENV        EXTERNAL           LSAME, ILAENV
 *     ..  *     ..
 *     .. External Subroutines ..  *     .. External Subroutines ..
       EXTERNAL           DGEMM, DPOTF2, DSYRK, DTRSM, XERBLA        EXTERNAL           DGEMM, DPOTRF2, DSYRK, DTRSM, XERBLA
 *     ..  *     ..
 *     .. Intrinsic Functions ..  *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN        INTRINSIC          MAX, MIN
Line 171 Line 171
 *  *
 *        Use unblocked code.  *        Use unblocked code.
 *  *
          CALL DPOTF2( UPLO, N, A, LDA, INFO )           CALL DPOTRF2( UPLO, N, A, LDA, INFO )
       ELSE        ELSE
 *  *
 *        Use blocked code.  *        Use blocked code.
Line 188 Line 188
                JB = MIN( NB, N-J+1 )                 JB = MIN( NB, N-J+1 )
                CALL DSYRK( 'Upper', 'Transpose', JB, J-1, -ONE,                 CALL DSYRK( 'Upper', 'Transpose', JB, J-1, -ONE,
      $                     A( 1, J ), LDA, ONE, A( J, J ), LDA )       $                     A( 1, J ), LDA, ONE, A( J, J ), LDA )
                CALL DPOTF2( 'Upper', JB, A( J, J ), LDA, INFO )                 CALL DPOTRF2( '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 216 Line 216
                JB = MIN( NB, N-J+1 )                 JB = MIN( NB, N-J+1 )
                CALL DSYRK( 'Lower', 'No transpose', JB, J-1, -ONE,                 CALL DSYRK( '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 DPOTF2( 'Lower', JB, A( J, J ), LDA, INFO )                 CALL DPOTRF2( '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>