diff options
| author | Julien Langou <julien.langou@ucdenver.edu> | 2016-12-22 11:28:32 +0100 |
|---|---|---|
| committer | Julien Langou <julien.langou@ucdenver.edu> | 2016-12-22 11:28:32 +0100 |
| commit | 4229de15d062b2e28fa935af186119ff50afda78 (patch) | |
| tree | 78a86238a5cc0d82630f1322267dbe8a782ef46e | |
| parent | 8f51f38d5c651b471474b7dc430613cb088c4f4e (diff) | |
somewhat better description of the minimal workspace needed for a workspace query
| -rw-r--r-- | SRC/cgejsv.f | 6 | ||||
| -rw-r--r-- | SRC/zgejsv.f | 6 |
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
|