From b8fc1edb2ccde053ba69110475ef15802b410fe5 Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Sun, 5 Oct 2008 13:45:43 +0000
Subject: [PATCH] Sparse module: enable support for incomplete cholesky
 factorization in CHOLMOD backend.

---
 Eigen/src/Sparse/CholmodSupport.h | 14 +++++++-------
 1 file changed, 7 insertions(+), 7 deletions(-)

diff --git a/Eigen/src/Sparse/CholmodSupport.h b/Eigen/src/Sparse/CholmodSupport.h
index 9a13ae883..0a57badc7 100644
--- a/Eigen/src/Sparse/CholmodSupport.h
+++ b/Eigen/src/Sparse/CholmodSupport.h
@@ -104,14 +104,14 @@ void SparseCholesky<MatrixType>::computeUsingCholmod(const MatrixType& a)
   cholmod_common c;
   cholmod_start(&c);
   cholmod_sparse A = const_cast<MatrixType&>(a).asCholmodMatrix();
-  std::vector<int> perm(a.cols());
-  for (int i=0; i<a.cols(); ++i)
-    perm[i] = i;
-  c.nmethods = 1;
-  c.method [0].ordering = CHOLMOD_NATURAL;
-  c.postorder = 0;
+  if (!(m_flags&CholPartial))
+  {
+    c.nmethods = 1;
+    c.method [0].ordering = CHOLMOD_NATURAL;
+    c.postorder = 0;
+  }
   c.final_ll = 1;
-  cholmod_factor *L = cholmod_analyze_p(&A, &perm[0], &perm[0], a.cols(), &c);
+  cholmod_factor *L = cholmod_analyze(&A, &c);
   cholmod_factorize(&A, L, &c);
   cholmod_sparse* cmRes = cholmod_factor_to_sparse(L, &c);
   m_matrix = CholMatrixType::Map(*cmRes);