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/zlarz.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/zlarz.f')
| -rw-r--r-- | SRC/zlarz.f | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/SRC/zlarz.f b/SRC/zlarz.f index d235191f..f2a90d04 100644 --- a/SRC/zlarz.f +++ b/SRC/zlarz.f @@ -21,13 +21,13 @@ * M-by-N 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. * * H is a product of k elementary reflectors as returned by ZTZRZF. @@ -105,7 +105,7 @@ CALL ZCOPY( N, C, LDC, WORK, 1 ) CALL ZLACGV( N, WORK, 1 ) * -* w( 1:n ) = conjg( w( 1:n ) + C( m-l+1:m, 1:n )' * v( 1:l ) ) +* w( 1:n ) = conjg( w( 1:n ) + C( m-l+1:m, 1:n )**H * v( 1:l ) ) * CALL ZGEMV( 'Conjugate transpose', L, N, ONE, C( M-L+1, 1 ), $ LDC, V, INCV, ONE, WORK, 1 ) @@ -116,7 +116,7 @@ CALL ZAXPY( N, -TAU, WORK, 1, C, LDC ) * * C( m-l+1:m, 1:n ) = C( m-l+1:m, 1:n ) - ... -* tau * v( 1:l ) * conjg( w( 1:n )' ) +* tau * v( 1:l ) * w( 1:n )**H * CALL ZGERU( L, N, -TAU, V, INCV, WORK, 1, C( M-L+1, 1 ), $ LDC ) @@ -142,7 +142,7 @@ CALL ZAXPY( M, -TAU, WORK, 1, C, 1 ) * * C( 1:m, n-l+1:n ) = C( 1:m, n-l+1:n ) - ... -* tau * w( 1:m ) * v( 1:l )' +* tau * w( 1:m ) * v( 1:l )**H * CALL ZGERC( M, L, -TAU, WORK, 1, V, INCV, C( 1, N-L+1 ), $ LDC ) |