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/slarft.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/slarft.f')
| -rw-r--r-- | SRC/slarft.f | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/SRC/slarft.f b/SRC/slarft.f index f4eaab07..d2801834 100644 --- a/SRC/slarft.f +++ b/SRC/slarft.f @@ -27,12 +27,12 @@ * If STOREV = 'C', the vector which defines the elementary reflector * H(i) is stored in the i-th column of the array V, and * -* H = I - V * T * V' +* H = I - V * T * V**T * * If STOREV = 'R', the vector which defines the elementary reflector * H(i) is stored in the i-th row of the array V, and * -* H = I - V' * T * V +* H = I - V**T * T * V * * Arguments * ========= @@ -150,7 +150,7 @@ END DO J = MIN( LASTV, PREVLASTV ) * -* T(1:i-1,i) := - tau(i) * V(i:j,1:i-1)' * V(i:j,i) +* T(1:i-1,i) := - tau(i) * V(i:j,1:i-1)**T * V(i:j,i) * CALL SGEMV( 'Transpose', J-I+1, I-1, -TAU( I ), $ V( I, 1 ), LDV, V( I, I ), 1, ZERO, @@ -162,7 +162,7 @@ END DO J = MIN( LASTV, PREVLASTV ) * -* T(1:i-1,i) := - tau(i) * V(1:i-1,i:j) * V(i,i:j)' +* T(1:i-1,i) := - tau(i) * V(1:i-1,i:j) * V(i,i:j)**T * CALL SGEMV( 'No transpose', I-1, J-I+1, -TAU( I ), $ V( 1, I ), LDV, V( I, I ), LDV, ZERO, @@ -207,7 +207,7 @@ J = MAX( LASTV, PREVLASTV ) * * T(i+1:k,i) := -* - tau(i) * V(j:n-k+i,i+1:k)' * V(j:n-k+i,i) +* - tau(i) * V(j:n-k+i,i+1:k)**T * V(j:n-k+i,i) * CALL SGEMV( 'Transpose', N-K+I-J+1, K-I, -TAU( I ), $ V( J, I+1 ), LDV, V( J, I ), 1, ZERO, @@ -223,7 +223,7 @@ J = MAX( LASTV, PREVLASTV ) * * T(i+1:k,i) := -* - tau(i) * V(i+1:k,j:n-k+i) * V(i,j:n-k+i)' +* - tau(i) * V(i+1:k,j:n-k+i) * V(i,j:n-k+i)**T * CALL SGEMV( 'No transpose', K-I, N-K+I-J+1, $ -TAU( I ), V( I+1, J ), LDV, V( I, J ), LDV, |