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/dormr3.f | |
| parent | 53b71f5605f83d116ab6bcf477bfb6d2ca757de1 (diff) | |
First pass to homgenize notation for transpose (**T) and conjugate transpose (**H)
Corresponds to bug0024
Diffstat (limited to 'SRC/dormr3.f')
| -rw-r--r-- | SRC/dormr3.f | 18 |
1 files changed, 9 insertions, 9 deletions
diff --git a/SRC/dormr3.f b/SRC/dormr3.f index 94b960fa..cf33553e 100644 --- a/SRC/dormr3.f +++ b/SRC/dormr3.f @@ -21,11 +21,11 @@ * * Q * C if SIDE = 'L' and TRANS = 'N', or * -* Q'* C if SIDE = 'L' and TRANS = 'T', or +* Q**T* C if SIDE = 'L' and TRANS = 'C', or * * C * Q if SIDE = 'R' and TRANS = 'N', or * -* C * Q' if SIDE = 'R' and TRANS = 'T', +* C * Q**T if SIDE = 'R' and TRANS = 'C', * * where Q is a real orthogonal matrix defined as the product of k * elementary reflectors @@ -39,12 +39,12 @@ * ========= * * SIDE (input) CHARACTER*1 -* = 'L': apply Q or Q' from the Left -* = 'R': apply Q or Q' from the Right +* = 'L': apply Q or Q**T from the Left +* = 'R': apply Q or Q**T from the Right * * TRANS (input) CHARACTER*1 * = 'N': apply Q (No transpose) -* = 'T': apply Q' (Transpose) +* = 'T': apply Q**T (Transpose) * * M (input) INTEGER * The number of rows of the matrix C. M >= 0. @@ -80,7 +80,7 @@ * * C (input/output) DOUBLE PRECISION array, dimension (LDC,N) * On entry, the m-by-n matrix C. -* On exit, C is overwritten by Q*C or Q'*C or C*Q' or C*Q. +* On exit, C is overwritten by Q*C or Q**T*C or C*Q**T or C*Q. * * LDC (input) INTEGER * The leading dimension of the array C. LDC >= max(1,M). @@ -181,19 +181,19 @@ DO 10 I = I1, I2, I3 IF( LEFT ) THEN * -* H(i) or H(i)' is applied to C(i:m,1:n) +* H(i) or H(i)**T is applied to C(i:m,1:n) * MI = M - I + 1 IC = I ELSE * -* H(i) or H(i)' is applied to C(1:m,i:n) +* H(i) or H(i)**T is applied to C(1:m,i:n) * NI = N - I + 1 JC = I END IF * -* Apply H(i) or H(i)' +* Apply H(i) or H(i)**T * CALL DLARZ( SIDE, MI, NI, L, A( I, JA ), LDA, TAU( I ), $ C( IC, JC ), LDC, WORK ) |