diff options
author | julie <julielangou@users.noreply.github.com> | 2011-01-31 16:14:51 +0000 |
---|---|---|
committer | julie <julielangou@users.noreply.github.com> | 2011-01-31 16:14:51 +0000 |
commit | 932dcade9c3a76520522f8a0d2ad35893fc6a849 (patch) | |
tree | 0cbbb626e1664d1c07e9306e85a39801448e18a9 /SRC/cuncsd.f | |
parent | 318f35554ef359cb5caafd5fcfcec4d30ce030d8 (diff) | |
download | lapack-932dcade9c3a76520522f8a0d2ad35893fc6a849.tar.gz lapack-932dcade9c3a76520522f8a0d2ad35893fc6a849.tar.bz2 lapack-932dcade9c3a76520522f8a0d2ad35893fc6a849.zip |
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.
Diffstat (limited to 'SRC/cuncsd.f')
-rw-r--r-- | SRC/cuncsd.f | 2 |
1 files changed, 1 insertions, 1 deletions
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. |