From 565f4c4f53f80105e4bddf967c033d1365d276a0 Mon Sep 17 00:00:00 2001
From: nico <nico.pietroni@gmail.com>
Date: Mon, 6 May 2019 18:17:17 +1000
Subject: [PATCH] added k-neighbors instead of radius to make it more stable

---
 wrap/igl/smooth_field.h | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/wrap/igl/smooth_field.h b/wrap/igl/smooth_field.h
index 34e514ca..e507ceb9 100644
--- a/wrap/igl/smooth_field.h
+++ b/wrap/igl/smooth_field.h
@@ -406,7 +406,8 @@ public:
 
         Eigen::MatrixXd PD1,PD2,PV1,PV2;
         MeshToMatrix<MeshType>::GetTriMeshData(mesh,F,V);
-        igl::principal_curvature(V,F,PD1,PD2,PV1,PV2,Nring);
+
+        igl::principal_curvature(V,F,PD1,PD2,PV1,PV2,Nring,true);
 
         //then copy curvature per vertex
         for (size_t i=0;i<mesh.vert.size();i++)
@@ -497,7 +498,6 @@ public:
 //             (SParam.curv_thr>0))
 
         InitByCurvature(mesh,SParam.curvRing);
-
         SelectConstraints(mesh,SParam);
         //then do the actual smooth
         SmoothDirections(mesh,SParam.Ndir,SParam.SmoothM,true,SParam.alpha_curv);