songsenand
/
cgal
mirror of https://github.com/CGAL/cgal
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

avoid insertion of points inside protecting balls, during initialization step

This commit is contained in:
Jane Tournois 2024-03-14 16:40:28 +01:00
parent 4e81649dd9
commit 80fd8759dd
1 changed files with 10 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
Mesh_3/include/CGAL/make_mesh_3.h
View File

@ -45,9 +45,11 @@ init_c3t3(C3T3& c3t3, const MeshDomain& domain, const MeshCriteria&,
{ {
typedef typename MeshDomain::Point_3 Point_3; typedef typename MeshDomain::Point_3 Point_3;
typedef typename MeshDomain::Index Index; typedef typename MeshDomain::Index Index;
typedef std::vector<std::pair<Point_3, Index> > Initial_points_vector; typedef typename std::pair<Point_3, Index> PI;
typedef std::vector<PI> Initial_points_vector;
typedef typename Initial_points_vector::iterator Ipv_iterator; typedef typename Initial_points_vector::iterator Ipv_iterator;
typedef typename C3T3::Vertex_handle Vertex_handle; typedef typename C3T3::Vertex_handle Vertex_handle;
typedef CGAL::Mesh_3::Triangulation_helpers<typename C3T3::Triangulation> Th;
// Mesh initialization : get some points and add them to the mesh // Mesh initialization : get some points and add them to the mesh
Initial_points_vector initial_points; Initial_points_vector initial_points;
@ -61,17 +63,19 @@ init_c3t3(C3T3& c3t3, const MeshDomain& domain, const MeshCriteria&,
c3t3.triangulation().geom_traits().construct_weighted_point_3_object(); c3t3.triangulation().geom_traits().construct_weighted_point_3_object();
// Insert points and set their index and dimension // Insert points and set their index and dimension
for ( Ipv_iterator it = initial_points.begin() ; const Vertex_handle vhint = c3t3.triangulation().finite_vertices_begin();
it != initial_points.end() ; for (const PI& pi : initial_points)
++it )
{ {
Vertex_handle v = c3t3.triangulation().insert(cwp(it->first)); if(Th().inside_protecting_balls(c3t3.triangulation(), vhint, pi.first))
continue;
Vertex_handle v = c3t3.triangulation().insert(cwp(pi.first));
// v could be null if point is hidden // v could be null if point is hidden
if ( v != Vertex_handle() ) if ( v != Vertex_handle() )
{ {
c3t3.set_dimension(v,2); // by construction, points are on surface c3t3.set_dimension(v,2); // by construction, points are on surface
c3t3.set_index(v,it->second); c3t3.set_index(v, pi.second);
} }
} }
} }