I am using MTL (a c++ interface to lapack), so I am not sure this is a Lapack or a MTL problem. In any case I ask this user group hoping in some hints.
Assume I try to solve the linear system Ax=b with
obviously I get an error (A is singular and with this b there is no solution). In fact mtl method getsv, which interfaces to lapack, returns:
> Factorization failed with INFO = 2
However if I change A to
I get a solution!:
Thank you for any help,
This is clearly an error. Did anybody noticed it before?