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

Remove CGAL_NEF3_TRIANGULATE_FACETS code

This commit is contained in:
Giles Bathgate 2021-05-31 17:05:45 +01:00
parent 5e9f9d313f
commit f30eeedda4
3 changed files with 0 additions and 576 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

217
Nef_3/include/CGAL/Nef_3/K3_tree.h
View File

@ -24,14 +24,6 @@
#include <CGAL/Cartesian.h> #include <CGAL/Cartesian.h>
#include <boost/container/deque.hpp> #include <boost/container/deque.hpp>
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
#include <CGAL/Constrained_triangulation_2.h>
#include <CGAL/Triangulation_data_structure_2.h>
#include <CGAL/Projection_traits_xy_3.h>
#include <CGAL/Projection_traits_yz_3.h>
#include <CGAL/Projection_traits_xz_3.h>
#include <CGAL/Constrained_triangulation_face_base_2.h>
#endif
#include <deque> #include <deque>
#include <sstream> #include <sstream>
@ -44,43 +36,6 @@
namespace CGAL { namespace CGAL {
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
template <typename Triangle_3>
void sort_triangle_by_lexicographically_smaller_vertex
(const Triangle_3& t, int& a, int& b, int& c) {
typedef typename Triangle_3::R Kernel;
a = 0;
for( int i = 1; i < 3; ++i) {
if( compare_xyz<Kernel>( t[a], t[i]) == SMALLER)
a = i;
}
b = (a + 1) % 3;
c = (b + 1) % 3;
if( compare_xyz<Kernel>( t[b], t[c]) == LARGER)
std::swap( b, c);
return;
}
template <typename Triangle_3>
struct Compare_triangle_3 {
typedef typename Triangle_3::R Kernel;
bool operator()( const Triangle_3& t1, const Triangle_3& t2) const {
int v1[3], v2[3];
sort_triangle_by_lexicographically_smaller_vertex
( t1, v1[0], v1[1], v1[2]);
sort_triangle_by_lexicographically_smaller_vertex
( t2, v2[0], v2[1], v2[2]);
for( int i = 0; i < 3; ++i) {
Comparison_result c = compare_xyz<Kernel>( t1[v1[i]], t2[v2[i]]);
if( c == SMALLER)
return true;
else if( c == LARGER)
return false;
}
return false; // the two triangles are equivalent
}
};
#endif
template <class Traits> template <class Traits>
class K3_tree class K3_tree
@ -161,89 +116,6 @@ private:
Coordinate coord; Coordinate coord;
}; };
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
template<typename SNC_structure, typename Kernel>
class Triangulation_handler {
typedef typename CGAL::Triangulation_vertex_base_2<Kernel> Vb;
typedef typename CGAL::Constrained_triangulation_face_base_2<Kernel> Fb;
typedef typename CGAL::Triangulation_data_structure_2<Vb,Fb> TDS;
typedef typename CGAL::No_constraint_intersection_requiring_constructions_tag Itag;
typedef typename CGAL::Constrained_triangulation_2<Kernel,TDS,Itag> CT;
typedef typename CT::Face_handle Face_handle;
typedef typename CT::Finite_faces_iterator Finite_face_iterator;
typedef typename CT::Edge Edge;
typedef typename SNC_structure::Halffacet_cycle_iterator
Halffacet_cycle_iterator;
typedef typename SNC_structure::SHalfedge_around_facet_circulator
SHalfedge_around_facet_circulator;
CT ct;
CGAL::Unique_hash_map<Face_handle, bool> visited;
Finite_face_iterator fi;
public:
template<typename Halffacet_handle>
Triangulation_handler(Halffacet_handle f) : visited(false) {
typedef typename SNC_structure::Halffacet_cycle_iterator
Halffacet_cycle_iterator;
typedef typename SNC_structure::SHalfedge_around_facet_circulator
SHalfedge_around_facet_circulator;
Halffacet_cycle_iterator fci;
for(fci=f->facet_cycles_begin(); fci!=f->facet_cycles_end(); ++fci) {
if(fci.is_shalfedge()) {
SHalfedge_around_facet_circulator sfc(fci), send(sfc);
CGAL_For_all(sfc,send) {
ct.insert_constraint(sfc->source()->source()->point(),
sfc->source()->twin()->source()->point());
}
}
}
CGAL_assertion(ct.is_valid());
typename CT::Face_handle infinite = ct.infinite_face();
typename CT::Vertex_handle ctv = infinite->vertex(1);
if(ct.is_infinite(ctv)) ctv = infinite->vertex(2);
CGAL_assertion(!ct.is_infinite(ctv));
typename CT::Face_handle opposite;
typename CT::Face_circulator vc(ctv,infinite);
do { opposite = vc++;
} while(!ct.is_constrained(CT::Edge(vc,vc->index(opposite))));
typename CT::Face_handle first = vc;
traverse_triangulation(first, first->index(opposite));
fi = ct.finite_faces_begin();
}
void traverse_triangulation(Face_handle f, int parent) {
visited[f] = true;
if(!ct.is_constrained(Edge(f,ct.cw(parent))) && !visited[f->neighbor(ct.cw(parent))]) {
Face_handle child(f->neighbor(ct.cw(parent)));
traverse_triangulation(child, child->index(f));
}
if(!ct.is_constrained(Edge(f,ct.cw(parent))) && !visited[f->neighbor(ct.cw(parent))]) {
Face_handle child(f->neighbor(ct.cw(parent)));
traverse_triangulation(child, child->index(f));
}
}
template<typename Triangle_3>
bool get_next_triangle(Triangle_3& tr) {
if(fi == ct.finite_faces_end()) return false;
++fi;
while(fi != ct.finite_faces_end() && visited[fi] == false) ++fi;
if(fi == ct.finite_faces_end()) return false;
tr = Triangle_3(fi->vertex(0)->point(), fi->vertex(1)->point(), fi->vertex(2)->point());
return true;
}
};
#endif
public: public:
friend class Objects_along_ray; friend class Objects_along_ray;
@ -257,9 +129,6 @@ typedef typename Traits::Infimaximal_box Infimaximal_box;
typedef typename Traits::Vertex_handle Vertex_handle; typedef typename Traits::Vertex_handle Vertex_handle;
typedef typename Traits::Halfedge_handle Halfedge_handle; typedef typename Traits::Halfedge_handle Halfedge_handle;
typedef typename Traits::Halffacet_handle Halffacet_handle; typedef typename Traits::Halffacet_handle Halffacet_handle;
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
typedef typename Traits::Halffacet_triangle_handle Halffacet_triangle_handle;
#endif
typedef typename Traits::Object_handle Object_handle; typedef typename Traits::Object_handle Object_handle;
typedef std::vector<Object_handle> Object_list; typedef std::vector<Object_handle> Object_list;
typedef typename Object_list::const_iterator Object_const_iterator; typedef typename Object_list::const_iterator Object_const_iterator;
@ -318,16 +187,6 @@ public:
left_node->transform(t); left_node->transform(t);
right_node->transform(t); right_node->transform(t);
splitting_plane = splitting_plane.transform(t); splitting_plane = splitting_plane.transform(t);
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
} else {
Halffacet_triangle_handle tri;
typename Object_list::iterator o;
for(o = object_list.begin(); o != object_list.end(); ++o)
if(assign(tri,*o)) {
tri.transform(t);
*o = make_object(tri);
}
#endif // CGAL_NEF3_TRIANGULATE_FACETS
} }
} }
@ -802,55 +661,7 @@ public:
CGAL_forall_edges( e, *W) CGAL_forall_edges( e, *W)
objects.push_back(make_object(Halfedge_handle(e))); objects.push_back(make_object(Halfedge_handle(e)));
CGAL_forall_facets( f, *W) { CGAL_forall_facets( f, *W) {
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
typedef typename SNC_structure::SHalfedge_around_facet_circulator
SHalfedge_around_facet_circulator;
typedef typename SNC_structure::Halffacet_cycle_iterator
Halffacet_cycle_iterator;
Halffacet_cycle_iterator fci = f->facet_cycles_begin();
CGAL_assertion(fci.is_shalfedge());
SHalfedge_around_facet_circulator safc(fci);
if(circulator_size(safc) > 10) {
typedef typename CGAL::Projection_traits_xy_3<Kernel> XY;
typedef typename CGAL::Projection_traits_yz_3<Kernel> YZ;
typedef typename CGAL::Projection_traits_xz_3<Kernel> XZ;
Triangle_3 tr;
Vector_3 orth = f->plane().orthogonal_vector();
int c = CGAL::abs(orth[0]) > CGAL::abs(orth[1]) ? 0 : 1;
c = CGAL::abs(orth[2]) > CGAL::abs(orth[c]) ? 2 : c;
std::list<Triangle_3> triangles;
if(c == 0) {
Triangulation_handler<SNC_structure, YZ> th(f);
while(th.get_next_triangle(tr)) {
triangles.push_front(tr);
Halffacet_triangle_handle th( f, *triangles.begin());
objects.push_back(make_object(th));
}
} else if(c == 1) {
Triangulation_handler<SNC_structure, XZ> th(f);
while(th.get_next_triangle(tr)) {
triangles.push_front(tr);
Halffacet_triangle_handle th( f, *triangles.begin());
objects.push_back(make_object(th));
}
} else if(c == 2) {
Triangulation_handler<SNC_structure, XY> th(f);
while(th.get_next_triangle(tr)) {
triangles.push_front(tr);
Halffacet_triangle_handle th( f, *triangles.begin());
objects.push_back(make_object(th));
}
} else
CGAL_error_msg( "wrong value");
} else
objects.push_back(make_object(Halffacet_handle(f)));
#else
objects.push_back(make_object(Halffacet_handle(f))); objects.push_back(make_object(Halffacet_handle(f)));
#endif // CGAL_NEF3_TRIANGULATE_FACETS
} }
Object_iterator oli=objects.begin()+v_end; Object_iterator oli=objects.begin()+v_end;
root = build_kdtree( objects, oli, 0); root = build_kdtree( objects, oli, 0);
@ -887,9 +698,6 @@ typename Object_list::difference_type n_vertices = std::distance(objects.begin()
Unique_hash_map< Vertex_handle, bool> v_mark(false); Unique_hash_map< Vertex_handle, bool> v_mark(false);
Unique_hash_map< Halfedge_handle, bool> e_mark(false); Unique_hash_map< Halfedge_handle, bool> e_mark(false);
Unique_hash_map< Halffacet_handle, bool> f_mark(false); Unique_hash_map< Halffacet_handle, bool> f_mark(false);
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
std::map< Triangle_3, bool, Compare_triangle_3<Triangle_3> > t_mark;
#endif
for( typename Objects_around_segment::Iterator oar = objects.begin(); for( typename Objects_around_segment::Iterator oar = objects.begin();
oar != objects.end(); ++oar) { oar != objects.end(); ++oar) {
for( typename Object_list::const_iterator o = (*oar).begin(); for( typename Object_list::const_iterator o = (*oar).begin();
@ -897,9 +705,6 @@ typename Object_list::difference_type n_vertices = std::distance(objects.begin()
Vertex_handle v; Vertex_handle v;
Halfedge_handle e; Halfedge_handle e;
Halffacet_handle f; Halffacet_handle f;
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
Halffacet_triangle_handle t;
#endif
if( CGAL::assign( v, *o)) { if( CGAL::assign( v, *o)) {
if( !v_mark[v]) { if( !v_mark[v]) {
O.push_back(*o); O.push_back(*o);
@ -918,15 +723,6 @@ typename Object_list::difference_type n_vertices = std::distance(objects.begin()
f_mark[f] = true; f_mark[f] = true;
} }
} }
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
else if(CGAL::assign(t, *o)) {
Triangle_3 tr = t.get_triangle();
if( !t_mark[tr]) {
O.push_back(*o);
t_mark[tr] = true;
}
}
#endif
else else
CGAL_error_msg( "wrong handle"); CGAL_error_msg( "wrong handle");
} }
@ -1013,10 +809,6 @@ std::string dump_object_list( const Object_list& O, int level = 0) {
Vertex_handle v; Vertex_handle v;
Halfedge_handle e; Halfedge_handle e;
Halffacet_handle f; Halffacet_handle f;
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
typename Object_list::size_type t_count = 0;
Halffacet_triangle_handle t;
#endif
for( o = O.begin(); o != O.end(); ++o) { for( o = O.begin(); o != O.end(); ++o) {
if( CGAL::assign( v, *o)) { if( CGAL::assign( v, *o)) {
if( level) os << v->point() << std::endl; if( level) os << v->point() << std::endl;
@ -1031,19 +823,10 @@ std::string dump_object_list( const Object_list& O, int level = 0) {
if( level) os << "facet" << std::endl; if( level) os << "facet" << std::endl;
++f_count; ++f_count;
} }
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
else if( CGAL::assign(t, *o)) {
if( level) os << "triangle" << std::endl;
++t_count;
}
#endif
else else
CGAL_error_msg( "wrong handle"); CGAL_error_msg( "wrong handle");
} }
os << v_count << "v " << e_count << "e " << f_count << "f "; os << v_count << "v " << e_count << "e " << f_count << "f ";
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
os << t_count << "t ";
#endif
return os.str(); return os.str();
} }

76
Nef_3/include/CGAL/Nef_3/SNC_k3_tree_traits.h
View File

@ -108,10 +108,6 @@ public:
typedef typename Decorator_traits::Halfedge_handle Halfedge_handle; typedef typename Decorator_traits::Halfedge_handle Halfedge_handle;
typedef typename Decorator_traits::Halffacet_handle Halffacet_handle; typedef typename Decorator_traits::Halffacet_handle Halffacet_handle;
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
typedef typename SNC_structure::Halffacet_triangle_handle
Halffacet_triangle_handle;
#endif
typedef typename SNC_structure::Object_handle Object_handle; typedef typename SNC_structure::Object_handle Object_handle;
typedef typename Decorator_traits::Halffacet_cycle_iterator typedef typename Decorator_traits::Halffacet_cycle_iterator
@ -145,10 +141,6 @@ public:
( const Point_3& pop, Halfedge_handle e, Depth depth); ( const Point_3& pop, Halfedge_handle e, Depth depth);
template<typename Depth> Oriented_side operator() template<typename Depth> Oriented_side operator()
( const Point_3& pop, Halffacet_handle f, Depth depth); ( const Point_3& pop, Halffacet_handle f, Depth depth);
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
template<typename Depth> Oriented_side operator()
( const Point_3& pop, Halffacet_triangle_handle f, Depth depth);
#endif
#ifdef CGAL_NEF_EXPLOIT_REFERENCE_COUNTING #ifdef CGAL_NEF_EXPLOIT_REFERENCE_COUNTING
bool reference_counted; bool reference_counted;
#endif #endif
@ -287,10 +279,6 @@ public:
typedef typename Decorator_traits::Vertex_handle Vertex_handle; typedef typename Decorator_traits::Vertex_handle Vertex_handle;
typedef typename Decorator_traits::Halfedge_handle Halfedge_handle; typedef typename Decorator_traits::Halfedge_handle Halfedge_handle;
typedef typename Decorator_traits::Halffacet_handle Halffacet_handle; typedef typename Decorator_traits::Halffacet_handle Halffacet_handle;
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
typedef typename SNC_structure::Halffacet_triangle_handle
Halffacet_triangle_handle;
#endif
typedef typename SNC_structure::Object_handle Object_handle; typedef typename SNC_structure::Object_handle Object_handle;
typedef std::vector<Object_handle> Object_list; typedef std::vector<Object_handle> Object_list;
@ -332,19 +320,12 @@ Side_of_plane<SNC_decorator>::operator()
Vertex_handle v; Vertex_handle v;
Halfedge_handle e; Halfedge_handle e;
Halffacet_handle f; Halffacet_handle f;
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
Halffacet_triangle_handle t;
#endif
if( CGAL::assign( v, o)) if( CGAL::assign( v, o))
return (*this)(pop, v, depth); return (*this)(pop, v, depth);
else if( CGAL::assign( e, o)) else if( CGAL::assign( e, o))
return (*this)(pop, e, depth); return (*this)(pop, e, depth);
else if( CGAL::assign( f, o)) else if( CGAL::assign( f, o))
return (*this)(pop, f, depth); return (*this)(pop, f, depth);
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
else if( CGAL::assign( t, o))
return (*this)(pop, t, depth);
#endif
else else
CGAL_error_msg( "wrong handle"); CGAL_error_msg( "wrong handle");
@ -461,63 +442,6 @@ Side_of_plane<SNC_decorator>::operator()
return ON_ORIENTED_BOUNDARY; return ON_ORIENTED_BOUNDARY;
} }
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
template <typename SNC_decorator>
template <typename Depth>
Oriented_side
Side_of_plane<SNC_decorator>::operator()
( const Point_3& pop, Halffacet_triangle_handle t, Depth depth) {
bool on_positive_side = false, on_negative_side = false;
Triangle_3 tr(t.get_triangle());
for( int i = 0; i < 3; ++i) {
Oriented_side side = ON_ORIENTED_BOUNDARY;
Comparison_result cr;
if(
#ifdef CGAL_NEF_EXPLOIT_REFERENCE_COUNTING
!reference_counted ||
#endif
!OnSideMapRC.is_defined(&(tr[i].hw()))) {
switch(depth%3) {
case 0:
cr = CGAL::compare_x(tr[i], pop);
side = cr == LARGER ? ON_POSITIVE_SIDE :
cr == SMALLER ? ON_NEGATIVE_SIDE : ON_ORIENTED_BOUNDARY;
break;
case 1:
cr = CGAL::compare_y(tr[i], pop);
side = cr == LARGER ? ON_POSITIVE_SIDE :
cr == SMALLER ? ON_NEGATIVE_SIDE : ON_ORIENTED_BOUNDARY;
break;
case 2:
cr = CGAL::compare_z(tr[i], pop);
side = cr == LARGER ? ON_POSITIVE_SIDE :
cr == SMALLER ? ON_NEGATIVE_SIDE : ON_ORIENTED_BOUNDARY;
break;
default: CGAL_error_msg( "wrong value");
}
#ifdef CGAL_NEF_EXPLOIT_REFERENCE_COUNTING
if(reference_counted)
OnSideMapRC[&(tr[i].hw())] = side;
} else if(reference_counted)
side = OnSideMapRC[&(tr[i].hw())];
#endif
if( side == ON_POSITIVE_SIDE)
on_positive_side = true;
else if( side == ON_NEGATIVE_SIDE)
on_negative_side = true;
}
if( on_positive_side && on_negative_side)
return ON_ORIENTED_BOUNDARY;
if( !on_positive_side && !on_negative_side)
return ON_ORIENTED_BOUNDARY;
if( on_positive_side) {
CGAL_assertion( !on_negative_side);
return ON_POSITIVE_SIDE;
}
CGAL_assertion( on_negative_side);
return ON_NEGATIVE_SIDE;
}
#endif
/* /*
As for the edges, if a facet is tanget to the plane it is not considered as As for the edges, if a facet is tanget to the plane it is not considered as

283
Nef_3/include/CGAL/Nef_3/SNC_point_locator.h
View File

@ -24,14 +24,6 @@
#include <CGAL/Timer.h> #include <CGAL/Timer.h>
#include <string> #include <string>
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
#include <CGAL/Constrained_triangulation_2.h>
#include <CGAL/Triangulation_data_structure_2.h>
#include <CGAL/Projection_traits_xy_3.h>
#include <CGAL/Projection_traits_yz_3.h>
#include <CGAL/Projection_traits_xz_3.h>
#include <CGAL/Constrained_triangulation_face_base_2.h>
#endif
#undef CGAL_NEF_DEBUG #undef CGAL_NEF_DEBUG
#define CGAL_NEF_DEBUG 509 #define CGAL_NEF_DEBUG 509
@ -164,10 +156,6 @@ public:
typedef K3_tree SNC_candidate_provider; typedef K3_tree SNC_candidate_provider;
typedef typename SNC_structure::Object_handle Object_handle; typedef typename SNC_structure::Object_handle Object_handle;
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
typedef typename Decorator_traits::Halffacet_triangle_handle
Halffacet_triangle_handle;
#endif
typedef typename SNC_structure::Point_3 Point_3; typedef typename SNC_structure::Point_3 Point_3;
typedef typename SNC_structure::Plane_3 Plane_3; typedef typename SNC_structure::Plane_3 Plane_3;
typedef typename SNC_structure::Segment_3 Segment_3; typedef typename SNC_structure::Segment_3 Segment_3;
@ -206,90 +194,6 @@ public:
SNC_point_locator_by_spatial_subdivision() : SNC_point_locator_by_spatial_subdivision() :
initialized(false), candidate_provider(0) {} initialized(false), candidate_provider(0) {}
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
template<typename Kernel>
class Triangulation_handler {
typedef typename CGAL::Triangulation_vertex_base_2<Kernel> Vb;
typedef typename CGAL::Constrained_triangulation_face_base_2<Kernel> Fb;
typedef typename CGAL::Triangulation_data_structure_2<Vb,Fb> TDS;
typedef typename CGAL::No_constraint_intersection_requiring_constructions_tag Itag;
typedef typename CGAL::Constrained_triangulation_2<Kernel,TDS,Itag> CT;
typedef typename CT::Face_handle Face_handle;
typedef typename CT::Finite_faces_iterator Finite_face_iterator;
typedef typename CT::Edge Edge;
CT ct;
CGAL::Unique_hash_map<Face_handle, bool> visited;
Finite_face_iterator fi;
public:
template<typename Halffacet_handle>
Triangulation_handler(Halffacet_handle f) : visited(false) {
Halffacet_cycle_iterator fci;
for(fci=f->facet_cycles_begin(); fci!=f->facet_cycles_end(); ++fci) {
if(fci.is_shalfedge()) {
SHalfedge_around_facet_circulator sfc(fci), send(sfc);
CGAL_For_all(sfc,send) {
CGAL_NEF_TRACEN(" insert constraint" << sfc->source()->source()->point()
<< "->" << sfc->source()->twin()->source()->point());
ct.insert_constraint(sfc->source()->source()->point(),
sfc->source()->twin()->source()->point());
}
}
}
CGAL_assertion(ct.is_valid());
CGAL_NEF_TRACEN("number of finite triangles " << ct.number_of_faces());
typename CT::Face_handle infinite = ct.infinite_face();
typename CT::Vertex_handle ctv = infinite->vertex(1);
if(ct.is_infinite(ctv)) ctv = infinite->vertex(2);
CGAL_assertion(!ct.is_infinite(ctv));
typename CT::Face_handle opposite;
typename CT::Face_circulator vc(ctv,infinite);
do { opposite = vc++;
} while(!ct.is_constrained(CT::Edge(vc,vc->index(opposite))));
typename CT::Face_handle first = vc;
CGAL_assertion(!ct.is_infinite(first));
traverse_triangulation(first, first->index(opposite));
/*
for(fi = ct.finite_faces_begin(); fi != ct.finite_faces_end(); ++fi)
CGAL_NEF_TRACEN(" finite face "
<< Triangle_3(fi->vertex(0)->point(), fi->vertex(1)->point(), fi->vertex(2)->point())
<< "was visited " << visited[fi]);
*/
fi = ct.finite_faces_begin();
}
void traverse_triangulation(Face_handle f, int parent) {
visited[f] = true;
if(!ct.is_constrained(Edge(f,ct.cw(parent))) && !visited[f->neighbor(ct.cw(parent))]) {
Face_handle child(f->neighbor(ct.cw(parent)));
traverse_triangulation(child, child->index(f));
}
if(!ct.is_constrained(Edge(f,ct.ccw(parent))) && !visited[f->neighbor(ct.ccw(parent))]) {
Face_handle child(f->neighbor(ct.ccw(parent)));
traverse_triangulation(child, child->index(f));
}
}
template<typename Triangle_3>
bool get_next_triangle(Triangle_3& tr) {
while(fi != ct.finite_faces_end() && visited[fi] == false) ++fi;
if(fi == ct.finite_faces_end()) return false;
tr = Triangle_3(fi->vertex(0)->point(), fi->vertex(1)->point(), fi->vertex(2)->point());
++fi;
return true;
}
};
#endif
virtual void initialize(SNC_structure* W) { virtual void initialize(SNC_structure* W) {
#ifdef CGAL_NEF_LIST_OF_TRIANGLES #ifdef CGAL_NEF_LIST_OF_TRIANGLES
@ -298,11 +202,7 @@ public:
#else // CGAL_NEF_LIST_OF_TRIANGLES #else // CGAL_NEF_LIST_OF_TRIANGLES
CGAL_NEF_TIMER(ct_t.start()); CGAL_NEF_TIMER(ct_t.start());
this->version_ = std::string("Point Locator by Spatial Subdivision (tm)"); this->version_ = std::string("Point Locator by Spatial Subdivision (tm)");
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
CGAL_NEF_CLOG(version()<<" (with triangulated facets)");
#else
CGAL_NEF_CLOG(version()); CGAL_NEF_CLOG(version());
#endif
CGAL_assertion( W != nullptr); CGAL_assertion( W != nullptr);
// (Base) *this = SNC_decorator(*W); // (Base) *this = SNC_decorator(*W);
this->set_snc(*W); this->set_snc(*W);
@ -316,62 +216,7 @@ public:
CGAL_forall_edges( e, *this->sncp()) CGAL_forall_edges( e, *this->sncp())
objects.push_back(make_object(Halfedge_handle(e))); objects.push_back(make_object(Halfedge_handle(e)));
CGAL_forall_facets( f, *this->sncp()) { CGAL_forall_facets( f, *this->sncp()) {
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
#ifndef CGAL_NEF3_TRIANGULATION_MINIMUM
#define CGAL_NEF3_TRIANGULATION_MINIMUM 25
#endif
Halffacet_cycle_iterator fci = f->facet_cycles_begin();
CGAL_assertion(fci.is_shalfedge());
SHalfedge_around_facet_circulator safc(fci), send(safc);
int length = 0;
int stop = CGAL_NEF3_TRIANGULATION_MINIMUM;
while(++length < stop && ++safc != send);
if(length >= stop) {
CGAL_NEF_TRACEN("triangulate facet " << f->plane());
typedef typename CGAL::Projection_traits_xy_3<Kernel> XY;
typedef typename CGAL::Projection_traits_yz_3<Kernel> YZ;
typedef typename CGAL::Projection_traits_xz_3<Kernel> XZ;
Triangle_3 tr;
Vector_3 orth = f->plane().orthogonal_vector();
int c = CGAL::abs(orth[0]) > CGAL::abs(orth[1]) ? 0 : 1;
c = CGAL::abs(orth[2]) > CGAL::abs(orth[c]) ? 2 : c;
if(c == 0) {
Triangulation_handler<YZ> th(f);
while(th.get_next_triangle(tr)) {
Halffacet_triangle_handle th( f, tr);
objects.push_back(make_object(th));
CGAL_NEF_TRACEN("add triangle " << tr);
}
} else if(c == 1) {
Triangulation_handler<XZ> th(f);
while(th.get_next_triangle(tr)) {
Halffacet_triangle_handle th( f, tr);
objects.push_back(make_object(th));
CGAL_NEF_TRACEN("add triangle " << tr);
}
} else if(c == 2) {
Triangulation_handler<XY> th(f);
while(th.get_next_triangle(tr)) {
Halffacet_triangle_handle th( f, tr);
objects.push_back(make_object(th));
CGAL_NEF_TRACEN("add triangle " << tr);
}
} else
CGAL_error_msg( "wrong value");
} else {
CGAL_NEF_TRACEN("add facet " << f->plane());
objects.push_back(make_object(Halffacet_handle(f)));
}
#else
objects.push_back(make_object(Halffacet_handle(f))); objects.push_back(make_object(Halffacet_handle(f)));
#endif
} }
if(initialized) if(initialized)
delete candidate_provider; delete candidate_provider;
@ -538,9 +383,6 @@ public:
Vertex_handle v; Vertex_handle v;
Halfedge_handle e; Halfedge_handle e;
Halffacet_handle f; Halffacet_handle f;
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
Halffacet_triangle_handle t;
#endif
bool hit = false; bool hit = false;
Point_3 eor = CGAL::ORIGIN; // 'end of ray', the latest ray's hit point Point_3 eor = CGAL::ORIGIN; // 'end of ray', the latest ray's hit point
Objects_along_ray objects = candidate_provider->objects_along_ray(ray); Objects_along_ray objects = candidate_provider->objects_along_ray(ray);
@ -605,32 +447,6 @@ public:
_CGAL_NEF_TRACEN("the facet becomes the new hit object"); _CGAL_NEF_TRACEN("the facet becomes the new hit object");
} }
} }
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
else if( CGAL::assign( t, *o) && ((mask&8) != 0)) {
Point_3 q;
Triangle_3 tr = t.get_triangle();
_CGAL_NEF_TRACEN("trying triangle "<<tr);
if( is.does_intersect( ray, tr, q) && normalized(q) != ray.source()) {
_CGAL_NEF_TRACEN("ray intersect triangle on "<< normalized(q));
_CGAL_NEF_TRACEN("prev. intersection? "<<hit);
if( hit) _CGAL_NEF_TRACEN("prev. intersection on "<< normalized(eor));
if( hit && !has_smaller_distance_to_point( ray.source(), q, eor))
continue;
_CGAL_NEF_TRACEN("is the intersection point on the boundary of the facet? "<<
is.does_contain_on_boundary( t, q));
if( is.does_contain_on_boundary( t, q))
continue;
_CGAL_NEF_TRACEN("is the intersection point on the current cell? "<<
candidate_provider->is_point_on_cell(q,objects_iterator));
if( !candidate_provider->is_point_on_cell( q, objects_iterator))
continue;
eor = q;
result = make_object(Halffacet_handle(t));
hit = true;
_CGAL_NEF_TRACEN("the facet becomes the new hit object");
}
}
#endif
else if((mask&15) == 15) else if((mask&15) == 15)
CGAL_error_msg( "wrong handle"); CGAL_error_msg( "wrong handle");
} }
@ -650,9 +466,6 @@ public:
Vertex_handle v; Vertex_handle v;
Halfedge_handle e; Halfedge_handle e;
Halffacet_handle f; Halffacet_handle f;
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
Halffacet_triangle_handle t;
#endif
Object_list candidates = candidate_provider->objects_around_point(p); Object_list candidates = candidate_provider->objects_around_point(p);
Object_list_iterator o = candidates.begin(); Object_list_iterator o = candidates.begin();
bool found = false; bool found = false;
@ -678,20 +491,6 @@ public:
found = true; found = true;
} }
} }
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
else if( CGAL::assign( t, *o)) {
Triangle_3 tr = t.get_triangle();
if( tr.has_on(p)) {
_CGAL_NEF_TRACEN("found on triangle "<<tr);
if(is.does_contain_on_boundary( t, p)) {
_CGAL_NEF_TRACEN("but located on the facet's boundary");
continue;
}
result = make_object(Halffacet_handle(t));
found = true;
}
}
#endif
o++; o++;
} }
if( !found) { if( !found) {
@ -715,9 +514,6 @@ public:
Vertex_handle v; Vertex_handle v;
Halfedge_handle e; Halfedge_handle e;
Halffacet_handle f; Halffacet_handle f;
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
Halffacet_triangle_handle t;
#endif
Object_list candidates = candidate_provider->objects_around_point(p); Object_list candidates = candidate_provider->objects_around_point(p);
Object_list_iterator o = candidates.begin(); Object_list_iterator o = candidates.begin();
@ -818,21 +614,6 @@ public:
s = Segment_3(p, normalized(ip)); s = Segment_3(p, normalized(ip));
result = make_object(f); result = make_object(f);
} }
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
} else if( CGAL::assign(t, *o)) {
CGAL_NEF_TRACEN("test triangle of facet " << t->plane());
Triangle_3 tr = t.get_triangle();
CGAL_NEF_TRACEN("trying triangle "<<tr);
if(tr.has_on(p)) {
f = t;
return make_object(f);
}
if( is.does_intersect( s, tr, ip)) {
CGAL_assertion(!is.does_contain_on_boundary( t, ip));
s = Segment_3(p, normalized(ip));
result = make_object(t);
}
#endif
} }
else CGAL_error_msg( "wrong handle type"); else CGAL_error_msg( "wrong handle type");
} }
@ -895,14 +676,6 @@ public:
if(f->plane().oriented_side(p) == ON_NEGATIVE_SIDE) if(f->plane().oriented_side(p) == ON_NEGATIVE_SIDE)
f = f->twin(); f = f->twin();
return make_object(f->incident_volume()); return make_object(f->incident_volume());
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
} else if( CGAL::assign(t, result)) {
f = t;
_CGAL_NEF_TRACEN("facet hit, obtaining volume...");
if(f->plane().oriented_side(p) == ON_NEGATIVE_SIDE)
f = f->twin();
return make_object(f->incident_volume());
#endif
} else if( CGAL::assign(e, result)) { } else if( CGAL::assign(e, result)) {
SM_decorator SD(&*e->source()); SM_decorator SD(&*e->source());
if( SD.is_isolated(e)) if( SD.is_isolated(e))
@ -922,17 +695,11 @@ public:
_CGAL_NEF_TRACEN( "intersecting edge: "<<&*e0<<' '<<Segment_3(e0->source()->point(), _CGAL_NEF_TRACEN( "intersecting edge: "<<&*e0<<' '<<Segment_3(e0->source()->point(),
e0->twin()->source()->point())); e0->twin()->source()->point()));
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
Unique_hash_map< Halffacet_triangle_handle, bool> f_mark(false);
#endif // CGAL_NEF3_TRIANGULATE_FACETS
Segment_3 s(Segment_3(e0->source()->point(),e0->twin()->source()->point())); Segment_3 s(Segment_3(e0->source()->point(),e0->twin()->source()->point()));
Vertex_handle v; Vertex_handle v;
Halfedge_handle e; Halfedge_handle e;
Halffacet_handle f; Halffacet_handle f;
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
Halffacet_triangle_handle t;
#endif
Object_list_iterator o; Object_list_iterator o;
Object_list objects = candidate_provider->objects_around_segment(s); Object_list objects = candidate_provider->objects_around_segment(s);
CGAL_for_each( o, objects) { CGAL_for_each( o, objects) {
@ -966,23 +733,6 @@ public:
_CGAL_NEF_TRACEN("edge intersects facet on plane "<<f->plane()<<" on "<<q); _CGAL_NEF_TRACEN("edge intersects facet on plane "<<f->plane()<<" on "<<q);
} }
} }
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
else if( CGAL::assign( t, *o)) {
Point_3 q;
Triangle_3 tr = t.get_triangle();
if( f_mark[t])
continue;
_CGAL_NEF_TRACEN("trying with triangle "<<tr);
if( is.does_intersect( s, tr, q) &&
!is.does_contain_on_boundary( t, q)) {
q = normalized(q);
call_back( e0, make_object(Halffacet_handle(t)), q);
_CGAL_NEF_TRACEN("edge intersects facet triangle on plane "<<t->plane()<<" on "<<q);
f_mark[t] = true;
}
}
#endif // CGAL_NEF3_TRIANGULATE_FACETS
else else
CGAL_error_msg( "wrong handle"); CGAL_error_msg( "wrong handle");
} }
@ -999,9 +749,6 @@ public:
Vertex_handle v; Vertex_handle v;
Halfedge_handle e; Halfedge_handle e;
Halffacet_handle f; Halffacet_handle f;
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
Halffacet_triangle_handle t;
#endif
Object_list_iterator o; Object_list_iterator o;
Object_list objects = candidate_provider->objects_around_segment(s); Object_list objects = candidate_provider->objects_around_segment(s);
CGAL_for_each( o, objects) { CGAL_for_each( o, objects) {
@ -1026,11 +773,6 @@ public:
else if( CGAL::assign( f, *o)) { else if( CGAL::assign( f, *o)) {
/* do nothing */ /* do nothing */
} }
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
else if( CGAL::assign( t, *o)) {
/* do nothing */
}
#endif
else else
CGAL_error_msg( "wrong handle"); CGAL_error_msg( "wrong handle");
} }
@ -1043,16 +785,10 @@ public:
CGAL_assertion( initialized); CGAL_assertion( initialized);
_CGAL_NEF_TRACEN( "intersecting edge: "<< Segment_3(e0->source()->point(), _CGAL_NEF_TRACEN( "intersecting edge: "<< Segment_3(e0->source()->point(),
e0->twin()->source()->point())); e0->twin()->source()->point()));
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
Unique_hash_map< Halffacet_triangle_handle, bool> f_mark(false);
#endif // CGAL_NEF3_TRIANGULATE_FACETS
Segment_3 s(Segment_3(e0->source()->point(),e0->twin()->source()->point())); Segment_3 s(Segment_3(e0->source()->point(),e0->twin()->source()->point()));
Vertex_handle v; Vertex_handle v;
Halfedge_handle e; Halfedge_handle e;
Halffacet_handle f; Halffacet_handle f;
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
Halffacet_triangle_handle t;
#endif
Object_list_iterator o; Object_list_iterator o;
Object_list objects = candidate_provider->objects_around_segment(s); Object_list objects = candidate_provider->objects_around_segment(s);
CGAL_for_each( o, objects) { CGAL_for_each( o, objects) {
@ -1075,22 +811,6 @@ public:
_CGAL_NEF_TRACEN("edge intersects facet on plane "<<f->plane()<<" on "<<q); _CGAL_NEF_TRACEN("edge intersects facet on plane "<<f->plane()<<" on "<<q);
} }
} }
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
else if( CGAL::assign( t, *o)) {
Point_3 q;
Triangle_3 tr = t.get_triangle();
if( f_mark[t])
continue;
_CGAL_NEF_TRACEN("trying with triangle "<<tr);
if( is.does_intersect( s, tr, q) &&
!is.does_contain_on_boundary( t, q)) {
q = normalized(q);
call_back( e0, make_object(Halffacet_handle(t)), q);
_CGAL_NEF_TRACEN("edge intersects facet triangle on plane "<<t->plane()<<" on "<<q);
f_mark[t] = true;
}
}
#endif // CGAL_NEF3_TRIANGULATE_FACETS
else else
CGAL_error_msg( "wrong handle"); CGAL_error_msg( "wrong handle");
} }
@ -1148,9 +868,6 @@ private:
bool initialized; bool initialized;
SNC_candidate_provider* candidate_provider; SNC_candidate_provider* candidate_provider;
SNC_intersection is; SNC_intersection is;
#ifdef CGAL_NEF3_TRIANGULATE_FACETS
std::list<Halffacet_triangle_handle> triangulation;
#endif
}; };
#ifdef CGAL_NEF3_POINT_LOCATOR_NAIVE #ifdef CGAL_NEF3_POINT_LOCATOR_NAIVE