Fix bug0070 : IWORK dimension is incorrect

  bug report sent by nmozarto on Jan 27th
  see forum topic 2156

From Brian:
 The correct dimension is (M-R), in which R is the smallest of P, M-P, Q, and M-Q.
 An M-by-M matrix is partitioned into a 2-by-2 block structure. The
 dimensions of the blocks are P-by-Q, P-by-(M-Q), (M-P)-by-Q, and
 (M-P)-by-(M-Q). IWORK is involved in the permutation of rows or
 columns of these blocks, and it is the largest block, in terms of
 number of rows or columns, that determines the size of IWORK.

diff --git a/SRC/cuncsd.f b/SRC/cuncsd.f
index ff224cc3603b29f0f2fe91d589460d419d0bed06..070b45cca235696b6c51c4af946f163e91576871 100644 (file)
--- a/SRC/cuncsd.f
+++ b/SRC/cuncsd.f
@@ -157,7 +157,7 @@
 *          this value as the first entry of the work array, and no error
 *          message related to LRWORK is issued by XERBLA.
 *
-*  IWORK   (workspace) INTEGER array, dimension (M-Q)
+*  IWORK   (workspace) INTEGER array, dimension (M-MIN(P, M-P, Q, M-Q))
 *
 *  INFO    (output) INTEGER
 *          = 0:  successful exit.

diff --git a/SRC/dorcsd.f b/SRC/dorcsd.f
index e845acd038a6aad96ad6ea2c1327fdfaf40fbabc..d1c4377c77c9e4d571fc5358b83a90f807599aa5 100644 (file)
--- a/SRC/dorcsd.f
+++ b/SRC/dorcsd.f
@@ -144,7 +144,7 @@
 *          this value as the first entry of the work array, and no error
 *          message related to LWORK is issued by XERBLA.
 *
-*  IWORK   (workspace) INTEGER array, dimension (M-Q)
+*  IWORK   (workspace) INTEGER array, dimension (M-MIN(P, M-P, Q, M-Q))
 *
 *  INFO    (output) INTEGER
 *          = 0:  successful exit.

diff --git a/SRC/sorcsd.f b/SRC/sorcsd.f
index d78051d6cd45a97874cb3acb5de2d2d085277bd7..36809c2447ae7de97c60374d12f9058fd5cc3f43 100644 (file)
--- a/SRC/sorcsd.f
+++ b/SRC/sorcsd.f
@@ -144,7 +144,7 @@
 *          this value as the first entry of the work array, and no error
 *          message related to LWORK is issued by XERBLA.
 *
-*  IWORK   (workspace) INTEGER array, dimension (M-Q)
+*  IWORK   (workspace) INTEGER array, dimension (M-MIN(P, M-P, Q, M-Q))
 *
 *  INFO    (output) INTEGER
 *          = 0:  successful exit.

diff --git a/SRC/zuncsd.f b/SRC/zuncsd.f
index 6694f4e303f69836fedf5c21c34894f75c36b4e1..1d8b87b1d7217f06cbeb52bcdae6169278cf60f7 100644 (file)
--- a/SRC/zuncsd.f
+++ b/SRC/zuncsd.f
@@ -157,7 +157,7 @@
 *          this value as the first entry of the work array, and no error
 *          message related to LRWORK is issued by XERBLA.
 *
-*  IWORK   (workspace) INTEGER array, dimension (M-Q)
+*  IWORK   (workspace) INTEGER array, dimension (M-MIN(P, M-P, Q, M-Q))
 *
 *  INFO    (output) INTEGER
 *          = 0:  successful exit.