diff options
author | deaglanhalligan <deaglanhalligan@8a072113-8704-0410-8d35-dd094bca7971> | 2009-02-11 06:06:54 +0000 |
---|---|---|
committer | deaglanhalligan <deaglanhalligan@8a072113-8704-0410-8d35-dd094bca7971> | 2009-02-11 06:06:54 +0000 |
commit | 1f0792306441cdcd175199ad6720671df249dc47 (patch) | |
tree | 3a1c16a498b614128bc3977416ee7b3090d03dae /SRC/zla_hercond_x.f | |
parent | e9454e224e5fa4d61076493be3916a64f274a5ae (diff) | |
download | lapack-1f0792306441cdcd175199ad6720671df249dc47.tar.gz lapack-1f0792306441cdcd175199ad6720671df249dc47.tar.bz2 lapack-1f0792306441cdcd175199ad6720671df249dc47.zip |
Added Jasons patches to bring work and rwork down to 2*N in the complex refinement routines.
Diffstat (limited to 'SRC/zla_hercond_x.f')
-rw-r--r-- | SRC/zla_hercond_x.f | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/SRC/zla_hercond_x.f b/SRC/zla_hercond_x.f index ea031ced..25129156 100644 --- a/SRC/zla_hercond_x.f +++ b/SRC/zla_hercond_x.f @@ -91,7 +91,7 @@ TMP = TMP + CABS1( A( I, J ) * X( J ) ) END IF END DO - RWORK( 2*N+I ) = TMP + RWORK( I ) = TMP ANORM = MAX( ANORM, TMP ) END DO ELSE @@ -104,7 +104,7 @@ TMP = TMP + CABS1( A( I, J ) * X( J ) ) END IF END DO - RWORK( 2*N+I ) = TMP + RWORK( I ) = TMP ANORM = MAX( ANORM, TMP ) END DO END IF @@ -131,7 +131,7 @@ * Multiply by R. * DO I = 1, N - WORK( I ) = WORK( I ) * RWORK( 2*N+I ) + WORK( I ) = WORK( I ) * RWORK( I ) END DO * IF ( UP ) THEN @@ -166,7 +166,7 @@ * Multiply by R. * DO I = 1, N - WORK( I ) = WORK( I ) * RWORK( 2*N+I ) + WORK( I ) = WORK( I ) * RWORK( I ) END DO END IF GO TO 10 |