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

within a macro block called CGAL_SDG_CHECK_INCIRCLE_CONSISTENCY, 

we check if the results of the old and new incircle tests are the same
This commit is contained in:
Menelaos Karavelas 2010-03-08 19:24:42 +00:00
parent d781e8698a
commit d672394d1a
1 changed files with 28 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
Segment_Delaunay_graph_2/include/CGAL/Segment_Delaunay_graph_2/Vertex_conflict_C2.h
View File

@ -26,6 +26,11 @@
#include <CGAL/Segment_Delaunay_graph_2/Are_same_points_C2.h>
#include <CGAL/Segment_Delaunay_graph_2/Are_same_segments_C2.h>
#ifdef CGAL_SDG_CHECK_INCIRCLE_CONSISTENCY
#include <CGAL/Segment_Delaunay_graph_2/Incircle_sqrt_field_C2.h>
#include <CGAL/Segment_Delaunay_graph_2/Incircle_operator_sqrt_field_C2.h>
#endif
CGAL_BEGIN_NAMESPACE
CGAL_SEGMENT_DELAUNAY_GRAPH_2_BEGIN_NAMESPACE
@ -479,9 +484,32 @@ public:
}
#endif
#ifdef CGAL_SDG_CHECK_INCIRCLE_CONSISTENCY
typedef Incircle_sqrt_field_C2<K> Alternative1_Voronoi_vertex_2;
typedef Incircle_operator_sqrt_field_C2<K> Alternative2_Voronoi_vertex_2;
Voronoi_vertex_2 v(p, q, r);
Alternative1_Voronoi_vertex_2 valt1(p, q, r);
Alternative2_Voronoi_vertex_2 valt2(p, q, r);
Sign s = v.incircle(t);
Sign salt1 = valt1.incircle(t);
Sign salt2 = valt2.incircle(t);
if ( s != salt1 || s != salt2 ) {
std::cerr << "different results" << std::endl;
std::cerr << p << std::endl;
std::cerr << q << std::endl;
std::cerr << r << std::endl;
std::cerr << t << std::endl;
}
return s;
#else
Voronoi_vertex_2 v(p, q, r);
return v.incircle(t);
#endif
}