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

Diff for /rpl/lapack/lapack/ztgsyl.f between versions 1.7 and 1.8

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

CVSweb interface <joel.bertrand@systella.fr>