CGAL users discussion list

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

Le Monday 13 January 2014 02:42:05 Marek a écrit :
> Hi
> I've implemented it as you said some time ago and everything looked alright
> but now i've stumbled upon situation in which some points are neglected in
> final mesh.
> 
> I'm adding them like this:
> 
> C3t3 c3t3;
> CGAL::internal::Mesh_3::init_c3t3_with_features(c3t3, domain, criteria);
> Tr& tr = c3t3.triangulation();
> 
> Then, for each point I want to preserve:
>           Vertex_handle vh = tr.insert(v->point());
>           c3t3.add_to_complex(vh,1);
> 
> where v are of type Polyhedron::Vertex_handle
> 
> It looks like at least some of inserted points are stored as verices in
> triangulation but are not part of any facet/cell. Instead new points are
> generated fairly close to ones I want to preserve.
> Is there something wrong with my method of insertion to c3t3? How can I fix
> it?

Indeed, there is no guaranty that there are facets or cells in the c3t3 
object 
that will be incident to the points you have inserted that way. Maybe your 
domain, and its featured curves, is not coherent with the points you have 
added.

Please explain what you really want to do. Probably the forced insertion of 
vertices in the c3t3 is not really want you need.

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