First pass to homgenize notation for transpose (**T) and conjugate transpose (**H)

Corresponds to bug0024

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 )