CGAL users discussion list

[Andreas Fabri <>]

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