Return more information if matrix decomposition fails

Submitted by Christoph Hertzberg @chhtz

Assigned to Nobody

Link to original bugzilla bug (#407)
Version: 3.1

Description

Created attachment 249
Check for non-positve definiteness

Currently SimplicialLLt results in NaN results for non-positive matrices, i.e. it takes the square root of a negative number.

Simple patch which checks for negative values and gracefully exits is attached.

Attachment 249, "Check for non-positve definiteness":
LLtPositiveDefinite.patch

Blocking

#631 (closed)

Edited by Christoph Hertzberg