summaryrefslogtreecommitdiff
path: root/SRC
diff options
context:
space:
mode:
authorJulien Langou <julien.langou@ucdenver.edu>2016-12-22 11:28:32 +0100
committerJulien Langou <julien.langou@ucdenver.edu>2016-12-22 11:28:32 +0100
commit4229de15d062b2e28fa935af186119ff50afda78 (patch)
tree78a86238a5cc0d82630f1322267dbe8a782ef46e /SRC
parent8f51f38d5c651b471474b7dc430613cb088c4f4e (diff)
downloadlapack-4229de15d062b2e28fa935af186119ff50afda78.tar.gz
lapack-4229de15d062b2e28fa935af186119ff50afda78.tar.bz2
lapack-4229de15d062b2e28fa935af186119ff50afda78.zip
somewhat better description of the minimal workspace needed for a workspace query
Diffstat (limited to 'SRC')
-rw-r--r--SRC/cgejsv.f6
-rw-r--r--SRC/zgejsv.f6
2 files changed, 6 insertions, 6 deletions
diff --git a/SRC/cgejsv.f b/SRC/cgejsv.f
index 02794332..ceff5e6b 100644
--- a/SRC/cgejsv.f
+++ b/SRC/cgejsv.f
@@ -270,7 +270,7 @@
*>
*> \param[out] CWORK
*> \verbatim
-*> CWORK is COMPLEX array, dimension at least LWORK.
+*> CWORK is COMPLEX array, dimension (MAX(2,LWORK))
*> If the call to CGEJSV is a workspace query (indicated by LWORK=-1 or
*> LRWORK=-1), then on exit CWORK(1) contains the required length of
*> CWORK for the job parameters used in the call.
@@ -350,7 +350,7 @@
*>
*> \param[out] RWORK
*> \verbatim
-*> RWORK is REAL array, dimension at least LRWORK.
+*> RWORK is REAL array, dimension (MAX(7,LWORK))
*> On exit,
*> RWORK(1) = Determines the scaling factor SCALE = RWORK(2) / RWORK(1)
*> such that SCALE*SVA(1:N) are the computed singular values
@@ -959,7 +959,7 @@
* Quick return for void matrix (Y3K safe)
* #:)
IF ( ( M .EQ. 0 ) .OR. ( N .EQ. 0 ) ) THEN
- IWORK(1:3) = 0
+ IWORK(1:4) = 0
RWORK(1:7) = 0
RETURN
ENDIF
diff --git a/SRC/zgejsv.f b/SRC/zgejsv.f
index f8b4ba9a..ebc8183f 100644
--- a/SRC/zgejsv.f
+++ b/SRC/zgejsv.f
@@ -270,7 +270,7 @@
*>
*> \param[out] CWORK
*> \verbatim
-*> CWORK is COMPLEX*16 array, dimension at least LWORK.
+*> CWORK is COMPLEX*16 array, dimension (MAX(2,LWORK))
*> If the call to ZGEJSV is a workspace query (indicated by LWORK=-1 or
*> LRWORK=-1), then on exit CWORK(1) contains the required length of
*> CWORK for the job parameters used in the call.
@@ -349,7 +349,7 @@
*>
*> \param[out] RWORK
*> \verbatim
-*> RWORK is DOUBLE PRECISION array, dimension at least LRWORK.
+*> RWORK is DOUBLE PRECISION array, dimension (MAX(7,LWORK))
*> On exit,
*> RWORK(1) = Determines the scaling factor SCALE = RWORK(2) / RWORK(1)
*> such that SCALE*SVA(1:N) are the computed singular values
@@ -962,7 +962,7 @@
* Quick return for void matrix (Y3K safe)
* #:)
IF ( ( M .EQ. 0 ) .OR. ( N .EQ. 0 ) ) THEN
- IWORK(1:3) = 0
+ IWORK(1:4) = 0
RWORK(1:7) = 0
RETURN
ENDIF