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/zgebd2.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/zgebd2.f')
-rw-r--r-- | SRC/zgebd2.f | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/SRC/zgebd2.f b/SRC/zgebd2.f index 5effb30f..2d2db958 100644 --- a/SRC/zgebd2.f +++ b/SRC/zgebd2.f @@ -17,7 +17,7 @@ * ======= * * ZGEBD2 reduces a complex general m by n matrix A to upper or lower -* real bidiagonal form B by a unitary transformation: Q' * A * P = B. +* real bidiagonal form B by a unitary transformation: Q**H * A * P = B. * * If m >= n, B is upper bidiagonal; if m < n, B is lower bidiagonal. * @@ -87,7 +87,7 @@ * * Each H(i) and G(i) has the form: * -* H(i) = I - tauq * v * v' and G(i) = I - taup * u * u' +* H(i) = I - tauq * v * v**H and G(i) = I - taup * u * u**H * * where tauq and taup are complex scalars, and v and u are complex * vectors; v(1:i-1) = 0, v(i) = 1, and v(i+1:m) is stored on exit in @@ -100,7 +100,7 @@ * * Each H(i) and G(i) has the form: * -* H(i) = I - tauq * v * v' and G(i) = I - taup * u * u' +* H(i) = I - tauq * v * v**H and G(i) = I - taup * u * u**H * * where tauq and taup are complex scalars, v and u are complex vectors; * v(1:i) = 0, v(i+1) = 1, and v(i+2:m) is stored on exit in A(i+2:m,i); @@ -170,7 +170,7 @@ D( I ) = ALPHA A( I, I ) = ONE * -* Apply H(i)' to A(i:m,i+1:n) from the left +* Apply H(i)**H to A(i:m,i+1:n) from the left * IF( I.LT.N ) $ CALL ZLARF( 'Left', M-I+1, N-I, A( I, I ), 1, @@ -233,7 +233,7 @@ E( I ) = ALPHA A( I+1, I ) = ONE * -* Apply H(i)' to A(i+1:m,i+1:n) from the left +* Apply H(i)**H to A(i+1:m,i+1:n) from the left * CALL ZLARF( 'Left', M-I, N-I, A( I+1, I ), 1, $ DCONJG( TAUQ( I ) ), A( I+1, I+1 ), LDA, |