Ondra,
Thanks for the bug report. I have added your bug report to the LAPACK
Errata file. See:
http://www.netlib.org/lapack/Errata/ , bug0048.
I confirmed the problem and we will fix it before next release.
Best wishes,
Julie
On Mar 2, 2010, at 1:24 AM, Ondra Kamenik wrote:
Dear lapack developers,
I believe I found a minor little small bug in xDGGES, which should
return info=n+3 if xTGSEN fails to swap eigenvalues at all, and
return info=n+2 if swapping slightly changes eigenvalues so that the
selected eigenvalues are not only in the upper part.
I am referening to dgges.f in lapack 3.2.1 which i downloaded recently.
Indeed, on line 444 the routine sets info=n+3 upon error in dtgsen.
However, then it checks whether the selected eigenvalues are only in
the upper part. If dtgsen has failed, it is almost certain that it is
the case that the selected eigenvalues are not only in the upper
part, and info=n+2 is set. In this way user does not get info=n+3,
but only info=n+2.
Many thanks for your great job.
Ondra K.
_______________________________________________
Lapack mailing list
Lapack@Domain.Removed <mailto:Lapack@Domain.Removed>
http://lists.eecs.utk.edu/mailman/listinfo/lapack
**********************************************
Julie Langou; Research Associate in Computer Science
Innovative Computing Laboratory;
University of Tennessee from Denver, Colorado ;-)
julie@Domain.Removed
<mailto:julie@Domain.Removed>; http://www.cs.utk.edu/~julie/
<http://www.cs.utk.edu/%7Ejulie/>
------------------------------------------------------------------------
_______________________________________________
Lapack mailing list
Lapack@Domain.Removed
http://lists.eecs.utk.edu/mailman/listinfo/lapack