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

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

version 1.12, 2012/12/14 14:22:28 version 1.13, 2014/01/27 09:24:34
Line 121 Line 121
 *> \author Univ. of Colorado Denver   *> \author Univ. of Colorado Denver 
 *> \author NAG Ltd.   *> \author NAG Ltd. 
 *  *
 *> \date November 2011  *> \date November 2013
 *  *
 *> \ingroup doubleGEcomputational  *> \ingroup doubleGEcomputational
 *  *
Line 160 Line 160
 *  =====================================================================  *  =====================================================================
       SUBROUTINE DGEBAL( JOB, N, A, LDA, ILO, IHI, SCALE, INFO )        SUBROUTINE DGEBAL( JOB, N, A, LDA, ILO, IHI, SCALE, INFO )
 *  *
 *  -- LAPACK computational routine (version 3.4.0) --  *  -- LAPACK computational routine (version 3.5.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 2013
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          JOB        CHARACTER          JOB
Line 192 Line 192
 *     .. External Functions ..  *     .. External Functions ..
       LOGICAL            DISNAN, LSAME        LOGICAL            DISNAN, LSAME
       INTEGER            IDAMAX        INTEGER            IDAMAX
       DOUBLE PRECISION   DLAMCH        DOUBLE PRECISION   DLAMCH, DNRM2
       EXTERNAL           DISNAN, LSAME, IDAMAX, DLAMCH        EXTERNAL           DISNAN, LSAME, IDAMAX, DLAMCH, DNRM2
 *     ..  *     ..
 *     .. External Subroutines ..  *     .. External Subroutines ..
       EXTERNAL           DSCAL, DSWAP, XERBLA        EXTERNAL           DSCAL, DSWAP, XERBLA
Line 312 Line 312
       SFMAX1 = ONE / SFMIN1        SFMAX1 = ONE / SFMIN1
       SFMIN2 = SFMIN1*SCLFAC        SFMIN2 = SFMIN1*SCLFAC
       SFMAX2 = ONE / SFMIN2        SFMAX2 = ONE / SFMIN2
   *
   140 CONTINUE    140 CONTINUE
       NOCONV = .FALSE.        NOCONV = .FALSE.
 *  *
       DO 200 I = K, L        DO 200 I = K, L
          C = ZERO  
          R = ZERO  
 *  *
          DO 150 J = K, L           C = DNRM2( L-K+1, A( K, I ), 1 )
             IF( J.EQ.I )           R = DNRM2( L-K+1, A( I, K ), LDA )
      $         GO TO 150  
             C = C + ABS( A( J, I ) )  
             R = R + ABS( A( I, J ) )  
   150    CONTINUE  
          ICA = IDAMAX( L, A( 1, I ), 1 )           ICA = IDAMAX( L, A( 1, I ), 1 )
          CA = ABS( A( ICA, I ) )           CA = ABS( A( ICA, I ) )
          IRA = IDAMAX( N-K+1, A( I, K ), LDA )           IRA = IDAMAX( N-K+1, A( I, K ), LDA )
Removed from v.1.12  
changed lines
  Added in v.1.13

CVSweb interface <joel.bertrand@systella.fr>