LAPACK Archives

[Lapack] Bug: xSTEDC eigenvalues not sorted

Dear Christoph
Thank you for the bug report.
Your bug is listed as bug 139 in
We will take a look at it very shortly.
On Oct 6, 2015, at 4:20 AM, Christoph Conrads <cconrads@Domain.Removed> wrote:


according to the LAPACK Doxygen 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

As a fix, I propose always sorting the eigenvalues (patch attached).

The bug can be triggered in practice (see compressed file in matrix
market format).

Christoph Conrads
Lapack mailing list

<Prev in Thread] Current Thread [Next in Thread>

For additional information you may use the LAPACK/ScaLAPACK Forum.
Or one of the mailing lists, or