Fix bug 139 submitted by Christoph Conrads on Oct 6th 2015
authorjulie <julielangou@users.noreply.github.com>
Thu, 22 Oct 2015 04:39:19 +0000 (04:39 +0000)
committerjulie <julielangou@users.noreply.github.com>
Thu, 22 Oct 2015 04:39:19 +0000 (04:39 +0000)
commit6d3c6b3d719c4fa38f8321e398a8ab142c124f22
treeffda9dbbfb739b23789a2f0e5056e8960acec387
parent2491c9367a9c5e827a935679a4755387b348e159
Fix bug 139 submitted by Christoph Conrads on Oct 6th 2015
"according to the LAPACK documentation, xSTEDC guarantees to
return the eigenvalues in ascending order (parameter D on exit). To
this end, it sorts the eigenvalues if necessary.

If there is only a single subproblem of size n, no sorting algorithm is
called (cf. {s,d}tedc.f, line 450). Furthermore, xLAED0
({s,d}tedc.f:400) does not guarantee to return eigenvalues in ascending
order. Thus, xSTEDC may return eigenvalues that are not in ascending
order."

Applied patch provided by Christoph Conrads: "always sorting the eigenvalues"
SRC/cstedc.f
SRC/dstedc.f
SRC/sstedc.f
SRC/zstedc.f