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

Diff for /rpl/lapack/lapack/dsterf.f between versions 1.7 and 1.8

version 1.7, 2010/12/21 13:53:38 version 1.8, 2011/07/22 07:38:11
Line 1 Line 1
       SUBROUTINE DSTERF( N, D, E, INFO )        SUBROUTINE DSTERF( N, D, E, INFO )
 *  *
 *  -- LAPACK routine (version 3.2) --  *  -- LAPACK routine (version 3.3.1) --
 *  -- 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 2006  *  -- April 2011                                                      --
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       INTEGER            INFO, N        INTEGER            INFO, N
Line 54 Line 54
      $                   NMAXIT       $                   NMAXIT
       DOUBLE PRECISION   ALPHA, ANORM, BB, C, EPS, EPS2, GAMMA, OLDC,        DOUBLE PRECISION   ALPHA, ANORM, BB, C, EPS, EPS2, GAMMA, OLDC,
      $                   OLDGAM, P, R, RT1, RT2, RTE, S, SAFMAX, SAFMIN,       $                   OLDGAM, P, R, RT1, RT2, RTE, S, SAFMAX, SAFMIN,
      $                   SIGMA, SSFMAX, SSFMIN       $                   SIGMA, SSFMAX, SSFMIN, RMAX
 *     ..  *     ..
 *     .. External Functions ..  *     .. External Functions ..
       DOUBLE PRECISION   DLAMCH, DLANST, DLAPY2        DOUBLE PRECISION   DLAMCH, DLANST, DLAPY2
Line 90 Line 90
       SAFMAX = ONE / SAFMIN        SAFMAX = ONE / SAFMIN
       SSFMAX = SQRT( SAFMAX ) / THREE        SSFMAX = SQRT( SAFMAX ) / THREE
       SSFMIN = SQRT( SAFMIN ) / EPS2        SSFMIN = SQRT( SAFMIN ) / EPS2
         RMAX = DLAMCH( 'O' )
 *  *
 *     Compute the eigenvalues of the tridiagonal matrix.  *     Compute the eigenvalues of the tridiagonal matrix.
 *  *
Line 128 Line 129
 *  *
 *     Scale submatrix in rows and columns L to LEND  *     Scale submatrix in rows and columns L to LEND
 *  *
       ANORM = DLANST( 'I', LEND-L+1, D( L ), E( L ) )        ANORM = DLANST( 'M', LEND-L+1, D( L ), E( L ) )
       ISCALE = 0        ISCALE = 0
       IF( ANORM.GT.SSFMAX ) THEN        IF( ANORM.EQ.ZERO )
        $   GO TO 10      
         IF( (ANORM.GT.SSFMAX) ) THEN
          ISCALE = 1           ISCALE = 1
          CALL DLASCL( 'G', 0, 0, ANORM, SSFMAX, LEND-L+1, 1, D( L ), N,           CALL DLASCL( 'G', 0, 0, ANORM, SSFMAX, LEND-L+1, 1, D( L ), N,
      $                INFO )       $                INFO )
Removed from v.1.7  
changed lines
  Added in v.1.8

CVSweb interface <joel.bertrand@systella.fr>