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/dggqrf.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/dggqrf.f')
-rw-r--r-- | SRC/dggqrf.f | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/SRC/dggqrf.f b/SRC/dggqrf.f index edcdb6a7..c4660326 100644 --- a/SRC/dggqrf.f +++ b/SRC/dggqrf.f @@ -40,9 +40,9 @@ * In particular, if B is square and nonsingular, the GQR factorization * of A and B implicitly gives the QR factorization of inv(B)*A: * -* inv(B)*A = Z'*(inv(T)*R) +* inv(B)*A = Z**T*(inv(T)*R) * -* 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 @@ -119,7 +119,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(1:i-1) = 0 and v(i) = 1; v(i+1:n) is stored on exit in A(i+1:n,i), @@ -133,7 +133,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(p-k+i+1:p) = 0 and v(p-k+i) = 1; v(1:p-k+i-1) is stored on exit in @@ -194,7 +194,7 @@ CALL DGEQRF( N, M, A, LDA, TAUA, WORK, LWORK, INFO ) LOPT = WORK( 1 ) * -* Update B := Q'*B. +* Update B := Q**T*B. * CALL DORMQR( 'Left', 'Transpose', N, P, MIN( N, M ), A, LDA, TAUA, $ B, LDB, WORK, LWORK, INFO ) |