Skip to Content.
Sympa Menu

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

Subject: CGAL users discussion list

List archive

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


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

Hi Andreas,

thanks for your prompt reply, and I apologize if the email just sent is also covered by this workaround.

I'll give it a try and report back.

Best wishes,

Michel

On Mon, Jun 20, 2011 at 3:35 PM, Andreas Fabri <> wrote:

Hi Michel,

Some days ago we found a bug in the code that removes slivers.
The workaround is to define the macro
CGAL_MESH_3_USE_RELAXED_HEAP

Can you report back if it helped in your case.

andreas



On 20/06/2011 20:58, Michel Audette wrote:
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.



--
Andreas Fabri, PhD
Chief Officer, GeometryFactory
Editor, The CGAL Project

phone: +33.492.954.912    skype: andreas.fabri

--
You are currently subscribed to cgal-discuss.
To unsubscribe or access the archives, go to
https://lists-sop.inria.fr/wws/info/cgal-discuss




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




Archive powered by MHonArc 2.6.16.

Top of Page