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/dsptrd.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/dsptrd.f')
| -rw-r--r-- | SRC/dsptrd.f | 16 |
1 files changed, 8 insertions, 8 deletions
diff --git a/SRC/dsptrd.f b/SRC/dsptrd.f index f91adbdc..331dd2af 100644 --- a/SRC/dsptrd.f +++ b/SRC/dsptrd.f @@ -73,7 +73,7 @@ * * Each H(i) has the form * -* H(i) = I - tau * v * v' +* H(i) = I - tau * v * v**T * * where tau is a real scalar, and v is a real vector with * v(i+1:n) = 0 and v(i) = 1; v(1:i-1) is stored on exit in AP, @@ -86,7 +86,7 @@ * * Each H(i) has the form * -* H(i) = I - tau * v * v' +* H(i) = I - tau * v * v**T * * where tau is a real scalar, and v is a real vector with * v(1:i) = 0 and v(i+1) = 1; v(i+2:n) is stored on exit in AP, @@ -141,7 +141,7 @@ I1 = N*( N-1 ) / 2 + 1 DO 10 I = N - 1, 1, -1 * -* Generate elementary reflector H(i) = I - tau * v * v' +* Generate elementary reflector H(i) = I - tau * v * v**T * to annihilate A(1:i-1,i+1) * CALL DLARFG( I, AP( I1+I-1 ), AP( I1 ), 1, TAUI ) @@ -158,13 +158,13 @@ CALL DSPMV( UPLO, I, TAUI, AP, AP( I1 ), 1, ZERO, TAU, $ 1 ) * -* Compute w := y - 1/2 * tau * (y'*v) * v +* Compute w := y - 1/2 * tau * (y**T *v) * v * ALPHA = -HALF*TAUI*DDOT( I, TAU, 1, AP( I1 ), 1 ) CALL DAXPY( I, ALPHA, AP( I1 ), 1, TAU, 1 ) * * Apply the transformation as a rank-2 update: -* A := A - v * w' - w * v' +* A := A - v * w**T - w * v**T * CALL DSPR2( UPLO, I, -ONE, AP( I1 ), 1, TAU, 1, AP ) * @@ -184,7 +184,7 @@ DO 20 I = 1, N - 1 I1I1 = II + N - I + 1 * -* Generate elementary reflector H(i) = I - tau * v * v' +* Generate elementary reflector H(i) = I - tau * v * v**T * to annihilate A(i+2:n,i) * CALL DLARFG( N-I, AP( II+1 ), AP( II+2 ), 1, TAUI ) @@ -201,14 +201,14 @@ CALL DSPMV( UPLO, N-I, TAUI, AP( I1I1 ), AP( II+1 ), 1, $ ZERO, TAU( I ), 1 ) * -* Compute w := y - 1/2 * tau * (y'*v) * v +* Compute w := y - 1/2 * tau * (y**T *v) * v * ALPHA = -HALF*TAUI*DDOT( N-I, TAU( I ), 1, AP( II+1 ), $ 1 ) CALL DAXPY( N-I, ALPHA, AP( II+1 ), 1, TAU( I ), 1 ) * * Apply the transformation as a rank-2 update: -* A := A - v * w' - w * v' +* A := A - v * w' - w * v**T * CALL DSPR2( UPLO, N-I, -ONE, AP( II+1 ), 1, TAU( I ), 1, $ AP( I1I1 ) ) |