From eedd28ee8d97e407f4b551e92dce4b193bfa674a Mon Sep 17 00:00:00 2001 From: langou Date: Wed, 4 Nov 2015 19:54:31 +0000 Subject: [PATCH] Thanks to Lawrence Mulholland (NAG) for reporting the problem. While the array E should really be of size N-1, xSYCONV() is not ready for this and so, in the comment at the header level, we now require E to be of size N. Also xSYTRS2() is using E(N) = ZERO, so if we make E of size N-1 in xSYCONV(), then we need to see the impact on xSYTRS2(). This is long term cleaning, so for now, we revert back to requesting an array E of size N. Thanks Lawrence for noticing this. --- SRC/csyconv.f | 2 +- SRC/dsyconv.f | 2 +- SRC/ssyconv.f | 2 +- SRC/zsyconv.f | 2 +- 4 files changed, 4 insertions(+), 4 deletions(-) diff --git a/SRC/csyconv.f b/SRC/csyconv.f index 79a5793..f98c084 100644 --- a/SRC/csyconv.f +++ b/SRC/csyconv.f @@ -87,7 +87,7 @@ *> *> \param[out] E *> \verbatim -*> E is COMPLEX array, dimension (N-1) +*> E is COMPLEX array, dimension (N) *> E stores the supdiagonal/subdiagonal of the symmetric 1-by-1 *> or 2-by-2 block diagonal matrix D in LDLT. *> \endverbatim diff --git a/SRC/dsyconv.f b/SRC/dsyconv.f index 768fd51..a08368b 100644 --- a/SRC/dsyconv.f +++ b/SRC/dsyconv.f @@ -87,7 +87,7 @@ *> *> \param[out] E *> \verbatim -*> E is DOUBLE PRECISION array, dimension (N-1) +*> E is DOUBLE PRECISION array, dimension (N) *> E stores the supdiagonal/subdiagonal of the symmetric 1-by-1 *> or 2-by-2 block diagonal matrix D in LDLT. *> \endverbatim diff --git a/SRC/ssyconv.f b/SRC/ssyconv.f index b636037..f7384ed 100644 --- a/SRC/ssyconv.f +++ b/SRC/ssyconv.f @@ -87,7 +87,7 @@ *> *> \param[out] E *> \verbatim -*> E is REAL array, dimension (N-1) +*> E is REAL array, dimension (N) *> E stores the supdiagonal/subdiagonal of the symmetric 1-by-1 *> or 2-by-2 block diagonal matrix D in LDLT. *> \endverbatim diff --git a/SRC/zsyconv.f b/SRC/zsyconv.f index 26d1c27..e62ed01 100644 --- a/SRC/zsyconv.f +++ b/SRC/zsyconv.f @@ -87,7 +87,7 @@ *> *> \param[out] E *> \verbatim -*> E is COMPLEX*16 array, dimension (N-1) +*> E is COMPLEX*16 array, dimension (N) *> E stores the supdiagonal/subdiagonal of the symmetric 1-by-1 *> or 2-by-2 block diagonal matrix D in LDLT. *> \endverbatim -- 2.7.4