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

Diff for /rpl/lapack/lapack/dstein.f between versions 1.12 and 1.13

version 1.12, 2014/01/27 09:28:27 version 1.13, 2015/11/26 11:44:20
Line 166 Line 166
 *> \author Univ. of Colorado Denver   *> \author Univ. of Colorado Denver 
 *> \author NAG Ltd.   *> \author NAG Ltd. 
 *  *
 *> \date November 2011  *> \date November 2015
 *  *
 *> \ingroup doubleOTHERcomputational  *> \ingroup doubleOTHERcomputational
 *  *
Line 174 Line 174
       SUBROUTINE DSTEIN( N, D, E, M, W, IBLOCK, ISPLIT, Z, LDZ, WORK,        SUBROUTINE DSTEIN( N, D, E, M, W, IBLOCK, ISPLIT, Z, LDZ, WORK,
      $                   IWORK, IFAIL, INFO )       $                   IWORK, IFAIL, 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 ..
       INTEGER            INFO, LDZ, M, N        INTEGER            INFO, LDZ, M, N
Line 297 Line 297
          BLKSIZ = BN - B1 + 1           BLKSIZ = BN - B1 + 1
          IF( BLKSIZ.EQ.1 )           IF( BLKSIZ.EQ.1 )
      $      GO TO 60       $      GO TO 60
          GPIND = B1           GPIND = J1
 *  *
 *        Compute reorthogonalization criterion and stopping criterion.  *        Compute reorthogonalization criterion and stopping criterion.
 *  *
Line 370 Line 370
 *  *
 *           Normalize and scale the righthand side vector Pb.  *           Normalize and scale the righthand side vector Pb.
 *  *
               JMAX = IDAMAX( BLKSIZ, WORK( INDRV1+1 ), 1 )
             SCL = BLKSIZ*ONENRM*MAX( EPS,              SCL = BLKSIZ*ONENRM*MAX( EPS,
      $            ABS( WORK( INDRV4+BLKSIZ ) ) ) /       $            ABS( WORK( INDRV4+BLKSIZ ) ) ) /
      $            DASUM( BLKSIZ, WORK( INDRV1+1 ), 1 )       $            ABS( WORK( INDRV1+JMAX ) )
             CALL DSCAL( BLKSIZ, SCL, WORK( INDRV1+1 ), 1 )              CALL DSCAL( BLKSIZ, SCL, WORK( INDRV1+1 ), 1 )
 *  *
 *           Solve the system LU = Pb.  *           Solve the system LU = Pb.
Removed from v.1.12  
changed lines
  Added in v.1.13

CVSweb interface <joel.bertrand@systella.fr>