From 4e7ec01ed2d33ae755bd780d49ef90d831719983 Mon Sep 17 00:00:00 2001
From: Daniele Panozzo <daniele.panozzo@gmail.com>
Date: Tue, 25 Oct 2016 21:40:20 -0400
Subject: [PATCH] improves performance of 3D slim by ~2x

---
 include/igl/slim.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/include/igl/slim.cpp b/include/igl/slim.cpp
index 437f9c9fd..a8d7e7387 100644
--- a/include/igl/slim.cpp
+++ b/include/igl/slim.cpp
@@ -361,7 +361,7 @@ void WeightedGlobalLocal::solve_weighted_arap(const Eigen::MatrixXd& V, const Ei
   } else { // seems like CG performs much worse for 2D and way better for 3D
     Eigen::VectorXd guess(uv.rows()*dim);
     for (int i = 0; i < dim; i++) for (int j = 0; j < dim; j++) guess(uv.rows()*i + j) = uv(i,j); // flatten vector
-    ConjugateGradient<SparseMatrix<double>, Eigen::Upper> solver;
+    ConjugateGradient<SparseMatrix<double>, Eigen::Lower|Upper> solver;
     solver.setTolerance(1e-8);
     Uc = solver.compute(L).solveWithGuess(rhs,guess);
   }