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

first solution to correctly re-constrain edges refined

This commit is contained in:
Sébastien Loriot 2018-04-29 13:13:58 +02:00
parent 95d76709e8
commit ab9d668f52
1 changed files with 12 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
Polygon_mesh_processing/include/CGAL/Polygon_mesh_processing/internal/Isotropic_remeshing/remesh_impl.h
View File

@ -287,6 +287,7 @@ namespace internal {
typedef typename boost::property_map<
PM, CGAL::dynamic_halfedge_property_t<Halfedge_status> >::type Halfedge_status_pmap;
typedef boost::unordered_set<edge_descriptor> Edge_is_constrained_set;
public:
Incremental_remesher(PolygonMesh& pmesh
@ -314,6 +315,13 @@ namespace internal {
BOOST_FOREACH(halfedge_descriptor h, halfedges(mesh_))
put(halfedge_status_pmap_, h, MESH);
/// \todo the following is not working in case the default map is passed (we need to have a way to detect it)
BOOST_FOREACH(edge_descriptor e, edges(pmesh))
{
if ( get(ecmap_, e) )
edge_is_constrained_set_.insert(e);
}
CGAL_assertion(CGAL::is_triangle_mesh(mesh_));
}
@ -413,6 +421,7 @@ namespace internal {
Point refinement_point = this->midpoint(he);
halfedge_descriptor hnew = CGAL::Euler::split_edge(he, mesh_);
CGAL_assertion(he == next(hnew, mesh_));
if (edge_is_constrained_set_.count(edge(he, mesh_))!=0) edge_is_constrained_set_.insert(edge(hnew, mesh_));
++nb_splits;
//move refinement point
@ -511,6 +520,7 @@ namespace internal {
Point refinement_point = this->midpoint(he);
halfedge_descriptor hnew = CGAL::Euler::split_edge(he, mesh_);
CGAL_assertion(he == next(hnew, mesh_));
if (edge_is_constrained_set_.count(edge(he, mesh_))!=0) edge_is_constrained_set_.insert(edge(hnew, mesh_));
++nb_splits;
//move refinement point
@ -1911,8 +1921,7 @@ private:
{
BOOST_FOREACH(edge_descriptor e, edges(mesh_))
{
if (is_on_patch_border(halfedge(e, mesh_))
|| is_on_patch_border(opposite(halfedge(e, mesh_), mesh_)))
if (edge_is_constrained_set_.count(e))
put(ecmap_, e, true);
else
put(ecmap_, e, false);
@ -1929,6 +1938,7 @@ private:
Triangle_list input_triangles_;
Patch_id_list input_patch_ids_;
Halfedge_status_pmap halfedge_status_pmap_;
Edge_is_constrained_set edge_is_constrained_set_;
bool protect_constraints_;
FacePatchMap patch_ids_map_;
EdgeIsConstrainedMap ecmap_;