1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
|
#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/Exact_predicates_exact_constructions_kernel.h>
#include <CGAL/Delaunay_triangulation_3.h>
#include <CGAL/point_generators_3.h>
#include <cassert>
typedef CGAL::Exact_predicates_inexact_constructions_kernel EPIC;
typedef CGAL::Exact_predicates_exact_constructions_kernel EPEC;
typedef CGAL::Delaunay_triangulation_3< EPIC > DT3_epic;
typedef CGAL::Delaunay_triangulation_3< EPEC > DT3_epec;
typedef CGAL::Cartesian_converter<
CGAL::Exact_predicates_inexact_constructions_kernel,
CGAL::Exact_predicates_exact_constructions_kernel > Converter;
typedef CGAL::Creator_uniform_3<double,EPIC::Point_3> Creator;
struct Convert_vertex{
mutable bool first_vertex;
Convert_vertex():first_vertex(true) {}
DT3_epec::Vertex operator()(const DT3_epic::Vertex&) const { return DT3_epec::Vertex(); }
void operator()(const DT3_epic::Vertex& src,DT3_epec::Vertex& tgt) const
{
if (!first_vertex)
tgt.point() = Converter()( src.point() );
else
first_vertex=false;
}
};
struct Convert_cell{
DT3_epec::Cell operator()(const DT3_epic::Cell&) const { return DT3_epec::Cell(); }
void operator()(const DT3_epic::Cell&,DT3_epec::Cell&) const {}
};
int main()
{
std::vector< EPIC::Point_3> points;
CGAL::Random_points_in_sphere_3<EPIC::Point_3,Creator> g(1.0);
std::copy_n( g, 600, std::back_inserter(points) );
DT3_epic dt3_epic(points.begin(), points.end());
DT3_epec dt3_epec;
dt3_epec.set_infinite_vertex(
dt3_epec.tds().copy_tds( dt3_epic.tds(),dt3_epic.infinite_vertex(), Convert_vertex(), Convert_cell() ) );
assert( dt3_epec.is_valid() );
}
|
