From 932dcade9c3a76520522f8a0d2ad35893fc6a849 Mon Sep 17 00:00:00 2001 From: julie Date: Mon, 31 Jan 2011 16:14:51 +0000 Subject: 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. --- SRC/cuncsd.f | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'SRC/cuncsd.f') diff --git a/SRC/cuncsd.f b/SRC/cuncsd.f index ff224cc3..070b45cc 100644 --- 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. -- cgit v1.2.3