Modified {c,z}syequb to use 2*N WORK rather than 3*N.

diff --git a/SRC/csyequb.f b/SRC/csyequb.f
index 90d57d9..2d0315a 100644 (file)
--- a/SRC/csyequb.f
+++ b/SRC/csyequb.f
@@ -185,10 +185,10 @@
         AVG = AVG / N
 
         STD = 0.0
-        DO I = 2*N+1, 3*N
-           WORK( I ) = S( I-2*N ) * WORK( I-2*N ) - AVG
+        DO I = N+1, 2*N
+           WORK( I ) = S( I-N ) * WORK( I-N ) - AVG
         END DO
-        CALL CLASSQ( N, WORK( 2*N+1 ), 1, SCALE, SUMSQ )
+        CALL CLASSQ( N, WORK( N+1 ), 1, SCALE, SUMSQ )
         STD = SCALE * SQRT( SUMSQ / N )
 
         IF ( STD .LT. TOL * AVG ) GOTO 999

diff --git a/SRC/zsyequb.f b/SRC/zsyequb.f
index b46b760..5a749d2 100644 (file)
--- a/SRC/zsyequb.f
+++ b/SRC/zsyequb.f
@@ -185,10 +185,10 @@
         AVG = AVG / N
 
         STD = 0.0D+0
-        DO I = 2*N+1, 3*N
-           WORK( I ) = S( I-2*N ) * WORK( I-2*N ) - AVG
+        DO I = N+1, 2*N
+           WORK( I ) = S( I-N ) * WORK( I-N ) - AVG
         END DO
-        CALL ZLASSQ( N, WORK( 2*N+1 ), 1, SCALE, SUMSQ )
+        CALL ZLASSQ( N, WORK( N+1 ), 1, SCALE, SUMSQ )
         STD = SCALE * SQRT( SUMSQ / N )
 
         IF ( STD .LT. TOL * AVG ) GOTO 999