Reduce workspace size of xSYEQUB, xHEEQUB

- give correct minimum workspace size in the parameter description
- {s,d}SYEQUB, xHEEQUB: use WORK(N+1:2*N) instead of WORK(2N+1:3*N)

diff --git a/SRC/cheequb.f b/SRC/cheequb.f
index f5a721a..59b32c0 100644 (file)
--- a/SRC/cheequb.f
+++ b/SRC/cheequb.f
@@ -98,7 +98,7 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is COMPLEX array, dimension (3*N)
+*>          WORK is COMPLEX array, dimension (2*N)
 *> \endverbatim
 *>
 *> \param[out] INFO
@@ -271,10 +271,10 @@
          AVG = AVG / N
 
          STD = 0.0E0
-         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/csyequb.f b/SRC/csyequb.f
index 81a97cd..7b6de6f 100644 (file)
--- a/SRC/csyequb.f
+++ b/SRC/csyequb.f
@@ -98,7 +98,7 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is COMPLEX array, dimension (3*N)
+*>          WORK is COMPLEX array, dimension (2*N)
 *> \endverbatim
 *>
 *> \param[out] INFO

diff --git a/SRC/dsyequb.f b/SRC/dsyequb.f
index 3a5c87a..7a54950 100644 (file)
--- a/SRC/dsyequb.f
+++ b/SRC/dsyequb.f
@@ -97,7 +97,7 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is DOUBLE PRECISION array, dimension (3*N)
+*>          WORK is DOUBLE PRECISION array, dimension (2*N)
 *> \endverbatim
 *>
 *> \param[out] INFO
@@ -262,10 +262,10 @@
          AVG = AVG / N
 
          STD = 0.0D0
-         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 DLASSQ( N, WORK( 2*N+1 ), 1, SCALE, SUMSQ )
+         CALL DLASSQ( N, WORK( N+1 ), 1, SCALE, SUMSQ )
          STD = SCALE * SQRT( SUMSQ / N )
 
          IF ( STD .LT. TOL * AVG ) GOTO 999

diff --git a/SRC/ssyequb.f b/SRC/ssyequb.f
index 7316e92..56dfca8 100644 (file)
--- a/SRC/ssyequb.f
+++ b/SRC/ssyequb.f
@@ -97,7 +97,7 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is REAL array, dimension (3*N)
+*>          WORK is REAL array, dimension (2*N)
 *> \endverbatim
 *>
 *> \param[out] INFO
@@ -262,10 +262,10 @@
          AVG = AVG / N
 
          STD = 0.0E0
-         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 SLASSQ( N, WORK( 2*N+1 ), 1, SCALE, SUMSQ )
+         CALL SLASSQ( N, WORK( N+1 ), 1, SCALE, SUMSQ )
          STD = SCALE * SQRT( SUMSQ / N )
 
          IF ( STD .LT. TOL * AVG ) GOTO 999

diff --git a/SRC/zheequb.f b/SRC/zheequb.f
index eee3fdb..1e06299 100644 (file)
--- a/SRC/zheequb.f
+++ b/SRC/zheequb.f
@@ -98,7 +98,7 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is COMPLEX*16 array, dimension (3*N)
+*>          WORK is COMPLEX*16 array, dimension (2*N)
 *> \endverbatim
 *>
 *> \param[out] INFO
@@ -271,10 +271,10 @@
          AVG = AVG / N
 
          STD = 0.0D0
-         DO I = 2*N+1, 3*N
-            WORK( I ) = S( I-2*N ) * WORK( I-2*N ) - AVG
+         DO I = N+1, 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

diff --git a/SRC/zsyequb.f b/SRC/zsyequb.f
index 7f2e1e8..847f506 100644 (file)
--- a/SRC/zsyequb.f
+++ b/SRC/zsyequb.f
@@ -98,7 +98,7 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is COMPLEX*16 array, dimension (3*N)
+*>          WORK is COMPLEX*16 array, dimension (2*N)
 *> \endverbatim
 *>
 *> \param[out] INFO