diff options
Diffstat (limited to 'SRC/dtrsyl.f')
| -rw-r--r-- | SRC/dtrsyl.f | 26 |
1 files changed, 13 insertions, 13 deletions
diff --git a/SRC/dtrsyl.f b/SRC/dtrsyl.f index ec3e23d7..e63db3c5 100644 --- a/SRC/dtrsyl.f +++ b/SRC/dtrsyl.f @@ -355,17 +355,17 @@ * ELSE IF( .NOT.NOTRNA .AND. NOTRNB ) THEN * -* Solve A' *X + ISGN*X*B = scale*C. +* Solve A**T *X + ISGN*X*B = scale*C. * * The (K,L)th block of X is determined starting from * upper-left corner column by column by * -* A(K,K)'*X(K,L) + ISGN*X(K,L)*B(L,L) = C(K,L) - R(K,L) +* A(K,K)**T*X(K,L) + ISGN*X(K,L)*B(L,L) = C(K,L) - R(K,L) * * Where -* K-1 L-1 -* R(K,L) = SUM [A(I,K)'*X(I,L)] +ISGN*SUM [X(K,J)*B(J,L)] -* I=1 J=1 +* K-1 L-1 +* R(K,L) = SUM [A(I,K)**T*X(I,L)] +ISGN*SUM [X(K,J)*B(J,L)] +* I=1 J=1 * * Start column loop (index = L) * L1 (L2): column index of the first (last) row of X(K,L) @@ -530,17 +530,17 @@ * ELSE IF( .NOT.NOTRNA .AND. .NOT.NOTRNB ) THEN * -* Solve A'*X + ISGN*X*B' = scale*C. +* Solve A**T*X + ISGN*X*B**T = scale*C. * * The (K,L)th block of X is determined starting from * top-right corner column by column by * -* A(K,K)'*X(K,L) + ISGN*X(K,L)*B(L,L)' = C(K,L) - R(K,L) +* A(K,K)**T*X(K,L) + ISGN*X(K,L)*B(L,L)**T = C(K,L) - R(K,L) * * Where -* K-1 N -* R(K,L) = SUM [A(I,K)'*X(I,L)] + ISGN*SUM [X(K,J)*B(L,J)']. -* I=1 J=L+1 +* K-1 N +* R(K,L) = SUM [A(I,K)**T*X(I,L)] + ISGN*SUM [X(K,J)*B(L,J)**T]. +* I=1 J=L+1 * * Start column loop (index = L) * L1 (L2): column index of the first (last) row of X(K,L) @@ -714,16 +714,16 @@ * ELSE IF( NOTRNA .AND. .NOT.NOTRNB ) THEN * -* Solve A*X + ISGN*X*B' = scale*C. +* Solve A*X + ISGN*X*B**T = scale*C. * * The (K,L)th block of X is determined starting from * bottom-right corner column by column by * -* A(K,K)*X(K,L) + ISGN*X(K,L)*B(L,L)' = C(K,L) - R(K,L) +* A(K,K)*X(K,L) + ISGN*X(K,L)*B(L,L)**T = C(K,L) - R(K,L) * * Where * M N -* R(K,L) = SUM [A(K,I)*X(I,L)] + ISGN*SUM [X(K,J)*B(L,J)']. +* R(K,L) = SUM [A(K,I)*X(I,L)] + ISGN*SUM [X(K,J)*B(L,J)**T]. * I=K+1 J=L+1 * * Start column loop (index = L) |