diff options
author | langou <langou@users.noreply.github.com> | 2015-11-04 19:54:31 +0000 |
---|---|---|
committer | langou <langou@users.noreply.github.com> | 2015-11-04 19:54:31 +0000 |
commit | eedd28ee8d97e407f4b551e92dce4b193bfa674a (patch) | |
tree | e774385f027544c871a65c957b6dfd929421fd57 /SRC/csyconv.f | |
parent | 54547abdff9232c9d7cf1744a6e3b8f26be94543 (diff) | |
download | lapack-eedd28ee8d97e407f4b551e92dce4b193bfa674a.tar.gz lapack-eedd28ee8d97e407f4b551e92dce4b193bfa674a.tar.bz2 lapack-eedd28ee8d97e407f4b551e92dce4b193bfa674a.zip |
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.
Diffstat (limited to 'SRC/csyconv.f')
-rw-r--r-- | SRC/csyconv.f | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/SRC/csyconv.f b/SRC/csyconv.f index 79a5793f..f98c084d 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 |