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.

diff --git a/SRC/csyconv.f b/SRC/csyconv.f
index 79a5793..f98c084 100644 (file)
--- 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 (file)
--- 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 (file)
--- 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 (file)
--- 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