diff options
author | langou <langou@users.noreply.github.com> | 2009-12-21 21:25:30 +0000 |
---|---|---|
committer | langou <langou@users.noreply.github.com> | 2009-12-21 21:25:30 +0000 |
commit | 4252cf18cecf22991c4fe66bb4db58b9074a851b (patch) | |
tree | 8ad84f2265cfe73814fe77db52edf58835ab9211 | |
parent | f48dd07dbaa629e05ad96ca81b63fb1b896171b7 (diff) | |
download | lapack-4252cf18cecf22991c4fe66bb4db58b9074a851b.tar.gz lapack-4252cf18cecf22991c4fe66bb4db58b9074a851b.tar.bz2 lapack-4252cf18cecf22991c4fe66bb4db58b9074a851b.zip |
Darned!!! I messed up with my commit r709!!!!
This is
" ( 1 - CONJG( TAU ) ) * ( ALPHA ) = ABS( ALPHA ) "
that needs to be true so, indeed, the previous code was correct and r709 is a
mistake. (r710 is good.) So this commit rollbacks r709 and modifies a little
r710.
-rw-r--r-- | SRC/clarfp.f | 4 | ||||
-rw-r--r-- | SRC/zlarfp.f | 4 |
2 files changed, 4 insertions, 4 deletions
diff --git a/SRC/clarfp.f b/SRC/clarfp.f index a00f7ef4..c4522b5e 100644 --- a/SRC/clarfp.f +++ b/SRC/clarfp.f @@ -109,7 +109,7 @@ ELSE * Only "reflecting" the diagonal entry to be real and non-negative. XNORM = SLAPY2( ALPHR, ALPHI ) - TAU = CMPLX( ONE - ALPHR / XNORM, ALPHI / XNORM ) + TAU = CMPLX( ONE - ALPHR / XNORM, -ALPHI / XNORM ) DO J = 1, N-1 X( 1 + (J-1)*INCX ) = ZERO END DO @@ -176,7 +176,7 @@ END IF ELSE XNORM = SLAPY2( ALPHR, ALPHI ) - TAU = CMPLX( ONE - ALPHR / XNORM, ALPHI / XNORM ) + TAU = CMPLX( ONE - ALPHR / XNORM, -ALPHI / XNORM ) DO J = 1, N-1 X( 1 + (J-1)*INCX ) = ZERO END DO diff --git a/SRC/zlarfp.f b/SRC/zlarfp.f index cdb1e3ad..9c40ffc5 100644 --- a/SRC/zlarfp.f +++ b/SRC/zlarfp.f @@ -109,7 +109,7 @@ ELSE * Only "reflecting" the diagonal entry to be real and non-negative. XNORM = DLAPY2( ALPHR, ALPHI ) - TAU = DCMPLX( ONE - ALPHR / XNORM, ALPHI / XNORM ) + TAU = DCMPLX( ONE - ALPHR / XNORM, -ALPHI / XNORM ) DO J = 1, N-1 X( 1 + (J-1)*INCX ) = ZERO END DO @@ -176,7 +176,7 @@ END IF ELSE XNORM = DLAPY2( ALPHR, ALPHI ) - TAU = DCMPLX( ONE - ALPHR / XNORM, ALPHI / XNORM ) + TAU = DCMPLX( ONE - ALPHR / XNORM, -ALPHI / XNORM ) DO J = 1, N-1 X( 1 + (J-1)*INCX ) = ZERO END DO |