Dear openfliiper,<div><br></div><div>We think, In Curvature.cc, the calculation for Gaussian curvature is wrong.</div><div><br></div><div>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.</div>
<div><br></div><div>Below is the existing calculation.</div><div>It should judge whether the vertex is a boundary one or not!</div><div><br></div><div>--------------- Curvature.cc --------------</div><div><div>for(; voh_it; ++voh_it, ++n_voh_it)</div>
<div>  {</div><div>    typename MeshT::Point p1 = _mesh.point(_mesh.to_vertex_handle(   voh_it));</div><div>    typename MeshT::Point p2 = _mesh.point(_mesh.to_vertex_handle( n_voh_it));</div><div><br></div><div>    gauss_curv -= acos(OpenMesh::sane_aarg( ((p1-p0).normalize() | (p2-p0).normalize()) ));</div>
<div>  }</div></div><div><br></div><div>Best Regards,</div><div>Xianyong Liu.</div><div><br></div>