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
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
|
// Copyright (c) 2014
// INRIA Saclay-Ile de France (France)
//
// This file is part of CGAL (www.cgal.org)
//
// $URL: https://github.com/CGAL/cgal/blob/v6.1/NewKernel_d/include/CGAL/NewKernel_d/Kernel_d_interface.h $
// $Id: include/CGAL/NewKernel_d/Kernel_d_interface.h b26b07a1242 $
// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-Commercial
//
// Author(s) : Marc Glisse
#ifndef CGAL_KD_KERNEL_D_INTERFACE_H
#define CGAL_KD_KERNEL_D_INTERFACE_H
#include <CGAL/NewKernel_d/functor_tags.h>
#include <CGAL/transforming_iterator.h>
#include <CGAL/NewKernel_d/utils.h>
#include <CGAL/tuple.h>
#include <type_traits>
namespace CGAL {
template <class> struct Kernel_d_interface;
namespace internal::NewKernel_d {
// If the base kernel has approx/exact kernels, also wrap those, for packages like Frechet that want to dig through the internals
template <class K, class = void> struct Kernel_d_interface_filters_ : K { using K::K; };
template <class K> struct Kernel_d_interface_filters_<K, std::void_t<typename K::Approximate_kernel, typename K::Exact_kernel>> : K {
typedef CGAL::Kernel_d_interface<typename K::Approximate_kernel> Approximate_kernel;
typedef CGAL::Kernel_d_interface<typename K::Exact_kernel> Exact_kernel;
using K::K; // inherit constructors
};
}
template <class Base_> struct Kernel_d_interface : public internal::NewKernel_d::Kernel_d_interface_filters_<Base_> {
constexpr Kernel_d_interface(){}
constexpr Kernel_d_interface(int d)
: internal::NewKernel_d::Kernel_d_interface_filters_<Base_>(d) {}
typedef Base_ Base;
typedef Kernel_d_interface<Base> Kernel;
typedef Base_ R_; // for the macros
typedef typename Get_type<Base, RT_tag>::type RT;
typedef typename Get_type<Base, FT_tag>::type FT;
typedef typename Get_type<Base, Bool_tag>::type Boolean;
typedef typename Get_type<Base, Sign_tag>::type Sign;
typedef typename Get_type<Base, Comparison_result_tag>::type Comparison_result;
typedef typename Get_type<Base, Orientation_tag>::type Orientation;
typedef typename Get_type<Base, Oriented_side_tag>::type Oriented_side;
typedef typename Get_type<Base, Bounded_side_tag>::type Bounded_side;
typedef typename Get_type<Base, Angle_tag>::type Angle;
typedef typename Get_type<Base, Flat_orientation_tag>::type Flat_orientation_d;
typedef typename Get_type<Base, Point_tag>::type Point_d;
typedef typename Get_type<Base, Vector_tag>::type Vector_d;
typedef typename Get_type<Base, Segment_tag>::type Segment_d;
typedef typename Get_type<Base, Sphere_tag>::type Sphere_d;
typedef typename Get_type<Base, Hyperplane_tag>::type Hyperplane_d;
typedef Vector_d Direction_d;
typedef typename Get_type<Base, Line_tag>::type Line_d;
typedef typename Get_type<Base, Ray_tag>::type Ray_d;
typedef typename Get_type<Base, Iso_box_tag>::type Iso_box_d;
typedef typename Get_type<Base, Aff_transformation_tag>::type Aff_transformation_d;
typedef typename Get_type<Base, Weighted_point_tag>::type Weighted_point_d;
typedef typename Get_functor<Base, Compute_point_cartesian_coordinate_tag>::type Compute_coordinate_d;
typedef typename Get_functor<Base, Compare_lexicographically_tag>::type Compare_lexicographically_d;
typedef typename Get_functor<Base, Compare_squared_distance_tag>::type Compare_squared_distance_d;
typedef typename Get_functor<Base, Equal_points_tag>::type Equal_d;
typedef typename Get_functor<Base, Less_lexicographically_tag>::type Less_lexicographically_d;
typedef typename Get_functor<Base, Less_or_equal_lexicographically_tag>::type Less_or_equal_lexicographically_d;
// FIXME: and vectors?
typedef typename Get_functor<Base, Orientation_of_points_tag>::type Orientation_d;
typedef typename Get_functor<Base, Less_point_cartesian_coordinate_tag>::type Less_coordinate_d;
typedef typename Get_functor<Base, Point_dimension_tag>::type Point_dimension_d;
typedef typename Get_functor<Base, Side_of_oriented_sphere_tag>::type Side_of_oriented_sphere_d;
typedef typename Get_functor<Base, Power_side_of_power_sphere_tag>::type Power_side_of_power_sphere_d;
typedef typename Get_functor<Base, Power_side_of_bounded_power_circumsphere_tag>::type Power_side_of_bounded_power_sphere_d;
typedef typename Get_functor<Base, Power_center_tag>::type Construct_power_sphere_d;
typedef typename Get_functor<Base, Power_distance_tag>::type Compute_power_product_d;
typedef typename Get_functor<Base, Contained_in_affine_hull_tag>::type Contained_in_affine_hull_d;
typedef typename Get_functor<Base, Construct_flat_orientation_tag>::type Construct_flat_orientation_d;
typedef typename Get_functor<Base, In_flat_orientation_tag>::type In_flat_orientation_d;
typedef typename Get_functor<Base, In_flat_side_of_oriented_sphere_tag>::type In_flat_side_of_oriented_sphere_d;
typedef typename Get_functor<Base, In_flat_power_side_of_power_sphere_tag>::type In_flat_power_side_of_power_sphere_d;
typedef typename Get_functor<Base, Point_to_vector_tag>::type Point_to_vector_d;
typedef typename Get_functor<Base, Vector_to_point_tag>::type Vector_to_point_d;
typedef typename Get_functor<Base, Translated_point_tag>::type Translated_point_d;
typedef typename Get_functor<Base, Scaled_vector_tag>::type Scaled_vector_d;
typedef typename Get_functor<Base, Difference_of_vectors_tag>::type Difference_of_vectors_d;
typedef typename Get_functor<Base, Difference_of_points_tag>::type Difference_of_points_d;
//typedef typename Get_functor<Base, Construct_ttag<Point_tag> >::type Construct_point_d;
struct Construct_point_d : private Store_kernel<Kernel> {
typedef Kernel R_; // for the macro
CGAL_FUNCTOR_INIT_STORE(Construct_point_d)
typedef typename Get_functor<Base, Construct_ttag<Point_tag> >::type CP;
typedef Point_d result_type;
Point_d operator()(Weighted_point_d const&wp)const{
return typename Get_functor<Base, Point_drop_weight_tag>::type(this->kernel())(wp);
}
Point_d operator()(Weighted_point_d &wp)const{
return typename Get_functor<Base, Point_drop_weight_tag>::type(this->kernel())(wp);
}
Point_d operator()(Weighted_point_d &&wp)const{
return typename Get_functor<Base, Point_drop_weight_tag>::type(this->kernel())(std::move(wp));
}
Point_d operator()(Weighted_point_d const&&wp)const{
return typename Get_functor<Base, Point_drop_weight_tag>::type(this->kernel())(std::move(wp));
}
template<class...T>
decltype(auto)
operator()(T&&...t)const{
return CP(this->kernel())(std::forward<T>(t)...);
}
};
//typedef typename Get_functor<Base, Construct_ttag<Vector_tag> >::type Construct_vector_d;
struct Construct_vector_d : private Store_kernel<Kernel> {
typedef Kernel R_; // for the macro
CGAL_FUNCTOR_INIT_STORE(Construct_vector_d)
typedef typename Get_functor<Base, Construct_ttag<Vector_tag> >::type CV;
typedef Vector_d result_type;
Vector_d operator()(Point_d const& p, Point_d const& q)const{
return Difference_of_points_d(this->kernel())(q,p);
}
Vector_d operator()(Point_d & p, Point_d & q)const{
return Difference_of_points_d(this->kernel())(q,p);
}
template<class...T>
decltype(auto)
operator()(T&&...t)const{
return CV(this->kernel())(std::forward<T>(t)...);
}
};
typedef typename Get_functor<Base, Construct_ttag<Segment_tag> >::type Construct_segment_d;
typedef typename Get_functor<Base, Construct_ttag<Sphere_tag> >::type Construct_sphere_d;
typedef typename Get_functor<Base, Construct_ttag<Hyperplane_tag> >::type Construct_hyperplane_d;
typedef Construct_vector_d Construct_direction_d;
typedef typename Get_functor<Base, Construct_ttag<Line_tag> >::type Construct_line_d;
typedef typename Get_functor<Base, Construct_ttag<Ray_tag> >::type Construct_ray_d;
typedef typename Get_functor<Base, Construct_ttag<Iso_box_tag> >::type Construct_iso_box_d;
typedef typename Get_functor<Base, Construct_ttag<Aff_transformation_tag> >::type Construct_aff_transformation_d;
typedef typename Get_functor<Base, Construct_ttag<Weighted_point_tag> >::type Construct_weighted_point_d;
typedef typename Get_functor<Base, Midpoint_tag>::type Midpoint_d;
struct Component_accessor_d : private Store_kernel<Kernel> {
typedef Kernel R_; // for the macro
CGAL_FUNCTOR_INIT_STORE(Component_accessor_d)
int dimension(Point_d const&p){
return this->kernel().point_dimension_d_object()(p);
}
FT cartesian(Point_d const&p, int i){
return this->kernel().compute_coordinate_d_object()(p,i);
}
RT homogeneous(Point_d const&p, int i){
if (i == dimension(p))
return 1;
return cartesian(p, i);
}
};
struct Construct_cartesian_const_iterator_d : private Store_kernel<R_> {
CGAL_FUNCTOR_INIT_STORE(Construct_cartesian_const_iterator_d)
typedef typename Get_functor<Base, Construct_ttag<Point_cartesian_const_iterator_tag> >::type CPI;
typedef typename Get_functor<Base, Construct_ttag<Vector_cartesian_const_iterator_tag> >::type CVI;
// @mglisse shall we update that code ?
// FIXME: The following sometimes breaks compilation. The typedef below forces instantiation of this, which forces Point_d, which itself (in the wrapper) needs the derived kernel to tell it what the base kernel is, and that's a cycle. The exact circumstances are not clear, g++ and clang++ are ok in both C++03 and C++11, it is only clang in C++11 without CGAL_CXX11 that breaks. Relying on CPI::result_type is great for Epick_d but not Epeck_d.
// typedef typename CGAL::decay<typename boost::result_of<CPI(Point_d,CGAL::Begin_tag)>::type>::type result_type;
// typedef typename CGAL::decay<typename CPI::result_type>::type result_type;
// typedef decltype(std::declval<CPI>()(std::declval<Point_d>(),Begin_tag{})) result_type;
// HACK
typedef typename Base::Point_cartesian_const_iterator result_type;
// Kernel_d requires a common iterator type for points and vectors
// TODO: provide this mixed functor in preKernel?
// static_assert(std::is_same<typename CGAL::decay<typename boost::result_of<CVI(Vector_d,CGAL::Begin_tag)>::type>::type, result_type>::value);
// static_assert(std::is_same<typename CGAL::decay<typename CVI::result_type>::type, result_type>::value);
template <class Tag_>
auto operator()(Point_d const&p, Tag_ t)const{
return CPI(this->kernel())(p,t);
}
template <class Tag_>
auto operator()(typename First_if_different<Vector_d,Point_d>::Type const&v, Tag_ t)const{
return CVI(this->kernel())(v,t);
}
template <class Obj>
auto operator()(Obj const&o)const{
return operator()(o, Begin_tag());
}
auto operator()(Point_d const&p, int)const{
return operator()(p, End_tag());
}
auto operator()(typename First_if_different<Vector_d,Point_d>::Type const&v, int)const{
return operator()(v, End_tag());
}
};
struct Compute_squared_radius_d : private Store_kernel<Kernel> {
typedef Kernel R_; // for the macro
CGAL_FUNCTOR_INIT_STORE(Compute_squared_radius_d)
typedef FT result_type;
template<class S> FT operator()(S&& s)const{
return typename Get_functor<Base, Squared_radius_tag>::type(this->kernel())(std::forward<S>(s));
}
template<class I> FT operator()(I b, I e)const{
return typename Get_functor<Base, Squared_circumradius_tag>::type(this->kernel())(b,e);
}
};
struct Compute_squared_radius_smallest_orthogonal_sphere_d : private Store_kernel<Kernel> {
typedef Kernel R_; // for the macro
CGAL_FUNCTOR_INIT_STORE(Compute_squared_radius_smallest_orthogonal_sphere_d)
typedef FT result_type;
template<class I> FT operator()(I b, I e)const{
typename Get_functor<Base, Point_weight_tag>::type pw(this->kernel());
typename Get_functor<Base, Power_center_tag>::type pc(this->kernel());
return pw(pc(b,e));
}
};
typedef typename Construct_cartesian_const_iterator_d::result_type Cartesian_const_iterator_d;
typedef typename Get_functor<Base, Squared_distance_tag>::type Squared_distance_d;
typedef typename Get_functor<Base, Squared_length_tag>::type Squared_length_d;
typedef typename Get_functor<Base, Construct_bbox_tag>::type Construct_bbox_d;
typedef typename Get_functor<Base, Scalar_product_tag>::type Scalar_product_d;
typedef typename Get_functor<Base, Affine_rank_tag>::type Affine_rank_d;
typedef typename Get_functor<Base, Affinely_independent_tag>::type Affinely_independent_d;
typedef typename Get_functor<Base, Contained_in_linear_hull_tag>::type Contained_in_linear_hull_d;
typedef typename Get_functor<Base, Contained_in_simplex_tag>::type Contained_in_simplex_d;
typedef typename Get_functor<Base, Has_on_positive_side_tag>::type Has_on_positive_side_d;
typedef typename Get_functor<Base, Linear_rank_tag>::type Linear_rank_d;
typedef typename Get_functor<Base, Linearly_independent_tag>::type Linearly_independent_d;
typedef typename Get_functor<Base, Oriented_side_tag>::type Oriented_side_d;
typedef typename Get_functor<Base, Side_of_bounded_circumsphere_tag>::type Side_of_bounded_sphere_d;
typedef typename Get_functor<Base, Center_of_sphere_tag>::type Center_of_sphere_d;
typedef Center_of_sphere_d Construct_center_d; // RangeSearchTraits
typedef typename Get_functor<Base, Construct_circumcenter_tag>::type Construct_circumcenter_d;
typedef typename Get_functor<Base, Value_at_tag>::type Value_at_d;
typedef typename Get_functor<Base, Point_of_sphere_tag>::type Point_of_sphere_d;
typedef typename Get_functor<Base, Orthogonal_vector_tag>::type Orthogonal_vector_d;
typedef typename Get_functor<Base, Linear_base_tag>::type Linear_base_d;
typedef typename Get_functor<Base, Construct_min_vertex_tag>::type Construct_min_vertex_d;
typedef typename Get_functor<Base, Construct_max_vertex_tag>::type Construct_max_vertex_d;
typedef typename Get_functor<Base, Point_weight_tag>::type Compute_weight_d;
typedef typename Get_functor<Base, Point_drop_weight_tag>::type Point_drop_weight_d;
//TODO:
//typedef ??? Intersect_d;
Compute_coordinate_d compute_coordinate_d_object()const{ return Compute_coordinate_d(*this); }
Has_on_positive_side_d has_on_positive_side_d_object()const{ return Has_on_positive_side_d(*this); }
Compare_lexicographically_d compare_lexicographically_d_object()const{ return Compare_lexicographically_d(*this); }
Compare_squared_distance_d compare_squared_distance_d_object()const{ return Compare_squared_distance_d(*this); }
Equal_d equal_d_object()const{ return Equal_d(*this); }
Less_lexicographically_d less_lexicographically_d_object()const{ return Less_lexicographically_d(*this); }
Less_or_equal_lexicographically_d less_or_equal_lexicographically_d_object()const{ return Less_or_equal_lexicographically_d(*this); }
Less_coordinate_d less_coordinate_d_object()const{ return Less_coordinate_d(*this); }
Orientation_d orientation_d_object()const{ return Orientation_d(*this); }
Oriented_side_d oriented_side_d_object()const{ return Oriented_side_d(*this); }
Point_dimension_d point_dimension_d_object()const{ return Point_dimension_d(*this); }
Point_of_sphere_d point_of_sphere_d_object()const{ return Point_of_sphere_d(*this); }
Side_of_oriented_sphere_d side_of_oriented_sphere_d_object()const{ return Side_of_oriented_sphere_d(*this); }
Power_side_of_power_sphere_d power_side_of_power_sphere_d_object()const{ return Power_side_of_power_sphere_d(*this); }
Power_side_of_bounded_power_sphere_d power_side_of_bounded_power_sphere_d_object()const{ return Power_side_of_bounded_power_sphere_d(*this); }
Construct_power_sphere_d construct_power_sphere_d_object()const{ return Construct_power_sphere_d(*this); }
Compute_power_product_d compute_power_product_d_object()const{ return Compute_power_product_d(*this); }
Side_of_bounded_sphere_d side_of_bounded_sphere_d_object()const{ return Side_of_bounded_sphere_d(*this); }
Contained_in_affine_hull_d contained_in_affine_hull_d_object()const{ return Contained_in_affine_hull_d(*this); }
Contained_in_linear_hull_d contained_in_linear_hull_d_object()const{ return Contained_in_linear_hull_d(*this); }
Contained_in_simplex_d contained_in_simplex_d_object()const{ return Contained_in_simplex_d(*this); }
Construct_flat_orientation_d construct_flat_orientation_d_object()const{ return Construct_flat_orientation_d(*this); }
In_flat_orientation_d in_flat_orientation_d_object()const{ return In_flat_orientation_d(*this); }
In_flat_side_of_oriented_sphere_d in_flat_side_of_oriented_sphere_d_object()const{ return In_flat_side_of_oriented_sphere_d(*this); }
In_flat_power_side_of_power_sphere_d in_flat_power_side_of_power_sphere_d_object()const{ return In_flat_power_side_of_power_sphere_d(*this); }
Point_to_vector_d point_to_vector_d_object()const{ return Point_to_vector_d(*this); }
Vector_to_point_d vector_to_point_d_object()const{ return Vector_to_point_d(*this); }
Translated_point_d translated_point_d_object()const{ return Translated_point_d(*this); }
Scaled_vector_d scaled_vector_d_object()const{ return Scaled_vector_d(*this); }
Difference_of_vectors_d difference_of_vectors_d_object()const{ return Difference_of_vectors_d(*this); }
Difference_of_points_d difference_of_points_d_object()const{ return Difference_of_points_d(*this); }
Affine_rank_d affine_rank_d_object()const{ return Affine_rank_d(*this); }
Affinely_independent_d affinely_independent_d_object()const{ return Affinely_independent_d(*this); }
Linear_base_d linear_base_d_object()const{ return Linear_base_d(*this); }
Linear_rank_d linear_rank_d_object()const{ return Linear_rank_d(*this); }
Linearly_independent_d linearly_independent_d_object()const{ return Linearly_independent_d(*this); }
Midpoint_d midpoint_d_object()const{ return Midpoint_d(*this); }
Value_at_d value_at_d_object()const{ return Value_at_d(*this); }
/// Intersect_d intersect_d_object()const{ return Intersect_d(*this); }
Component_accessor_d component_accessor_d_object()const{ return Component_accessor_d(*this); }
Orthogonal_vector_d orthogonal_vector_d_object()const{ return Orthogonal_vector_d(*this); }
Construct_cartesian_const_iterator_d construct_cartesian_const_iterator_d_object()const{ return Construct_cartesian_const_iterator_d(*this); }
Construct_point_d construct_point_d_object()const{ return Construct_point_d(*this); }
Construct_vector_d construct_vector_d_object()const{ return Construct_vector_d(*this); }
Construct_segment_d construct_segment_d_object()const{ return Construct_segment_d(*this); }
Construct_sphere_d construct_sphere_d_object()const{ return Construct_sphere_d(*this); }
Construct_hyperplane_d construct_hyperplane_d_object()const{ return Construct_hyperplane_d(*this); }
Compute_squared_radius_d compute_squared_radius_d_object()const{ return Compute_squared_radius_d(*this); }
Compute_squared_radius_smallest_orthogonal_sphere_d compute_squared_radius_smallest_orthogonal_sphere_d_object()const{ return Compute_squared_radius_smallest_orthogonal_sphere_d(*this); }
Squared_distance_d squared_distance_d_object()const{ return Squared_distance_d(*this); }
Squared_length_d squared_length_d_object()const{ return Squared_length_d(*this); }
Construct_bbox_d construct_bbox_d_object()const{ return Construct_bbox_d(*this); }
Scalar_product_d scalar_product_d_object()const{ return Scalar_product_d(*this); }
Center_of_sphere_d center_of_sphere_d_object()const{ return Center_of_sphere_d(*this); }
Construct_circumcenter_d construct_circumcenter_d_object()const{ return Construct_circumcenter_d(*this); }
Construct_direction_d construct_direction_d_object()const{ return Construct_direction_d(*this); }
Construct_line_d construct_line_d_object()const{ return Construct_line_d(*this); }
Construct_ray_d construct_ray_d_object()const{ return Construct_ray_d(*this); }
Construct_iso_box_d construct_iso_box_d_object()const{ return Construct_iso_box_d(*this); }
Construct_aff_transformation_d construct_aff_transformation_d_object()const{ return Construct_aff_transformation_d(*this); }
Construct_min_vertex_d construct_min_vertex_d_object()const{ return Construct_min_vertex_d(*this); }
Construct_max_vertex_d construct_max_vertex_d_object()const{ return Construct_max_vertex_d(*this); }
Construct_weighted_point_d construct_weighted_point_d_object()const{ return Construct_weighted_point_d(*this); }
Compute_weight_d compute_weight_d_object()const{ return Compute_weight_d(*this); }
Point_drop_weight_d point_drop_weight_d_object()const{ return Point_drop_weight_d(*this); }
// Dummies for those required functors missing a concept.
typedef Null_functor Position_on_line_d;
Position_on_line_d position_on_line_d_object()const{return Null_functor();}
typedef Null_functor Barycentric_coordinates_d;
Barycentric_coordinates_d barycentric_coordinates_d_object()const{return Null_functor();}
/* Not provided because they don't make sense here:
Lift_to_paraboloid_d
Project_along_d_axis_d
*/
};
}
#endif // CGAL_KD_KERNEL_D_INTERFACE_H
|
