diff options
Diffstat (limited to 'SRC/stgsyl.f')
| -rw-r--r-- | SRC/stgsyl.f | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/SRC/stgsyl.f b/SRC/stgsyl.f index beff0470..f2c2a940 100644 --- a/SRC/stgsyl.f +++ b/SRC/stgsyl.f @@ -46,11 +46,11 @@ * Here Ik is the identity matrix of size k and X' is the transpose of * X. kron(X, Y) is the Kronecker product between the matrices X and Y. * -* If TRANS = 'T', STGSYL solves the transposed system Z'*y = scale*b, +* If TRANS = 'T', STGSYL solves the transposed system Z**T*y = scale*b, * which is equivalent to solve for R and L in * -* A' * R + D' * L = scale * C (3) -* R * B' + L * E' = scale * (-F) +* A**T * R + D**T * L = scale * C (3) +* R * B**T + L * E**T = scale * -F * * This case (TRANS = 'T') is used to compute an one-norm-based estimate * of Dif[(A,D), (B,E)], the separation between the matrix pairs (A,D) @@ -485,8 +485,8 @@ ELSE * * Solve transposed (I, J)-subsystem -* A(I, I)' * R(I, J) + D(I, I)' * L(I, J) = C(I, J) -* R(I, J) * B(J, J)' + L(I, J) * E(J, J)' = -F(I, J) +* A(I, I)**T * R(I, J) + D(I, I)**T * L(I, J) = C(I, J) +* R(I, J) * B(J, J)**T + L(I, J) * E(J, J)**T = -F(I, J) * for I = 1,2,..., P; J = Q, Q-1,..., 1 * SCALE = ONE |