CGAL users discussion list

[Michel Audette <>]

Dear CGAL developers,

In the past week or so, I've become interested in the new capacity of 3.8 to offer boundary-faithful and controlled-resolution tetrahedralization.

Returning to the example in


examples/Mesh_3/mesh_two_implicit_spheres_with_balls.cpp,
playing with the mesh size a little bit, 

  Mesh_criteria criteria(edge_size = 0.5, 
                         facet_angle = 25, facet_size = 0.5, 

                         cell_radius_edge_ratio = 2, cell_size = 0.5); 


  c3t3 = CGAL::make_mesh_3<C3t3>(domain, criteria, 
                                      odt(time_limit=30), 
                                      perturb(sliver_bound=10, time_limit=30), 

                                      exude(time_limit=30, sliver_bound=15));

and visualizing the interior with Paraview (by producing a vtk file), 
  vtkSmartPointer<vtkUnstructuredGrid>

    outputUnstructuredGrid = vtkSmartPointer<vtkUnstructuredGrid>::New(); 

  CGAL::output_c3t3_to_vtk_unstructured_grid(c3t3, outputUnstructuredGrid);
  outputUnstructuredGrid->Squeeze();
  
  vtkSmartPointer<vtkUnstructuredGridWriter> writer =

    vtkSmartPointer<vtkUnstructuredGridWriter>::New();

  writer->SetInput(outputUnstructuredGrid);
  writer->SetFileName("out.vtk");
  writer->Write();

I'm finding nonetheless that some slivers remain in the peanut... See the enclosed visualization. I'd be happy to help track down the problem. Is this a parameter of odt(), perturb() or exude() that I am not using properly? 

Best wishes, 

Michel

 
--
Michel Audette, Ph.D.
R & D Engineer,
Kitware Inc.,
Chapel Hill, N.C.

Attachment: CGAL3.8peanutSlivers.gif
Description: GIF image