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/sggrqf.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/sggrqf.f')
| -rw-r--r-- | SRC/sggrqf.f | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/SRC/sggrqf.f b/SRC/sggrqf.f index c752c7b0..aed20ed8 100644 --- a/SRC/sggrqf.f +++ b/SRC/sggrqf.f @@ -40,9 +40,9 @@ * In particular, if B is square and nonsingular, the GRQ factorization * of A and B implicitly gives the RQ factorization of A*inv(B): * -* A*inv(B) = (R*inv(T))*Z' +* A*inv(B) = (R*inv(T))*Z**T * -* where inv(B) denotes the inverse of the matrix B, and Z' denotes the +* where inv(B) denotes the inverse of the matrix B, and Z**T denotes the * transpose of the matrix Z. * * Arguments @@ -118,7 +118,7 @@ * * Each H(i) has the form * -* H(i) = I - taua * v * v' +* H(i) = I - taua * v * v**T * * where taua is a real scalar, and v is a real vector with * v(n-k+i+1:n) = 0 and v(n-k+i) = 1; v(1:n-k+i-1) is stored on exit in @@ -132,7 +132,7 @@ * * Each H(i) has the form * -* H(i) = I - taub * v * v' +* H(i) = I - taub * v * v**T * * where taub is a real scalar, and v is a real vector with * v(1:i-1) = 0 and v(i) = 1; v(i+1:p) is stored on exit in B(i+1:p,i), @@ -193,7 +193,7 @@ CALL SGERQF( M, N, A, LDA, TAUA, WORK, LWORK, INFO ) LOPT = WORK( 1 ) * -* Update B := B*Q' +* Update B := B*Q**T * CALL SORMRQ( 'Right', 'Transpose', P, N, MIN( M, N ), $ A( MAX( 1, M-N+1 ), 1 ), LDA, TAUA, B, LDB, WORK, |