Skip to Content.
Sympa Menu

cgal-discuss - [cgal-discuss] cutting the peanut in half: thin elements

Subject: CGAL users discussion list

List archive

[cgal-discuss] cutting the peanut in half: thin elements


Chronological Thread 
  • From: Michel Audette <>
  • To:
  • Subject: [cgal-discuss] cutting the peanut in half: thin elements
  • Date: Mon, 20 Jun 2011 14:58:54 -0400

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




Archive powered by MHonArc 2.6.16.

Top of Page