CGAL users discussion list

["Laurent Rineau (CGAL/GeometryFactory)" <>]

Le mardi 15 octobre 2013 23:35:28 Burak ER a écrit :

> Your solution seems to get all the facets either they are interface or

> a surface facets. How can I determine the surface facets? A zero

> neighbor value indicates it as a surface facet?

 

If your surface mesh is in a C2t3 object (from the Surface Mesh Generator package), change the loop to:

 

for(C2t3::Facets_iterator

fit = c2t3.facets_begin(),

end = c2t3.facets_end(); fit !=end; ++fit)

 

If it is in a C3t3 object (from the 3D Mesh Generator package), change the loop to:

 

for(C3t3::Facets_in_complex_iterator

fit = c2t3.facets_in_complex_begin(),

end = c2t3.facets_in_complex_end(); fit !=end; ++fit)

In all cases, the value type of the iterators is Tr::Facet, that is a typedef for `std::pair<Tr::Cell_handle, int>`.

 

--

Laurent Rineau, PhD

R&D Engineer at GeometryFactory http://www.geometryfactory.com/

Release Manager of the CGAL Project http://www.cgal.org/