diff --git a/include/Algo/Geometry/curvature.h b/include/Algo/Geometry/curvature.h
index 817bbb6ee77eef70acded326a393b6bd0ce4eec9..5ec978d27776fc50c0302a8c58807d1b9e4c6b46 100644
--- a/include/Algo/Geometry/curvature.h
+++ b/include/Algo/Geometry/curvature.h
@@ -79,7 +79,7 @@ void vertexQuadraticFitting(
 	typename PFP::MATRIX33& localFrame,
 	const VertexAttribute<typename PFP::VEC3, typename PFP::MAP>& position,
 	const VertexAttribute<typename PFP::VEC3, typename PFP::MAP>& normal,
-	float& a, float& b, float& c, float& d, float& e) ;
+    typename PFP::REAL& a, typename PFP::REAL& b, typename PFP::REAL& c, typename PFP::REAL& d, typename PFP::REAL& e) ;
 
 template <typename PFP>
 void quadraticFittingAddVertexPos(
diff --git a/include/Algo/Geometry/curvature.hpp b/include/Algo/Geometry/curvature.hpp
index 00a606e7d372284b9e8fd9aadfe52aa0b4a109bb..60bad980f9c245ec0506e4142ac2348e3db3cc5d 100644
--- a/include/Algo/Geometry/curvature.hpp
+++ b/include/Algo/Geometry/curvature.hpp
@@ -136,7 +136,7 @@ void vertexQuadraticFitting(
 	typename PFP::MATRIX33& localFrame,
 	const VertexAttribute<typename PFP::VEC3, typename PFP::MAP>& position,
 	const VertexAttribute<typename PFP::VEC3, typename PFP::MAP>& normal,
-	float& a, float& b, float& c, float& d, float& e)
+    typename PFP::REAL& a, typename PFP::REAL& b, typename PFP::REAL& c, typename PFP::REAL& d, typename PFP::REAL& e)
 {
 	typename PFP::VEC3 p = position[v] ;
 
@@ -148,7 +148,7 @@ void vertexQuadraticFitting(
 	nlBegin(NL_MATRIX) ;
 	foreach_adjacent2<EDGE>(map, v, [&] (Vertex it) {
 		typename PFP::VEC3 itp = position[it] ;
-//		quadraticFittingAddVertexPos<PFP>(v, itp, localFrame) ;
+        quadraticFittingAddVertexPos<PFP>(p, itp, localFrame) ;
 		typename PFP::VEC3 itn = normal[it] ;
 		quadraticFittingAddVertexNormal<PFP>(itp, itn, p, localFrame) ;
 	});
@@ -629,7 +629,7 @@ void normalCycles_SortTensor(Geom::Matrix<3,3,typename PFP::REAL>& tensor, unsig
 }
 
 template <typename PFP>
-void normalCycles_ProjectTensor(Geom::Matrix<3,3,typename PFP::REAL>& tensor, const typename PFP::VEC3& normal_vector, unsigned int thread)
+void normalCycles_ProjectTensor(Geom::Matrix<3,3,typename PFP::REAL>& tensor, const typename PFP::VEC3& normal_vector, unsigned int /*thread*/)
 {
 	Geom::Matrix<3,3,typename PFP::REAL> proj;
 	proj.identity();