Diff for /rpl/lapack/lapack/dlaein.f between versions 1.8 and 1.9

version 1.8, 2010/12/21 13:53:29 version 1.9, 2011/07/22 07:38:06
Line 1 Line 1
       SUBROUTINE DLAEIN( RIGHTV, NOINIT, N, H, LDH, WR, WI, VR, VI, B,        SUBROUTINE DLAEIN( RIGHTV, NOINIT, N, H, LDH, WR, WI, VR, VI, B,
      $                   LDB, WORK, EPS3, SMLNUM, BIGNUM, INFO )       $                   LDB, WORK, EPS3, SMLNUM, BIGNUM, INFO )
 *  *
 *  -- LAPACK auxiliary routine (version 3.2.2) --  *  -- LAPACK auxiliary 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..--
 *     June 2010  *  -- April 2011                                                      --
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       LOGICAL            NOINIT, RIGHTV        LOGICAL            NOINIT, RIGHTV
Line 232 Line 232
          DO 110 ITS = 1, N           DO 110 ITS = 1, N
 *  *
 *           Solve U*x = scale*v for a right eigenvector  *           Solve U*x = scale*v for a right eigenvector
 *             or U'*x = scale*v for a left eigenvector,  *             or U**T*x = scale*v for a left eigenvector,
 *           overwriting x on v.  *           overwriting x on v.
 *  *
             CALL DLATRS( 'Upper', TRANS, 'Nonunit', NORMIN, N, B, LDB,              CALL DLATRS( 'Upper', TRANS, 'Nonunit', NORMIN, N, B, LDB,
Line 427 Line 427
             VCRIT = BIGNUM              VCRIT = BIGNUM
 *  *
 *           Solve U*(xr,xi) = scale*(vr,vi) for a right eigenvector,  *           Solve U*(xr,xi) = scale*(vr,vi) for a right eigenvector,
 *             or U'*(xr,xi) = scale*(vr,vi) for a left eigenvector,  *             or U**T*(xr,xi) = scale*(vr,vi) for a left eigenvector,
 *           overwriting (xr,xi) on (vr,vi).  *           overwriting (xr,xi) on (vr,vi).
 *  *
             DO 250 I = I1, I2, I3              DO 250 I = I1, I2, I3

Removed from v.1.8  
changed lines
  Added in v.1.9


CVSweb interface <joel.bertrand@systella.fr>