From bd47060bcb3a470520622de69ac1426ca4186f5e Mon Sep 17 00:00:00 2001 From: Julie Date: Wed, 2 Nov 2016 18:18:22 -0700 Subject: [PATCH] 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)'. --- SRC/dgesvj.f | 2 +- SRC/sgesvj.f | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/SRC/dgesvj.f b/SRC/dgesvj.f index 2286b7b..3e5387e 100644 --- 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 --- 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. -- 2.7.4