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

added test with degenerate cases for HT2

This commit is contained in:
Iordan Iordanov 2019-02-13 23:36:19 +01:00
parent 5731a61554
commit 7eb7e44b1a
2 changed files with 92 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
Hyperbolic_triangulation_2/test/Hyperbolic_triangulation_2/CMakeLists.txt
View File

@ -20,6 +20,7 @@ if ( CGAL_FOUND )
create_single_source_cgal_program( "ht2_test_swap.cpp" )
create_single_source_cgal_program( "ht2_test_copy.cpp" )
create_single_source_cgal_program( "ht2_test_hyperbolic_circulator.cpp" )
create_single_source_cgal_program( "ht2_test_insert_degenerate.cpp" )
else()

91
Hyperbolic_triangulation_2/test/Hyperbolic_triangulation_2/ht2_test_insert_degenerate.cpp Normal file
View File

@ -0,0 +1,91 @@
#include <CGAL/Hyperbolic_Delaunay_triangulation_2.h>
#include <CGAL/Hyperbolic_Delaunay_triangulation_traits_2.h>
typedef CGAL::Hyperbolic_Delaunay_triangulation_traits_2<> Traits;
typedef CGAL::Hyperbolic_Delaunay_triangulation_2<Traits> HDTriangulation;
typedef HDTriangulation::Point Point;
typedef Traits::FT FT;
typedef Traits::Circle_2 Circle;
typedef Traits::Euclidean_line_2 ELine;
typedef Traits::Construct_hyperbolic_segment_2 CHSegment;
typedef Traits::Construct_intersection_2 CIntersection;
int main()
{
std::vector<Point> P(10), Pp(10);
std::vector<ELine> S(10), Sp(10);
P[0] = Point( FT(2)/FT(10), FT(8)/FT(10));
Pp[0] = Point( -FT(2)/FT(10), -FT(8)/FT(10));
P[1] = Point( FT(4)/FT(10), FT(8)/FT(10));
Pp[1] = Point( -FT(4)/FT(10), -FT(8)/FT(10));
P[2] = Point( FT(6)/FT(10), FT(6)/FT(10));
Pp[2] = Point( -FT(6)/FT(10), -FT(6)/FT(10));
P[3] = Point( FT(8)/FT(10), FT(4)/FT(10));
Pp[3] = Point( -FT(8)/FT(10), -FT(4)/FT(10));
P[4] = Point( FT(8)/FT(10), FT(2)/FT(10));
Pp[4] = Point( -FT(8)/FT(10), -FT(2)/FT(10));
P[5] = Point( FT(8)/FT(10), -FT(2)/FT(10));
Pp[5] = Point( -FT(8)/FT(10), FT(2)/FT(10));
P[6] = Point( FT(8)/FT(10), -FT(4)/FT(10));
Pp[6] = Point( -FT(8)/FT(10), FT(4)/FT(10));
P[7] = Point( FT(6)/FT(10), -FT(6)/FT(10));
Pp[7] = Point( -FT(6)/FT(10), FT(6)/FT(10));
P[8] = Point( FT(4)/FT(10), -FT(8)/FT(10));
Pp[9] = Point( -FT(4)/FT(10), FT(8)/FT(10));
P[9] = Point( FT(2)/FT(10), -FT(8)/FT(10));
Pp[9] = Point( -FT(2)/FT(10), FT(8)/FT(10));
Point O(FT(0), FT(0));
for (int i = 0; i < 10; i++) {
S[i] = ELine( P[i], Pp[i] );
}
HDTriangulation tri;
Circle c1( O, FT(1)/FT(3) );
for (int i = 0; i < 10; i++) {
std::cout << "--- c1, i = " << i << std::endl;
std::pair<Point, Point> res = CIntersection()(c1, S[i]);
tri.insert(res.first);
tri.insert(res.second);
}
Circle c2( O, FT(1)/FT(2) );
for (int i = 0; i < 10; i++) {
std::cout << "--- c2, i = " << i << std::endl;
std::pair<Point, Point> res = CIntersection()(c2, S[i]);
tri.insert(res.first);
tri.insert(res.second);
}
Circle c3( O, FT(2)/FT(3) );
for (int i = 0; i < 10; i++) {
std::cout << "--- c3, i = " << i << std::endl;
std::pair<Point, Point> res = CIntersection()(c3, S[i]);
tri.insert(res.first);
tri.insert(res.second);
}
std::cout << " -------- inserting --------" << std::endl;
std::cout << "Vertices in triangulation: " << tri.number_of_vertices() << std::endl;
std::cout << "Faces in triangulation: " << tri.number_of_hyperbolic_faces() << std::endl;
std::cout << "Dimension of triangulation: " << tri.dimension() << std::endl;
assert(tri.dimension() == 2);
assert(tri.is_valid());
std::cout << " -------- SUCCESS --------" << std::endl;
return 0;
}