[Lapack] INFO=N+2 in xDGGES

[James Demmel]

Thanks also for your bug report!
Do you have input data that caused this bug to occur,
that you could share with us to try to debug it?
Or did you just notice it by reading the code?
Or did you discover it another way?
Thanks,
Jim Demmel

julie langou wrote:
> 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
>