From f2953573ede24d7f8c01fdb18de48f65f00a9943 Mon Sep 17 00:00:00 2001
From: julie <julielangou@users.noreply.github.com>
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/dlarzb.f | 24 ++++++++++++------------
 1 file changed, 12 insertions(+), 12 deletions(-)

(limited to 'SRC/dlarzb.f')

diff --git a/SRC/dlarzb.f b/SRC/dlarzb.f
index 9fc839d4..fbe5a9a0 100644
--- a/SRC/dlarzb.f
+++ b/SRC/dlarzb.f
@@ -27,12 +27,12 @@
 *  =========
 *
 *  SIDE    (input) CHARACTER*1
-*          = 'L': apply H or H' from the Left
-*          = 'R': apply H or H' from the Right
+*          = 'L': apply H or H**T from the Left
+*          = 'R': apply H or H**T from the Right
 *
 *  TRANS   (input) CHARACTER*1
 *          = 'N': apply H (No transpose)
-*          = 'C': apply H' (Transpose)
+*          = 'C': apply H**T (Transpose)
 *
 *  DIRECT  (input) CHARACTER*1
 *          Indicates how H is formed from a product of elementary
@@ -77,7 +77,7 @@
 *
 *  C       (input/output) DOUBLE PRECISION array, dimension (LDC,N)
 *          On entry, the M-by-N matrix C.
-*          On exit, C is overwritten by H*C or H'*C or C*H or C*H'.
+*          On exit, C is overwritten by H*C or H**T*C or C*H or C*H**T.
 *
 *  LDC     (input) INTEGER
 *          The leading dimension of the array C. LDC >= max(1,M).
@@ -140,16 +140,16 @@
 *
       IF( LSAME( SIDE, 'L' ) ) THEN
 *
-*        Form  H * C  or  H' * C
+*        Form  H * C  or  H**T * C
 *
-*        W( 1:n, 1:k ) = C( 1:k, 1:n )'
+*        W( 1:n, 1:k ) = C( 1:k, 1:n )**T
 *
          DO 10 J = 1, K
             CALL DCOPY( N, C( J, 1 ), LDC, WORK( 1, J ), 1 )
    10    CONTINUE
 *
 *        W( 1:n, 1:k ) = W( 1:n, 1:k ) + ...
-*                        C( m-l+1:m, 1:n )' * V( 1:k, 1:l )'
+*                        C( m-l+1:m, 1:n )**T * V( 1:k, 1:l )**T
 *
          IF( L.GT.0 )
      $      CALL DGEMM( 'Transpose', 'Transpose', N, K, L, ONE,
@@ -160,7 +160,7 @@
          CALL DTRMM( 'Right', 'Lower', TRANST, 'Non-unit', N, K, ONE, T,
      $               LDT, WORK, LDWORK )
 *
-*        C( 1:k, 1:n ) = C( 1:k, 1:n ) - W( 1:n, 1:k )'
+*        C( 1:k, 1:n ) = C( 1:k, 1:n ) - W( 1:n, 1:k )**T
 *
          DO 30 J = 1, N
             DO 20 I = 1, K
@@ -169,7 +169,7 @@
    30    CONTINUE
 *
 *        C( m-l+1:m, 1:n ) = C( m-l+1:m, 1:n ) - ...
-*                            V( 1:k, 1:l )' * W( 1:n, 1:k )'
+*                            V( 1:k, 1:l )**T * W( 1:n, 1:k )**T
 *
          IF( L.GT.0 )
      $      CALL DGEMM( 'Transpose', 'Transpose', L, N, K, -ONE, V, LDV,
@@ -177,7 +177,7 @@
 *
       ELSE IF( LSAME( SIDE, 'R' ) ) THEN
 *
-*        Form  C * H  or  C * H'
+*        Form  C * H  or  C * H**T
 *
 *        W( 1:m, 1:k ) = C( 1:m, 1:k )
 *
@@ -186,13 +186,13 @@
    40    CONTINUE
 *
 *        W( 1:m, 1:k ) = W( 1:m, 1:k ) + ...
-*                        C( 1:m, n-l+1:n ) * V( 1:k, 1:l )'
+*                        C( 1:m, n-l+1:n ) * V( 1:k, 1:l )**T
 *
          IF( L.GT.0 )
      $      CALL DGEMM( 'No transpose', 'Transpose', M, K, L, ONE,
      $                  C( 1, N-L+1 ), LDC, V, LDV, ONE, WORK, LDWORK )
 *
-*        W( 1:m, 1:k ) = W( 1:m, 1:k ) * T  or  W( 1:m, 1:k ) * T'
+*        W( 1:m, 1:k ) = W( 1:m, 1:k ) * T  or  W( 1:m, 1:k ) * T**T
 *
          CALL DTRMM( 'Right', 'Lower', TRANS, 'Non-unit', M, K, ONE, T,
      $               LDT, WORK, LDWORK )
-- 
cgit v1.2.3