diff options
| -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 |