CGAL users discussion list

Thank you for your prompt reply!

Am I wrong if I suppose that you use: c2t3.triangulation().dual(f), where
f is a facet such as c2t3.is_in_complex(f) == true ?

Yes, this is the requirement I set for facets inside the volume and I
obtain the expected 3D-behavior, as you stress.
Conversely, when I select surface triangles (if
((*it_facet).first->is_facet_on_surface((*it_facet).second)==true)) I
don't obtain points, but only rays (also "choosing" finite facets).


I think I can't help more since I don't know the behaviour you expect. Are
you trying to retrieve the circumcenter of the surface triangles ?

This is the first quantity I need, are you right. Actually I try to
evaluate the Voronoi surface area around each vertex generate by the
Surface mesher. So I select finite incident facets to each vertex, I check
if the triangles are on the surface  and I search the circumcenters.
Best