LAPACK Archives

[Lapack] Clapack dgesv accuracy

Dear Sumant,

I suspect that (some of) your matrices are ill-conditioned, which would 
explain why sometimes dposv thinks your matrices are not positive 
definite.  If this is the case, then you cannot expect different solvers 
to produce the same results.

Perhaps you should think about how you are modelling your application 
and see if you can produce better conditioned matrices?

Best wishes,

Sven Hammarling.

<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