diff options
| author | julie <julielangou@users.noreply.github.com> | 2011-04-02 11:08:56 +0000 |
|---|---|---|
| committer | julie <julielangou@users.noreply.github.com> | 2011-04-02 11:08:56 +0000 |
| commit | f2953573ede24d7f8c01fdb18de48f65f00a9943 (patch) | |
| tree | 53172aa9083b9aa1abe2d6c130f7c173d8d8725b /SRC/ssytrs.f | |
| parent | 53b71f5605f83d116ab6bcf477bfb6d2ca757de1 (diff) | |
| download | lapack-f2953573ede24d7f8c01fdb18de48f65f00a9943.tar.gz lapack-f2953573ede24d7f8c01fdb18de48f65f00a9943.tar.bz2 lapack-f2953573ede24d7f8c01fdb18de48f65f00a9943.zip | |
First pass to homgenize notation for transpose (**T) and conjugate transpose (**H)
Corresponds to bug0024
Diffstat (limited to 'SRC/ssytrs.f')
| -rw-r--r-- | SRC/ssytrs.f | 16 |
1 files changed, 8 insertions, 8 deletions
diff --git a/SRC/ssytrs.f b/SRC/ssytrs.f index 9d2c8066..80600ea2 100644 --- a/SRC/ssytrs.f +++ b/SRC/ssytrs.f @@ -107,7 +107,7 @@ * IF( UPPER ) THEN * -* Solve A*X = B, where A = U*D*U'. +* Solve A*X = B, where A = U*D*U**T. * * First solve U*D*X = B, overwriting B with X. * @@ -178,7 +178,7 @@ GO TO 10 30 CONTINUE * -* Next solve U'*X = B, overwriting B with X. +* Next solve U**T *X = B, overwriting B with X. * * K is the main loop index, increasing from 1 to N in steps of * 1 or 2, depending on the size of the diagonal blocks. @@ -195,7 +195,7 @@ * * 1 x 1 diagonal block * -* Multiply by inv(U'(K)), where U(K) is the transformation +* Multiply by inv(U**T(K)), where U(K) is the transformation * stored in column K of A. * CALL SGEMV( 'Transpose', K-1, NRHS, -ONE, B, LDB, A( 1, K ), @@ -211,7 +211,7 @@ * * 2 x 2 diagonal block * -* Multiply by inv(U'(K+1)), where U(K+1) is the transformation +* Multiply by inv(U**T(K+1)), where U(K+1) is the transformation * stored in columns K and K+1 of A. * CALL SGEMV( 'Transpose', K-1, NRHS, -ONE, B, LDB, A( 1, K ), @@ -232,7 +232,7 @@ * ELSE * -* Solve A*X = B, where A = L*D*L'. +* Solve A*X = B, where A = L*D*L**T. * * First solve L*D*X = B, overwriting B with X. * @@ -306,7 +306,7 @@ GO TO 60 80 CONTINUE * -* Next solve L'*X = B, overwriting B with X. +* Next solve L**T *X = B, overwriting B with X. * * K is the main loop index, decreasing from N to 1 in steps of * 1 or 2, depending on the size of the diagonal blocks. @@ -323,7 +323,7 @@ * * 1 x 1 diagonal block * -* Multiply by inv(L'(K)), where L(K) is the transformation +* Multiply by inv(L**T(K)), where L(K) is the transformation * stored in column K of A. * IF( K.LT.N ) @@ -340,7 +340,7 @@ * * 2 x 2 diagonal block * -* Multiply by inv(L'(K-1)), where L(K-1) is the transformation +* Multiply by inv(L**T(K-1)), where L(K-1) is the transformation * stored in columns K-1 and K of A. * IF( K.LT.N ) THEN |