CGAL users discussion list

["Laurent Rineau (CGAL/GeometryFactory)" <>]

Le jeudi 17 octobre 2013 17:57:50 Garth N. Wells a écrit :
> I'm trying to create a volume mesh of a polyhedron with many faces, but
> I keep getting the error:
> 
> construct initial points:
> 8/8 initial point(s) found...
> Start surface scan...end scan. [Bad facets:12]
> 
> Refining Surface...
> Legende of the following line: (#vertices,#steps,#facets to refine,#tets
> to refine)
> (3286,3277,0,0) (2847.9 vertices/s))))
> Total refining surface time: 1.15069s
> 
> Start volume scan...end scan. [Bad tets:9713]
> 
> Refining...
> Legende of the following line: (#vertices,#steps,#facets to refine,#tets
> to refine)
> (7617,4330,0,37339) (11173.1 vertices/s)terminate called after throwing
> an instance of 'CGAL::Assertion_exception'
>    what():  CGAL ERROR: assertion violation!
> Expr: !r_tr_.is_infinite(*cell_it)
> File: /usr/include/CGAL/Mesh_3/Refine_cells_3.h
> Line: 385
> 
> 
> What I'm doing is surface meshing a sphere, converting the surface mesh
> to a polyhedron, and trying to generate a volume mesh of the polyhedron
> (the motivation is to add/subtract polyhedra to create more complicated
> domains). The code that reproduces the error (using CGAL 4.2 and GCC
> 4.8) is attached. The code is basically a concatenation of the demos in
> sections 48.3.1 and 51.3.2. Any suggestions to get this working?

The polyhedron generated by the intermediate step is probably not properly 
closed. Can you check?

-- 
Laurent Rineau, PhD
R&D Engineer at GeometryFactory           http://www.geometryfactory.com/
Release Manager of the CGAL Project       http://www.cgal.org/