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/dgglse.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/dgglse.f')
| -rw-r--r-- | SRC/dgglse.f | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/SRC/dgglse.f b/SRC/dgglse.f index 8b80ddeb..d10c75ad 100644 --- a/SRC/dgglse.f +++ b/SRC/dgglse.f @@ -183,9 +183,9 @@ * * Compute the GRQ factorization of matrices B and A: * -* B*Q' = ( 0 T12 ) P Z'*A*Q' = ( R11 R12 ) N-P -* N-P P ( 0 R22 ) M+P-N -* N-P P +* B*Q**T = ( 0 T12 ) P Z**T*A*Q**T = ( R11 R12 ) N-P +* N-P P ( 0 R22 ) M+P-N +* N-P P * * where T12 and R11 are upper triangular, and Q and Z are * orthogonal. @@ -194,8 +194,8 @@ $ WORK( P+MN+1 ), LWORK-P-MN, INFO ) LOPT = WORK( P+MN+1 ) * -* Update c = Z'*c = ( c1 ) N-P -* ( c2 ) M+P-N +* Update c = Z**T *c = ( c1 ) N-P +* ( c2 ) M+P-N * CALL DORMQR( 'Left', 'Transpose', M, 1, MN, A, LDA, WORK( P+1 ), $ C, MAX( 1, M ), WORK( P+MN+1 ), LWORK-P-MN, INFO ) @@ -254,7 +254,7 @@ CALL DAXPY( NR, -ONE, D, 1, C( N-P+1 ), 1 ) END IF * -* Backward transformation x = Q'*x +* Backward transformation x = Q**T*x * CALL DORMRQ( 'Left', 'Transpose', N, 1, P, B, LDB, WORK( 1 ), X, $ N, WORK( P+MN+1 ), LWORK-P-MN, INFO ) |