From 46e09e3c99a1c6aeedcde4281d7625a76f48d315 Mon Sep 17 00:00:00 2001 From: deaglanhalligan Date: Fri, 13 Feb 2009 23:58:46 +0000 Subject: Restructured loops in {c,z}la_yyrcond_{x,c} and {s,c,d,z}la_yyamv to do less comparisons. --- SRC/zla_porcond_x.f | 24 +++++++++++------------- 1 file changed, 11 insertions(+), 13 deletions(-) (limited to 'SRC/zla_porcond_x.f') diff --git a/SRC/zla_porcond_x.f b/SRC/zla_porcond_x.f index 2692dad6..3b93670e 100644 --- a/SRC/zla_porcond_x.f +++ b/SRC/zla_porcond_x.f @@ -33,7 +33,7 @@ * * WORK COMPLEX*16 workspace of size 2*N. * -* RWORK DOUBLE PRECISION workspace of size 3*N. +* RWORK DOUBLE PRECISION workspace of size N. * * ===================================================================== * @@ -83,12 +83,11 @@ IF ( UP ) THEN DO I = 1, N TMP = 0.0D+0 - DO J = 1, N - IF ( I.GT.J ) THEN - TMP = TMP + CABS1( A( J, I ) * X( J ) ) - ELSE - TMP = TMP + CABS1( A( I, J ) * X( J ) ) - END IF + DO J = 1, I + TMP = TMP + CABS1( A( J, I ) * X( J ) ) + END DO + DO J = I+1, N + TMP = TMP + CABS1( A( I, J ) * X( J ) ) END DO RWORK( I ) = TMP ANORM = MAX( ANORM, TMP ) @@ -96,12 +95,11 @@ ELSE DO I = 1, N TMP = 0.0D+0 - DO J = 1, N - IF ( I.LT.J ) THEN - TMP = TMP + CABS1( A( J, I ) * X( J ) ) - ELSE - TMP = TMP + CABS1( A( I, J ) * X( J ) ) - END IF + DO J = 1, I + TMP = TMP + CABS1( A( I, J ) * X( J ) ) + END DO + DO J = I+1, N + TMP = TMP + CABS1( A( J, I ) * X( J ) ) END DO RWORK( I ) = TMP ANORM = MAX( ANORM, TMP ) -- cgit v1.2.3