From f2953573ede24d7f8c01fdb18de48f65f00a9943 Mon Sep 17 00:00:00 2001 From: julie Date: Sat, 2 Apr 2011 11:08:56 +0000 Subject: First pass to homgenize notation for transpose (**T) and conjugate transpose (**H) Corresponds to bug0024 --- SRC/dggqrf.f | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) (limited to 'SRC/dggqrf.f') 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 ) -- cgit v1.2.3