Email from Edward Smyth from NAG. (Thanks Ed!)

In calls to {D,S}BDSVDX in the test programs {d,s}chkbd.f, insufficient space
may be provided for the array Z i.e. WORK(IWBZ). The correction is to change

 IWWORK = IWBZ + MNMIN*(MNMIN*2+1)

to

 IWWORK = IWBZ + 2*MNMIN*(MNMIN+1)

on line 1146 in both files.

diff --git a/TESTING/EIG/dchkbd.f b/TESTING/EIG/dchkbd.f
index cf0ad14fd88349574b85b922f7af851503977c25..3a70053ccf9e743c4a7e9363b353cfe7f2e5b4b5 100644 (file)
--- a/TESTING/EIG/dchkbd.f
+++ b/TESTING/EIG/dchkbd.f
@@ -1143,7 +1143,7 @@
             IWBD = IWBS + MNMIN
             IWBE = IWBD + MNMIN
             IWBZ = IWBE + MNMIN
-            IWWORK = IWBZ + MNMIN*(MNMIN*2+1)
+            IWWORK = IWBZ + 2*MNMIN*(MNMIN+1)
             MNMIN2 = MAX( 1,MNMIN*2 )
 *
             CALL DCOPY( MNMIN, BD, 1, WORK( IWBD ), 1 )

diff --git a/TESTING/EIG/schkbd.f b/TESTING/EIG/schkbd.f
index f7e883093e4cc9e3b3329bf7d775f3b68606b4d4..261567ed1fde7eb3f76b749ad3360338a2d0b9bc 100644 (file)
--- a/TESTING/EIG/schkbd.f
+++ b/TESTING/EIG/schkbd.f
@@ -1143,7 +1143,7 @@
             IWBD = IWBS + MNMIN
             IWBE = IWBD + MNMIN
             IWBZ = IWBE + MNMIN
-            IWWORK = IWBZ + MNMIN*(MNMIN*2+1)
+            IWWORK = IWBZ + 2*MNMIN*(MNMIN+1)
             MNMIN2 = MAX( 1,MNMIN*2 )
 *
             CALL SCOPY( MNMIN, BD, 1, WORK( IWBD ), 1 )