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
|
// Boost.Geometry
// Copyright (c) 2015-2020, Oracle and/or its affiliates.
// Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
// Distributed under the Boost Software License, Version 1.0.
// (See accompanying file LICENSE_1_0.txt or copy at
// http://www.boost.org/LICENSE_1_0.txt)
#ifndef BOOST_GEOMETRY_ALGORITHMS_DETAIL_EXPAND_EXPAND_BY_EPSILON_HPP
#define BOOST_GEOMETRY_ALGORITHMS_DETAIL_EXPAND_EXPAND_BY_EPSILON_HPP
#include <algorithm>
#include <cstddef>
#include <type_traits>
#include <boost/geometry/core/access.hpp>
#include <boost/geometry/core/coordinate_dimension.hpp>
#include <boost/geometry/core/coordinate_type.hpp>
#include <boost/geometry/util/math.hpp>
#include <boost/geometry/views/detail/indexed_point_view.hpp>
namespace boost { namespace geometry
{
#ifndef DOXYGEN_NO_DETAIL
namespace detail { namespace expand
{
template
<
typename Point,
template <typename> class PlusOrMinus,
std::size_t I = 0,
std::size_t D = dimension<Point>::value
>
struct corner_by_epsilon
{
static inline void apply(Point & point)
{
typedef typename coordinate_type<Point>::type coord_type;
coord_type const coord = get<I>(point);
coord_type const seps = math::scaled_epsilon(coord);
set<I>(point, PlusOrMinus<coord_type>()(coord, seps));
corner_by_epsilon<Point, PlusOrMinus, I+1>::apply(point);
}
static inline void apply(Point & point,
typename coordinate_type<Point>::type const& eps)
{
typedef typename coordinate_type<Point>::type coord_type;
coord_type const coord = get<I>(point);
coord_type const seps = math::scaled_epsilon(coord, eps);
set<I>(point, PlusOrMinus<coord_type>()(coord, seps));
corner_by_epsilon<Point, PlusOrMinus, I + 1>::apply(point);
}
};
template
<
typename Point,
template <typename> class PlusOrMinus,
std::size_t D
>
struct corner_by_epsilon<Point, PlusOrMinus, D, D>
{
static inline void apply(Point const&) {}
static inline void apply(Point const&, typename coordinate_type<Point>::type const&) {}
};
template
<
typename Box,
bool Enable = ! std::is_integral<typename coordinate_type<Box>::type>::value
>
struct expand_by_epsilon
{
static inline void apply(Box & box)
{
typedef detail::indexed_point_view<Box, min_corner> min_type;
min_type min_point(box);
corner_by_epsilon<min_type, std::minus>::apply(min_point);
typedef detail::indexed_point_view<Box, max_corner> max_type;
max_type max_point(box);
corner_by_epsilon<max_type, std::plus>::apply(max_point);
}
static inline void apply(Box & box,
typename coordinate_type<Box>::type const& eps)
{
typedef detail::indexed_point_view<Box, min_corner> min_type;
min_type min_point(box);
corner_by_epsilon<min_type, std::minus>::apply(min_point, eps);
typedef detail::indexed_point_view<Box, max_corner> max_type;
max_type max_point(box);
corner_by_epsilon<max_type, std::plus>::apply(max_point, eps);
}
};
template <typename Box>
struct expand_by_epsilon<Box, false>
{
static inline void apply(Box &) {}
static inline void apply(Box &, typename coordinate_type<Box>::type const&) {}
};
} // namespace expand
template <typename Box>
inline void expand_by_epsilon(Box & box)
{
expand::expand_by_epsilon<Box>::apply(box);
}
template <typename Box>
inline void expand_by_epsilon(Box & box,
typename coordinate_type<Box>::type const& eps)
{
expand::expand_by_epsilon<Box>::apply(box, eps);
}
} // namespace detail
#endif // DOXYGEN_NO_DETAIL
}} // namespace boost::geometry
#endif // BOOST_GEOMETRY_ALGORITHMS_DETAIL_EXPAND_EXPAND_BY_EPSILON_HPP
|
