summaryrefslogtreecommitdiff
path: root/SRC/dlarzb.f
diff options
context:
space:
mode:
authorjulie <julielangou@users.noreply.github.com>2011-04-02 11:08:56 +0000
committerjulie <julielangou@users.noreply.github.com>2011-04-02 11:08:56 +0000
commitf2953573ede24d7f8c01fdb18de48f65f00a9943 (patch)
tree53172aa9083b9aa1abe2d6c130f7c173d8d8725b /SRC/dlarzb.f
parent53b71f5605f83d116ab6bcf477bfb6d2ca757de1 (diff)
downloadlapack-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/dlarzb.f')
-rw-r--r--SRC/dlarzb.f24
1 files changed, 12 insertions, 12 deletions
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 )