Menu

#25 LaQRLinearSolveIP does not recognize info value > 0

open
nobody
7
2012-09-14
2012-01-11
No

Hi,

I am referring to Lapack++ 2.5.4:
In LaQRLinearSolveIP is a call to the Fortran function dgels. According to [1] the passed info-value can be 0 (which does not have to be handled), +1 or -1.
The info-value is checked in the following way:
if (info < 0)
throw(LaException(fname, "Internal error in LAPACK: SGELS()"));

On one hand the exception string is wrong (it should at least mean "Internal error in LAPACK: DGELS()") and on the other hand, info > 0 is not checked.
If I am right, it is not possible for the user to check if LaQRLinearSolveIP returned a correct solution of the system.

[1] http://www.netlib.org/lapack/explore-html/d7/d3b/group__double_g_esolve.html

Discussion

Anonymous
Anonymous

Add attachments
Cancel





Want the latest updates on software, tech news, and AI?
Get latest updates about software, tech news, and AI from SourceForge directly in your inbox once a month.