1 // Boost.Geometry
2
3 // Copyright (c) 2019, Oracle and/or its affiliates.
4
5 // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
6
7 // Licensed under the Boost Software License version 1.0.
8 // http://www.boost.org/users/license.html
9
10 #include "common.hpp"
11
12 #include <boost/geometry/algorithms/difference.hpp>
13 #include <boost/geometry/algorithms/intersection.hpp>
14 #include <boost/geometry/algorithms/sym_difference.hpp>
15 #include <boost/geometry/algorithms/union.hpp>
16
17 template <typename G1, typename G2, typename G3, typename S>
set_idsu(G1 const & g1,G2 const & g2,G3 & g3,S const & s)18 inline void set_idsu(G1 const& g1, G2 const& g2, G3 & g3, S const& s)
19 {
20 bg::intersection(g1, g2, g3, s);
21 bg::difference(g1, g2, g3, s);
22 bg::sym_difference(g1, g2, g3, s);
23 bg::union_(g1, g2, g3, s);
24 }
25
26 template <typename G1, typename G2, typename G3, typename S>
set_ids(G1 const & g1,G2 const & g2,G3 & g3,S const & s)27 inline void set_ids(G1 const& g1, G2 const& g2, G3 & g3, S const& s)
28 {
29 bg::intersection(g1, g2, g3, s);
30 bg::difference(g1, g2, g3, s);
31 bg::sym_difference(g1, g2, g3, s);
32 }
33
34 template <typename G1, typename G2, typename G3, typename S>
set_id(G1 const & g1,G2 const & g2,G3 & g3,S const & s)35 inline void set_id(G1 const& g1, G2 const& g2, G3 & g3, S const& s)
36 {
37 bg::intersection(g1, g2, g3, s);
38 bg::difference(g1, g2, g3, s);
39 }
40
41 template <typename G1, typename G2, typename G3, typename S>
set_i(G1 const & g1,G2 const & g2,G3 & g3,S const & s)42 inline void set_i(G1 const& g1, G2 const& g2, G3 & g3, S const& s)
43 {
44 bg::intersection(g1, g2, g3, s);
45 }
46
47 template <typename G1, typename G2, typename G3, typename S>
set_d(G1 const & g1,G2 const & g2,G3 & g3,S const & s)48 inline void set_d(G1 const& g1, G2 const& g2, G3 & g3, S const& s)
49 {
50 bg::difference(g1, g2, g3, s);
51 }
52
53 template <typename G1, typename G2, typename G3>
set_idsu_pp(G1 const & g1,G2 const & g2,G3 & g3)54 inline void set_idsu_pp(G1 const& g1, G2 const& g2, G3 & g3)
55 {
56 ::set_idsu(g1, g2, g3, bg::strategy::within::cartesian_point_point());
57 ::set_idsu(g1, g2, g3, bg::strategy::within::spherical_point_point());
58 }
59
60 template <typename G1, typename G2, typename G3>
set_idsu_ps(G1 const & g1,G2 const & g2,G3 & g3)61 inline void set_idsu_ps(G1 const& g1, G2 const& g2, G3 & g3)
62 {
63 typedef typename bg::point_type<G1>::type point_type;
64 ::set_idsu(g1, g2, g3, bg::strategy::within::cartesian_winding<point_type>());
65 ::set_idsu(g1, g2, g3, bg::strategy::within::spherical_winding<point_type>());
66 ::set_idsu(g1, g2, g3, bg::strategy::within::geographic_winding<point_type>());
67 }
68
69 template <typename G1, typename G2, typename G3>
set_idsu_ss(G1 const & g1,G2 const & g2,G3 & g3)70 inline void set_idsu_ss(G1 const& g1, G2 const& g2, G3 & g3)
71 {
72 ::set_idsu(g1, g2, g3, bg::strategy::intersection::cartesian_segments<>());
73 ::set_idsu(g1, g2, g3, bg::strategy::intersection::spherical_segments<>());
74 ::set_idsu(g1, g2, g3, bg::strategy::intersection::geographic_segments<>());
75 }
76
77 template <typename G1, typename G2, typename G3>
set_ids_pp(G1 const & g1,G2 const & g2,G3 & g3)78 inline void set_ids_pp(G1 const& g1, G2 const& g2, G3 & g3)
79 {
80 ::set_ids(g1, g2, g3, bg::strategy::within::cartesian_point_point());
81 ::set_ids(g1, g2, g3, bg::strategy::within::spherical_point_point());
82 }
83
84 template <typename G1, typename G2, typename G3>
set_ids_ps(G1 const & g1,G2 const & g2,G3 & g3)85 inline void set_ids_ps(G1 const& g1, G2 const& g2, G3 & g3)
86 {
87 typedef typename bg::point_type<G1>::type point_type;
88 ::set_ids(g1, g2, g3, bg::strategy::within::cartesian_winding<point_type>());
89 ::set_ids(g1, g2, g3, bg::strategy::within::spherical_winding<point_type>());
90 ::set_ids(g1, g2, g3, bg::strategy::within::geographic_winding<point_type>());
91 }
92
93 template <typename G1, typename G2, typename G3>
set_ids_ss(G1 const & g1,G2 const & g2,G3 & g3)94 inline void set_ids_ss(G1 const& g1, G2 const& g2, G3 & g3)
95 {
96 ::set_ids(g1, g2, g3, bg::strategy::intersection::cartesian_segments<>());
97 ::set_ids(g1, g2, g3, bg::strategy::intersection::spherical_segments<>());
98 ::set_ids(g1, g2, g3, bg::strategy::intersection::geographic_segments<>());
99 }
100
101 template <typename G1, typename G2, typename G3>
set_id_pp(G1 const & g1,G2 const & g2,G3 & g3)102 inline void set_id_pp(G1 const& g1, G2 const& g2, G3 & g3)
103 {
104 ::set_id(g1, g2, g3, bg::strategy::within::cartesian_point_point());
105 ::set_id(g1, g2, g3, bg::strategy::within::spherical_point_point());
106 }
107
108 template <typename G1, typename G2, typename G3>
set_id_ps(G1 const & g1,G2 const & g2,G3 & g3)109 inline void set_id_ps(G1 const& g1, G2 const& g2, G3 & g3)
110 {
111 typedef typename bg::point_type<G1>::type point_type;
112 ::set_id(g1, g2, g3, bg::strategy::within::cartesian_winding<point_type>());
113 ::set_id(g1, g2, g3, bg::strategy::within::spherical_winding<point_type>());
114 ::set_id(g1, g2, g3, bg::strategy::within::geographic_winding<point_type>());
115 }
116
117 template <typename G1, typename G2, typename G3>
set_id_ss(G1 const & g1,G2 const & g2,G3 & g3)118 inline void set_id_ss(G1 const& g1, G2 const& g2, G3 & g3)
119 {
120 ::set_id(g1, g2, g3, bg::strategy::intersection::cartesian_segments<>());
121 ::set_id(g1, g2, g3, bg::strategy::intersection::spherical_segments<>());
122 ::set_id(g1, g2, g3, bg::strategy::intersection::geographic_segments<>());
123 }
124
125 template <typename G1, typename G2, typename G3>
set_i_pp(G1 const & g1,G2 const & g2,G3 & g3)126 inline void set_i_pp(G1 const& g1, G2 const& g2, G3 & g3)
127 {
128 ::set_i(g1, g2, g3, bg::strategy::within::cartesian_point_point());
129 ::set_i(g1, g2, g3, bg::strategy::within::spherical_point_point());
130 }
131
132 template <typename G1, typename G2, typename G3>
set_i_ps(G1 const & g1,G2 const & g2,G3 & g3)133 inline void set_i_ps(G1 const& g1, G2 const& g2, G3 & g3)
134 {
135 typedef typename bg::point_type<G1>::type point_type;
136 ::set_i(g1, g2, g3, bg::strategy::within::cartesian_winding<point_type>());
137 ::set_i(g1, g2, g3, bg::strategy::within::spherical_winding<point_type>());
138 ::set_i(g1, g2, g3, bg::strategy::within::geographic_winding<point_type>());
139 }
140
141 template <typename G1, typename G2, typename G3>
set_i_ss(G1 const & g1,G2 const & g2,G3 & g3)142 inline void set_i_ss(G1 const& g1, G2 const& g2, G3 & g3)
143 {
144 ::set_i(g1, g2, g3, bg::strategy::intersection::cartesian_segments<>());
145 ::set_i(g1, g2, g3, bg::strategy::intersection::spherical_segments<>());
146 ::set_i(g1, g2, g3, bg::strategy::intersection::geographic_segments<>());
147 }
148
149 template <typename G1, typename G2, typename G3>
set_d_pp(G1 const & g1,G2 const & g2,G3 & g3)150 inline void set_d_pp(G1 const& g1, G2 const& g2, G3 & g3)
151 {
152 ::set_d(g1, g2, g3, bg::strategy::within::cartesian_point_point());
153 ::set_d(g1, g2, g3, bg::strategy::within::spherical_point_point());
154 }
155
156 template <typename G1, typename G2, typename G3>
set_d_ps(G1 const & g1,G2 const & g2,G3 & g3)157 inline void set_d_ps(G1 const& g1, G2 const& g2, G3 & g3)
158 {
159 typedef typename bg::point_type<G1>::type point_type;
160 ::set_d(g1, g2, g3, bg::strategy::within::cartesian_winding<point_type>());
161 ::set_d(g1, g2, g3, bg::strategy::within::spherical_winding<point_type>());
162 ::set_d(g1, g2, g3, bg::strategy::within::geographic_winding<point_type>());
163 }
164
165 template <typename G1, typename G2, typename G3>
set_d_ss(G1 const & g1,G2 const & g2,G3 & g3)166 inline void set_d_ss(G1 const& g1, G2 const& g2, G3 & g3)
167 {
168 ::set_d(g1, g2, g3, bg::strategy::intersection::cartesian_segments<>());
169 ::set_d(g1, g2, g3, bg::strategy::intersection::spherical_segments<>());
170 ::set_d(g1, g2, g3, bg::strategy::intersection::geographic_segments<>());
171 }
172
test_main(int,char * [])173 int test_main(int, char*[])
174 {
175 geom g;
176
177 // P/P->P
178 ::set_idsu_pp(g.pt, g.pt, g.mpt);
179 ::set_idsu_pp(g.pt, g.mpt, g.mpt);
180 ::set_idsu_pp(g.mpt, g.mpt, g.mpt);
181
182 // P/L->P
183 ::set_id_ps(g.pt, g.s, g.mpt);
184 ::set_id_ps(g.pt, g.ls, g.mpt);
185 ::set_id_ps(g.pt, g.mls, g.mpt);
186 ::set_id_ps(g.mpt, g.s, g.mpt);
187 ::set_id_ps(g.mpt, g.ls, g.mpt);
188 ::set_id_ps(g.mpt, g.mls, g.mpt);
189
190 // P/A->P
191 // no intersection nor difference
192 //::set_id_ps(g.pt, g.r, g.mpt);
193 //::set_id_ps(g.pt, g.po, g.mpt);
194 //::set_id_ps(g.pt, g.mpo, g.mpt);
195 //::set_id_ps(g.mpt, g.r, g.mpt);
196 //::set_id_ps(g.mpt, g.po, g.mpt);
197 //::set_id_ps(g.mpt, g.mpo, g.mpt);
198
199 // L/L->P
200 ::set_ids_ss(g.s, g.s, g.mpt);
201 //::set_i_ss(g.s, g.ls, g.mpt); // no intersection nor difference
202 //::set_i_ss(g.s, g.mls, g.mpt); // no intersection nor difference
203 //::set_i_ss(g.ls, g.s, g.mpt); // no intersection nor difference
204 ::set_ids_ss(g.ls, g.ls, g.mpt);
205 ::set_i_ss(g.ls, g.mls, g.mpt); // no difference nor sym_difference
206 //::set_i_ss(g.mls, g.s, g.mpt); // no intersection nor difference
207 ::set_i_ss(g.mls, g.ls, g.mpt); // no difference nor sym_difference
208 ::set_ids_ss(g.mls, g.mls, g.mpt);
209
210 // L/L->L
211 //::set_ids_ss(g.s, g.s, g.mls); // union not implemented, missing specialization
212 //::set_idsu_ss(g.s, g.ls, g.mls); // missing specialization
213 //::set_idsu_ss(g.s, g.mls, g.mls); // missing specialization
214 //::set_idsu_ss(g.ls, g.s, g.mls); // missing specialization
215 ::set_idsu_ss(g.ls, g.ls, g.mls);
216 ::set_idsu_ss(g.ls, g.mls, g.mls);
217 //::set_idsu_ss(g.mls, g.s, g.mls); // missing specialization
218 ::set_idsu_ss(g.mls, g.ls, g.mls);
219 ::set_idsu_ss(g.mls, g.mls, g.mls);
220
221 // S/B->L ?
222
223 // L/B->L ?
224
225 // L/A->P
226 //::set_ids_ss(g.s, g.r, g.mpt); // no intersection
227 //::set_ids_ss(g.s, g.po, g.mpt); // no intersection
228 //::set_ids_ss(g.s, g.mpo, g.mpt); // no intersection
229 ::set_ids_ss(g.ls, g.r, g.mpt);
230 ::set_ids_ss(g.ls, g.po, g.mpt);
231 ::set_ids_ss(g.ls, g.mpo, g.mpt);
232 ::set_ids_ss(g.mls, g.r, g.mpt);
233 ::set_ids_ss(g.mls, g.po, g.mpt);
234 ::set_ids_ss(g.mls, g.mpo, g.mpt);
235
236 // L/A->L
237 //::set_id_ss(g.s, g.r, g.mls); // no intersection
238 //::set_id_ss(g.s, g.po, g.mls); // no intersection
239 //::set_id_ss(g.s, g.mpo, g.mls); // no intersection
240 ::set_id_ss(g.ls, g.r, g.mls);
241 ::set_id_ss(g.ls, g.po, g.mls);
242 ::set_id_ss(g.ls, g.mpo, g.mls);
243 ::set_id_ss(g.mls, g.r, g.mls);
244 ::set_id_ss(g.mls, g.po, g.mls);
245 ::set_id_ss(g.mls, g.mpo, g.mls);
246
247 // A/A->P
248 ::set_i_ss(g.r, g.r, g.mpt);
249 ::set_i_ss(g.r, g.po, g.mpt);
250 ::set_i_ss(g.r, g.mpo, g.mpt);
251 ::set_i_ss(g.po, g.r, g.mpt);
252 ::set_i_ss(g.po, g.po, g.mpt);
253 ::set_i_ss(g.po, g.mpo, g.mpt);
254 ::set_i_ss(g.mpo, g.r, g.mpt);
255 ::set_i_ss(g.mpo, g.po, g.mpt);
256 ::set_i_ss(g.mpo, g.mpo, g.mpt);
257
258 // A/A->L
259 ::set_i_ss(g.r, g.r, g.mls);
260 ::set_i_ss(g.r, g.po, g.mls);
261 ::set_i_ss(g.r, g.mpo, g.mls);
262 ::set_i_ss(g.po, g.r, g.mls);
263 ::set_i_ss(g.po, g.po, g.mls);
264 ::set_i_ss(g.po, g.mpo, g.mls);
265 ::set_i_ss(g.mpo, g.r, g.mls);
266 ::set_i_ss(g.mpo, g.po, g.mls);
267 ::set_i_ss(g.mpo, g.mpo, g.mls);
268
269 // A/A->A
270 ::set_idsu_ss(g.r, g.r, g.mpo);
271 ::set_idsu_ss(g.r, g.po, g.mpo);
272 ::set_idsu_ss(g.r, g.mpo, g.mpo);
273 ::set_idsu_ss(g.po, g.r, g.mpo);
274 ::set_idsu_ss(g.po, g.po, g.mpo);
275 ::set_idsu_ss(g.po, g.mpo, g.mpo);
276 ::set_idsu_ss(g.mpo, g.r, g.mpo);
277 ::set_idsu_ss(g.mpo, g.po, g.mpo);
278 ::set_idsu_ss(g.mpo, g.mpo, g.mpo);
279
280 return 0;
281 }
282