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/clarf.f | |
| parent | 53b71f5605f83d116ab6bcf477bfb6d2ca757de1 (diff) | |
| download | lapack-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/clarf.f')
| -rw-r--r-- | SRC/clarf.f | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/SRC/clarf.f b/SRC/clarf.f index e4ee75ca..c9fd3ef0 100644 --- a/SRC/clarf.f +++ b/SRC/clarf.f @@ -22,13 +22,13 @@ * matrix C, from either the left or the right. H is represented in the * form * -* H = I - tau * v * v' +* H = I - tau * v * v**H * * where tau is a complex scalar and v is a complex vector. * * If tau = 0, then H is taken to be the unit matrix. * -* To apply H' (the conjugate transpose of H), supply conjg(tau) instead +* To apply H**H (the conjugate transpose of H), supply conjg(tau) instead * tau. * * Arguments @@ -126,12 +126,12 @@ * IF( LASTV.GT.0 ) THEN * -* w(1:lastc,1) := C(1:lastv,1:lastc)' * v(1:lastv,1) +* w(1:lastc,1) := C(1:lastv,1:lastc)**H * v(1:lastv,1) * CALL CGEMV( 'Conjugate transpose', LASTV, LASTC, ONE, $ C, LDC, V, INCV, ZERO, WORK, 1 ) * -* C(1:lastv,1:lastc) := C(...) - v(1:lastv,1) * w(1:lastc,1)' +* C(1:lastv,1:lastc) := C(...) - v(1:lastv,1) * w(1:lastc,1)**H * CALL CGERC( LASTV, LASTC, -TAU, V, INCV, WORK, 1, C, LDC ) END IF @@ -146,7 +146,7 @@ CALL CGEMV( 'No transpose', LASTC, LASTV, ONE, C, LDC, $ V, INCV, ZERO, WORK, 1 ) * -* C(1:lastc,1:lastv) := C(...) - w(1:lastc,1) * v(1:lastv,1)' +* C(1:lastc,1:lastv) := C(...) - w(1:lastc,1) * v(1:lastv,1)**H * CALL CGERC( LASTC, LASTV, -TAU, WORK, 1, V, INCV, C, LDC ) END IF |