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

version 1.12, 2012/12/14 14:22:44 version 1.13, 2014/01/27 09:24:36
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 complex16GEcomputational  *> \ingroup complex16GEcomputational
 *  *
Line 160 Line 160
 *  =====================================================================  *  =====================================================================
       SUBROUTINE ZGEBAL( JOB, N, A, LDA, ILO, IHI, SCALE, INFO )        SUBROUTINE ZGEBAL( 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 194 Line 194
 *     .. External Functions ..  *     .. External Functions ..
       LOGICAL            DISNAN, LSAME        LOGICAL            DISNAN, LSAME
       INTEGER            IZAMAX        INTEGER            IZAMAX
       DOUBLE PRECISION   DLAMCH        DOUBLE PRECISION   DLAMCH, DZNRM2
       EXTERNAL           DISNAN, LSAME, IZAMAX, DLAMCH        EXTERNAL           DISNAN, LSAME, IZAMAX, DLAMCH, DZNRM2
 *     ..  *     ..
 *     .. External Subroutines ..  *     .. External Subroutines ..
       EXTERNAL           XERBLA, ZDSCAL, ZSWAP        EXTERNAL           XERBLA, ZDSCAL, ZSWAP
Line 324 Line 324
       NOCONV = .FALSE.        NOCONV = .FALSE.
 *  *
       DO 200 I = K, L        DO 200 I = K, L
          C = ZERO  
          R = ZERO  
 *  *
          DO 150 J = K, L           C = DZNRM2( L-K+1, A( K, I ), 1 )
             IF( J.EQ.I )           R = DZNRM2( L-K+1, A( I, K ), LDA )
      $         GO TO 150  
             C = C + CABS1( A( J, I ) )  
             R = R + CABS1( A( I, J ) )  
   150    CONTINUE  
          ICA = IZAMAX( L, A( 1, I ), 1 )           ICA = IZAMAX( L, A( 1, I ), 1 )
          CA = ABS( A( ICA, I ) )           CA = ABS( A( ICA, I ) )
          IRA = IZAMAX( N-K+1, A( I, K ), LDA )           IRA = IZAMAX( 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>