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
|
// Boost.Geometry
// Copyright (c) 2015, 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 <cstddef>
#include <algorithm>
#include <boost/type_traits/is_floating_point.hpp>
#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,
bool Enable = boost::is_floating_point
<
typename coordinate_type<Point>::type
>::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 eps = math::scaled_epsilon(coord);
set<I>(point, PlusOrMinus<coord_type>()(coord, eps));
corner_by_epsilon<Point, PlusOrMinus, I+1>::apply(point);
}
};
template
<
typename Point,
template <typename> class PlusOrMinus,
std::size_t I,
std::size_t D
>
struct corner_by_epsilon<Point, PlusOrMinus, I, D, false>
{
static inline void apply(Point const&) {}
};
template
<
typename Point,
template <typename> class PlusOrMinus,
std::size_t D,
bool Enable
>
struct corner_by_epsilon<Point, PlusOrMinus, D, D, Enable>
{
static inline void apply(Point const&) {}
};
template
<
typename Point,
template <typename> class PlusOrMinus,
std::size_t D
>
struct corner_by_epsilon<Point, PlusOrMinus, D, D, false>
{
static inline void apply(Point const&) {}
};
} // namespace expand
template <typename Box>
inline void expand_by_epsilon(Box & box)
{
typedef detail::indexed_point_view<Box, min_corner> min_type;
min_type min_point(box);
expand::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);
expand::corner_by_epsilon<max_type, std::plus>::apply(max_point);
}
} // namespace detail
#endif // DOXYGEN_NO_DETAIL
}} // namespace boost::geometry
#endif // BOOST_GEOMETRY_ALGORITHMS_DETAIL_EXPAND_EXPAND_BY_EPSILON_HPP
|