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
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
|
// Copyright (c) 2009 INRIA Sophia-Antipolis (France).
// All rights reserved.
//
// This file is part of CGAL (www.cgal.org).
//
// $URL: https://github.com/CGAL/cgal/blob/v6.1/Mesh_3/include/CGAL/Mesh_3/Robust_intersection_kernel.h $
// $Id: include/CGAL/Mesh_3/Robust_intersection_kernel.h b26b07a1242 $
// SPDX-License-Identifier: GPL-3.0-or-later OR LicenseRef-Commercial
//
// Author(s) : Stephane Tayeb
//
//******************************************************************************
// File Description :
//******************************************************************************
#ifndef CGAL_MESH_3_ROBUST_INTERSECTION_KERNEL_3_H
#define CGAL_MESH_3_ROBUST_INTERSECTION_KERNEL_3_H
#include <CGAL/license/Mesh_3.h>
#include <CGAL/Cartesian_converter.h>
#include <CGAL/Exact_predicates_exact_constructions_kernel.h>
#include <CGAL/Simple_cartesian.h>
#include <CGAL/Gmpq.h>
namespace CGAL {
namespace Mesh_3 {
template < typename K_ >
class Robust_intersection_for_kernel_3
{
public:
typedef typename K_::Line_3 Line_3;
typedef typename K_::Plane_3 Plane_3;
typedef Object result_type;
typedef Robust_intersection_for_kernel_3<K_> Self;
typedef Exact_predicates_exact_constructions_kernel EK;
// typedef Simple_cartesian<Gmpq> EK;
typedef Cartesian_converter<typename K_::Kernel, EK> To_exact;
typedef Cartesian_converter<EK, typename K_::Kernel> Back_from_exact;
// forward to the base
template<class T1, class T2>
Object operator() (const T1& t, const T2& s) const
{
return CGAL::CommonKernelFunctors::Intersect_3<K_>()(t, s);
}
// exact computation only for Intersect_3()(Line_3, Plane_3)
Object operator() (const Line_3& line, const Plane_3& plane) const
{
// Switch to exact
To_exact to_exact;
Back_from_exact back_from_exact;
EK::Intersect_3 exact_intersection = EK().intersect_3_object();
Object object = exact_intersection(to_exact(line), to_exact(plane));
// std::exit(1);
if ( const EK::Point_3* p = object_cast<EK::Point_3>(&object) )
return make_object(back_from_exact(*p));
else if ( const EK::Segment_3* seg = object_cast<EK::Segment_3>(&object) )
return make_object(back_from_exact(*seg));
else
return Object();
}
}; // end template Robust_intersection_for_kernel_3
template <typename K_base, typename Kernel>
struct Robust_intersection_kernel_base
: public K_base::template Base<Kernel>::Type
{
// template < typename Kernel2 >
// struct Base {
// typedef typename K_base::template Base<Kernel2> K2;
// typedef Robust_intersection_kernel_base<K2> Type;
// };
typedef Robust_intersection_for_kernel_3<Kernel> Intersect_3;
Intersect_3
intersect_3_object() const
{
return Intersect_3();
}
}; // end template Robust_intersection_kernel_base
template <typename K>
struct Robust_intersection_kernel
: Type_equality_wrapper<Robust_intersection_kernel_base<K,
Robust_intersection_kernel<K> >,
Robust_intersection_kernel<K> >
{};
} // end namespace Mesh_3
} //namespace CGAL
#endif // CGAL_MESH_3_ROBUST_INTERSECTION_KERNEL_3_H
|
