CGAL users discussion list

[Andreas fabri <>]

The finite facets are all facets of the triangulation of the convex hull of the input. That is you only do not get facets with a vertex at infinity. 

Andreas

On 16 Aug 2019, at 06:22, Mukherjee, Soham <> wrote:

Hi All,

Is there any way to extract all the facets (not just the boundary facets) that belong to tetrahedra of C3T3::Triangulation? What I am doing is following:

1. Iterate over finite faces of C3T3::Triangulation object obtained from the function call refine_mesh(…..) 

C3t3 c3t3 = CGAL::make_mesh_3<C3t3>(domain, criteria);

C3t3::Triangulation tr = c3t3.triangulation();

for(auto fit=tr.finite_facets_begin();fit!=tr.finite_facets_end();fit++){

C3t3::Triangulation::Cell_handle c = fit->first;

int i = fit->second;

std::vector<unsigned> face_id;

face_id.push_back(vid[c->vertex((i+1)%4)]);

face_id.push_back(vid[c->vertex((i+2)%4)]);

face_id.push_back(vid[c->vertex((i+3)%4)]);

// Write this vertex indices of face i.e. connectivity information to file

}

But when I output the whole c3t3 to unstructuredgridVtk by calling output_to_vtu(c3t3) it is nicer (the way Delaunay Refined 3D Mesh should look like) and not like the mesh I obtained from manually outputting the info.

Am I missing something? That is needed to be taken care of?For sake of clarity I am attaching two pictures with this email.

Manual Output:

.<Manual_Output.png>

Output to vtu:

<Output_to_vtu.png>