--- rpl/lapack/lapack/ztgsyl.f	2010/12/21 13:53:57	1.7
+++ rpl/lapack/lapack/ztgsyl.f	2011/07/22 07:38:21	1.8
@@ -2,10 +2,10 @@
      $                   LDD, E, LDE, F, LDF, SCALE, DIF, WORK, LWORK,
      $                   IWORK, INFO )
 *
-*  -- LAPACK routine (version 3.2) --
+*  -- LAPACK 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..--
-*     January 2007
+*  -- April 2011                                                      --
 *
 *     .. Scalar Arguments ..
       CHARACTER          TRANS
@@ -39,18 +39,18 @@
 *  In matrix notation (1) is equivalent to solve Zx = scale*b, where Z
 *  is defined as
 *
-*         Z = [ kron(In, A)  -kron(B', Im) ]        (2)
-*             [ kron(In, D)  -kron(E', Im) ],
+*         Z = [ kron(In, A)  -kron(B**H, Im) ]        (2)
+*             [ kron(In, D)  -kron(E**H, Im) ],
 *
-*  Here Ix is the identity matrix of size x and X' is the conjugate
+*  Here Ix is the identity matrix of size x and X**H is the conjugate
 *  transpose of X. Kron(X, Y) is the Kronecker product between the
 *  matrices X and Y.
 *
-*  If TRANS = 'C', y in the conjugate transposed system Z'*y = scale*b
+*  If TRANS = 'C', y in the conjugate transposed system Z**H *y = scale*b
 *  is solved for, which is equivalent to solve for R and L in
 *
-*              A' * R + D' * L = scale * C           (3)
-*              R * B' + L * E' = scale * -F
+*              A**H * R + D**H * L = scale * C           (3)
+*              R * B**H + L * E**H = scale * -F
 *
 *  This case (TRANS = 'C') is used to compute an one-norm-based estimate
 *  of Dif[(A,D), (B,E)], the separation between the matrix pairs (A,D)
@@ -492,7 +492,7 @@
       ELSE
 *
 *        Solve transposed (I, J)-subsystem
-*            A(I, I)' * R(I, J) + D(I, I)' * L(I, J) = C(I, J)
+*            A(I, I)**H * R(I, J) + D(I, I)**H * L(I, J) = C(I, J)
 *            R(I, J) * B(J, J)  + L(I, J) * E(J, J) = -F(I, J)
 *        for I = 1,2,..., P; J = Q, Q-1,..., 1
 *