diff options
| author | julie <julielangou@users.noreply.github.com> | 2011-04-07 13:53:54 +0000 |
|---|---|---|
| committer | julie <julielangou@users.noreply.github.com> | 2011-04-07 13:53:54 +0000 |
| commit | 16973f0c15a4327f8a78b9652f178246895778c1 (patch) | |
| tree | 1d022faeffac8e1524f4f73517483f4cd828038d /SRC/dtgsy2.f | |
| parent | 8b5dd5433b205f27821b10d5389f55bcd4a1e87a (diff) | |
| download | lapack-16973f0c15a4327f8a78b9652f178246895778c1.tar.gz lapack-16973f0c15a4327f8a78b9652f178246895778c1.tar.bz2 lapack-16973f0c15a4327f8a78b9652f178246895778c1.zip | |
Second and hopefully last pass to homgenize notation for transpose (**T) and conjugate transpose (**H)
Corresponds to bug0024
Please take a look and let me know if you find some old notation of transpose.
I am going to close bug0024.
Julie
Diffstat (limited to 'SRC/dtgsy2.f')
| -rw-r--r-- | SRC/dtgsy2.f | 28 |
1 files changed, 14 insertions, 14 deletions
diff --git a/SRC/dtgsy2.f b/SRC/dtgsy2.f index 72dcf698..bfaa9ae5 100644 --- a/SRC/dtgsy2.f +++ b/SRC/dtgsy2.f @@ -38,19 +38,19 @@ * In matrix notation solving equation (1) corresponds to solve * Z*x = 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**T, Im) ] (2) +* [ kron(In, D) -kron(E**T, Im) ], * -* Ik is the identity matrix of size k and X' is the transpose of X. +* Ik is the identity matrix of size k and X**T is the transpose of X. * kron(X, Y) is the Kronecker product between the matrices X and Y. * In the process of solving (1), we solve a number of such systems * where Dim(In), Dim(In) = 1 or 2. * -* If TRANS = 'T', solve the transposed system Z'*y = scale*b for y, +* If TRANS = 'T', solve the transposed system Z**T*y = scale*b for y, * 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 is used to compute an estimate of Dif[(A, D), (B, E)] = * sigma_min(Z) using reverse communicaton with DLACON. @@ -649,7 +649,7 @@ ZDIM = MB*NB*2 IF( ( MB.EQ.1 ) .AND. ( NB.EQ.1 ) ) THEN * -* Build a 2-by-2 system Z' * x = RHS +* Build a 2-by-2 system Z**T * x = RHS * Z( 1, 1 ) = A( IS, IS ) Z( 2, 1 ) = -B( JS, JS ) @@ -661,7 +661,7 @@ RHS( 1 ) = C( IS, JS ) RHS( 2 ) = F( IS, JS ) * -* Solve Z' * x = RHS +* Solve Z**T * x = RHS * CALL DGETC2( ZDIM, Z, LDZ, IPIV, JPIV, IERR ) IF( IERR.GT.0 ) @@ -703,7 +703,7 @@ * ELSE IF( ( MB.EQ.1 ) .AND. ( NB.EQ.2 ) ) THEN * -* Build a 4-by-4 system Z' * x = RHS +* Build a 4-by-4 system Z**T * x = RHS * Z( 1, 1 ) = A( IS, IS ) Z( 2, 1 ) = ZERO @@ -732,7 +732,7 @@ RHS( 3 ) = F( IS, JS ) RHS( 4 ) = F( IS, JSP1 ) * -* Solve Z' * x = RHS +* Solve Z**T * x = RHS * CALL DGETC2( ZDIM, Z, LDZ, IPIV, JPIV, IERR ) IF( IERR.GT.0 ) @@ -775,7 +775,7 @@ * ELSE IF( ( MB.EQ.2 ) .AND. ( NB.EQ.1 ) ) THEN * -* Build a 4-by-4 system Z' * x = RHS +* Build a 4-by-4 system Z**T * x = RHS * Z( 1, 1 ) = A( IS, IS ) Z( 2, 1 ) = A( IS, ISP1 ) @@ -804,7 +804,7 @@ RHS( 3 ) = F( IS, JS ) RHS( 4 ) = F( ISP1, JS ) * -* Solve Z' * x = RHS +* Solve Z**T * x = RHS * CALL DGETC2( ZDIM, Z, LDZ, IPIV, JPIV, IERR ) IF( IERR.GT.0 ) @@ -846,7 +846,7 @@ * ELSE IF( ( MB.EQ.2 ) .AND. ( NB.EQ.2 ) ) THEN * -* Build an 8-by-8 system Z' * x = RHS +* Build an 8-by-8 system Z**T * x = RHS * CALL DLASET( 'F', LDZ, LDZ, ZERO, ZERO, Z, LDZ ) * @@ -898,7 +898,7 @@ 160 CONTINUE * * -* Solve Z' * x = RHS +* Solve Z**T * x = RHS * CALL DGETC2( ZDIM, Z, LDZ, IPIV, JPIV, IERR ) IF( IERR.GT.0 ) |