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.
openflipper@lists.rwth-aachen.de