Dear openfliiper,&nbsp;<div>Concerning the problem we sent last night.</div><div>We put two examples here. From the output,</div><div>we can see that for boundary vertex, the current way returns an unnecessary face(shown in bold font).</div>
<div><br></div><div>Besides, could you please help us how to iterate all incident vertices in a face in counterclockwise.</div><div><br></div><div><br></div><div>&gt;&gt;&gt;&gt;&gt;&gt;&gt; example one: 553 is a boundary vertex</div>
<div><br></div><div><div>-------------------------- retrieve by incident half edges</div><div>(553,517);(553,482)</div><div>(553,482);(553,516)</div><div>(553,516);(553,464)</div><div>(553,464);(553,603)</div><div><b>(553,603);(553,517)</b></div>
<div><br></div><div><div>--------------------------retrieve by incident faces</div><div>(553,517);(553,482)</div><div>(553,482);(553,516)</div><div>(553,516);(553,464)</div><div>(553,603);(553,464)</div></div><div><br></div>
<div><br></div><div>&gt;&gt;&gt;&gt;&gt;&gt;&gt; example two: 509 is an interior vertex
</div><div><br></div><div><div>-------------------------- retrieve by incident half edges
</div><div><div>(509,507);(509,534)</div><div>(509,534);(509,549)</div><div>(509,549);(509,535)</div><div>(509,535);(509,510)</div><div>(509,510);(509,492)</div><div>(509,492);(509,491)</div><div>(509,491);(509,507)</div>
</div></div><div><br></div><div>--------------------------retrieve by incident faces
</div><div><div>(509,507);(509,491)</div><div>(509,534);(509,507)</div><div>(509,549);(509,534)</div><div>(509,535);(509,549)</div><div>(509,535);(509,510)</div><div>(509,492);(509,510)</div><div>(509,492);(509,491)</div>
</div><div><br></div><div>Best Regards,</div><div>Xianyong Liu.</div><div><br></div><br><div class="gmail_quote">2012/10/24 ÁõÏÈÓÂ <span dir="ltr">&lt;<a href="mailto:xyliu80@gmail.com" target="_blank">xyliu80@gmail.com</a>&gt;</span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">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>---------------&nbsp;Curvature.cc --------------</div><div><div>for(; voh_it; ++voh_it, ++n_voh_it)</div>

<div>&nbsp; {</div><div>&nbsp; &nbsp; typename MeshT::Point p1 = _mesh.point(_mesh.to_vertex_handle( &nbsp; voh_it));</div><div>&nbsp; &nbsp; typename MeshT::Point p2 = _mesh.point(_mesh.to_vertex_handle( n_voh_it));</div><div><br></div><div>&nbsp; &nbsp; gauss_curv -= acos(OpenMesh::sane_aarg( ((p1-p0).normalize() | (p2-p0).normalize()) ));</div>

<div>&nbsp; }</div></div><div><br></div><div>Best Regards,</div><div>Xianyong Liu.</div><div><br></div>
</blockquote></div><br></div>