Fix #81 reported by @reeuwijk-altium

The documentation of 'work' in sgesvj states that it should be of size
max(4,M+N), whereas de documentation of 'lwork' in sgesvj states that
it should be 'length of WORK, WORK >= MAX(6,M+N)'.

diff --git a/SRC/dgesvj.f b/SRC/dgesvj.f
index 2286b7b..3e5387e 100644 (file)
--- a/SRC/dgesvj.f
+++ b/SRC/dgesvj.f
@@ -208,7 +208,7 @@
 *>
 *> \param[in,out] WORK
 *> \verbatim
-*>          WORK is DOUBLE PRECISION array, dimension max(4,M+N).
+*>          WORK is DOUBLE PRECISION array, dimension MAX(6,M+N).
 *>          On entry :
 *>          If JOBU .EQ. 'C' :
 *>          WORK(1) = CTOL, where CTOL defines the threshold for convergence.

diff --git a/SRC/sgesvj.f b/SRC/sgesvj.f
index e45af02..51270d0 100644 (file)
--- a/SRC/sgesvj.f
+++ b/SRC/sgesvj.f
@@ -208,7 +208,7 @@
 *>
 *> \param[in,out] WORK
 *> \verbatim
-*>          WORK is REAL array, dimension max(4,M+N).
+*>          WORK is REAL array, dimension MAX(6,M+N).
 *>          On entry,
 *>          If JOBU .EQ. 'C' :
 *>          WORK(1) = CTOL, where CTOL defines the threshold for convergence.