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/sggglm.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/sggglm.f')
-rw-r--r-- | SRC/sggglm.f | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/SRC/sggglm.f b/SRC/sggglm.f index a2b107d5..4b5eac46 100644 --- a/SRC/sggglm.f +++ b/SRC/sggglm.f @@ -186,9 +186,9 @@ * * Compute the GQR factorization of matrices A and B: * -* Q'*A = ( R11 ) M, Q'*B*Z' = ( T11 T12 ) M -* ( 0 ) N-M ( 0 T22 ) N-M -* M M+P-N N-M +* Q**T*A = ( R11 ) M, Q**T*B*Z**T = ( T11 T12 ) M +* ( 0 ) N-M ( 0 T22 ) N-M +* M M+P-N N-M * * where R11 and T22 are upper triangular, and Q and Z are * orthogonal. @@ -197,8 +197,8 @@ $ WORK( M+NP+1 ), LWORK-M-NP, INFO ) LOPT = WORK( M+NP+1 ) * -* Update left-hand-side vector d = Q'*d = ( d1 ) M -* ( d2 ) N-M +* Update left-hand-side vector d = Q**T*d = ( d1 ) M +* ( d2 ) N-M * CALL SORMQR( 'Left', 'Transpose', N, 1, M, A, LDA, WORK, D, $ MAX( 1, N ), WORK( M+NP+1 ), LWORK-M-NP, INFO ) @@ -245,7 +245,7 @@ CALL SCOPY( M, D, 1, X, 1 ) END IF * -* Backward transformation y = Z'*y +* Backward transformation y = Z**T *y * CALL SORMRQ( 'Left', 'Transpose', P, 1, NP, $ B( MAX( 1, N-P+1 ), 1 ), LDB, WORK( M+1 ), Y, |