Detect factorisation issues in sparse cholesky

pkhuong · Dec 19, 2013 · de30e45 · de30e45
1 parent f625715
commit de30e45
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/sparse-cholesky.lisp b/sparse-cholesky.lisp
@@ -382,7 +382,10 @@
                (b (* cholmod-dense) :local (make-dense-from-matlisp b))
                (factor (* cholmod-factor) :local (cholmod-analyze As
                                                                   common)))
+    (cholmod-set-status common 0)
     (cholmod-factorize As factor common)
+    (when (/= (cholmod-get-status common) 0)
+      (return-from solve-dense))
     (with-alien ((x (* cholmod-dense) :local (cholmod-solve 0
                                                             factor
                                                             b