CGAL users discussion list

["Mukherjee, Soham" <>]

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:

.

Output to vtu: