--- rpl/lapack/lapack/dlaein.f	2010/04/21 13:45:16	1.2
+++ rpl/lapack/lapack/dlaein.f	2011/07/22 07:38:06	1.9
@@ -1,10 +1,10 @@
       SUBROUTINE DLAEIN( RIGHTV, NOINIT, N, H, LDH, WR, WI, VR, VI, B,
      $                   LDB, WORK, EPS3, SMLNUM, BIGNUM, INFO )
 *
-*  -- LAPACK auxiliary routine (version 3.2) --
+*  -- LAPACK auxiliary routine (version 3.3.1) --
 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-*     November 2006
+*  -- April 2011                                                      --
 *
 *     .. Scalar Arguments ..
       LOGICAL            NOINIT, RIGHTV
@@ -26,11 +26,11 @@
 *  Arguments
 *  =========
 *
-*  RIGHTV   (input) LOGICAL
+*  RIGHTV  (input) LOGICAL
 *          = .TRUE. : compute right eigenvector;
 *          = .FALSE.: compute left eigenvector.
 *
-*  NOINIT   (input) LOGICAL
+*  NOINIT  (input) LOGICAL
 *          = .TRUE. : no initial vector supplied in (VR,VI).
 *          = .FALSE.: initial vector supplied in (VR,VI).
 *
@@ -70,7 +70,7 @@
 *  LDB     (input) INTEGER
 *          The leading dimension of the array B.  LDB >= N+1.
 *
-*  WORK   (workspace) DOUBLE PRECISION array, dimension (N)
+*  WORK    (workspace) DOUBLE PRECISION array, dimension (N)
 *
 *  EPS3    (input) DOUBLE PRECISION
 *          A small machine-dependent value which is used to perturb
@@ -232,7 +232,7 @@
          DO 110 ITS = 1, N
 *
 *           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.
 *
             CALL DLATRS( 'Upper', TRANS, 'Nonunit', NORMIN, N, B, LDB,
@@ -427,7 +427,7 @@
             VCRIT = BIGNUM
 *
 *           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).
 *
             DO 250 I = I1, I2, I3