@@ -840,9 +840,6 @@ Field Mesh::curvature(const CurvatureType type) const {
840840 case Principal: {
841841 curv->SetName (" principal curvature" );
842842
843- // returns maximal curvature value for each vertex
844- // igl::principal_curvature(V, F, PD1, PD2, PV1, PV2);
845-
846843 // returns minimal curvature value for each vertex
847844 igl::principal_curvature (V, F, PD1, PD2, C, PV2);
848845 break ;
@@ -855,18 +852,6 @@ Field Mesh::curvature(const CurvatureType type) const {
855852 case Mean: {
856853 curv->SetName (" mean curvature" );
857854
858- Eigen::MatrixXd HN;
859- Eigen::SparseMatrix<double > L, M, Minv;
860- igl::cotmatrix (V, F, L);
861- igl::massmatrix (V, F, igl::MASSMATRIX_TYPE_VORONOI, M);
862- igl::invert_diag (M, Minv);
863-
864- // Laplace-Beltrami of position
865- HN = -Minv * (L * V);
866-
867- // Extract magnitude as mean curvature
868- C = HN.rowwise ().norm ();
869-
870855 // Compute curvature directions via quadric fitting
871856 igl::principal_curvature (V, F, PD1, PD2, PV1, PV2);
872857
0 commit comments