
24 Oct
2012
24 Oct
'12
4:56 p.m.
Dear openfliiper, We think, In Curvature.cc, the calculation for Gaussian curvature is wrong. The fact is that the number of incident angles(facets) overs a boundary vertex is one less than then number of incident vertices (edges) over it. Below is the existing calculation. It should judge whether the vertex is a boundary one or not! --------------- Curvature.cc -------------- for(; voh_it; ++voh_it, ++n_voh_it) { typename MeshT::Point p1 = _mesh.point(_mesh.to_vertex_handle( voh_it)); typename MeshT::Point p2 = _mesh.point(_mesh.to_vertex_handle( n_voh_it)); gauss_curv -= acos(OpenMesh::sane_aarg( ((p1-p0).normalize() | (p2-p0).normalize()) )); } Best Regards, Xianyong Liu.