]> git.proxmox.com Git - ceph.git/blame - ceph/src/boost/libs/geometry/test/algorithms/relate/relate_areal_areal.cpp
projects / ceph.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
update ceph source to reef 18.1.2
[ceph.git] / ceph / src / boost / libs / geometry / test / algorithms / relate / relate_areal_areal.cpp
CommitLineData
7c673cae
FG		 1// Boost.Geometry (aka GGL, Generic Geometry Library)
2
3// Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands.
4
1e59de90
TL		 5// This file was modified by Oracle on 2013-2021.
6// Modifications copyright (c) 2013-2021 Oracle and/or its affiliates.
7c673cae
FG		 7// Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
8
9// Use, modification and distribution is subject to the Boost Software License,
10// Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
11// http://www.boost.org/LICENSE_1_0.txt)
12
13
14#include "test_relate.hpp"
15
7c673cae
FG		 16
17template <typename P>
18void test_polygon_polygon()
19{
1e59de90
TL		 20 using coord_t = typename bg::coordinate_type<P>::type;
21 using poly = bg::model::polygon<P>;
22 using ring = bg::model::ring<P>;
7c673cae
FG		 23
24 // touching
25 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
26 "POLYGON((10 0,10 10,20 10,20 0,10 0))",
27 "FF2F11212");
28 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
29 "POLYGON((0 -10,0 0,10 0,10 -10,0 -10))",
30 "FF2F11212");
31 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
32 "POLYGON((10 0,15 10,20 10,20 0,10 0))",
33 "FF2F01212");
34
35 // containing
36 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
37 "POLYGON((5 5,5 10,6 10,6 5,5 5))",
38 "212F11FF2");
39 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
40 "POLYGON((5 5,5 10,6 5,5 5))",
41 "212F01FF2");
42 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
43 "POLYGON((5 5,5 6,6 6,6 5,5 5))",
44 "212FF1FF2");
45
46 // fully containing
47 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
48 "POLYGON((5 5,5 9,6 9,6 5,5 5))",
49 "212FF1FF2");
50 // fully containing, with a hole
51 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(3 3,7 3,7 7,3 7,3 3))",
52 "POLYGON((1 1,1 9,9 9,9 1,1 1))",
53 "2121F12F2");
54 // fully containing, both with holes
55 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(3 3,7 3,7 7,3 7,3 3))",
56 "POLYGON((1 1,1 9,9 9,9 1,1 1),(2 2,8 2,8 8,2 8,2 2))",
57 "212FF1FF2");
58 // fully containing, both with holes
59 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(3 3,7 3,7 7,3 7,3 3))",
60 "POLYGON((1 1,1 9,9 9,9 1,1 1),(4 4,6 4,6 6,4 6,4 4))",
61 "2121F1212");
62
63 // overlapping
64 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
65 "POLYGON((5 0,5 10,20 10,20 0,5 0))",
66 "212111212");
67 test_geometry<ring, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
68 "POLYGON((5 0,5 10,20 10,20 0,5 0))",
69 "212111212");
70 test_geometry<ring, ring>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
71 "POLYGON((5 0,5 10,20 10,20 0,5 0))",
72 "212111212");
73 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,15 5,10 0,0 0))",
74 "POLYGON((10 0,5 5,10 10,20 10,20 0,10 0))",
75 "212101212");
76
77 // equal
78 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
79 "POLYGON((10 10,10 5,10 0,5 0,0 0,0 10,5 10,10 10))",
80 "2FFF1FFF2");
81 // hole-sized
82 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(5 5,6 5,6 6,5 6,5 5))",
83 "POLYGON((5 5,5 6,6 6,6 5,5 5))",
84 "FF2F112F2");
85
86 // disjoint
87 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
88 "POLYGON((0 20,0 30,10 30,10 20,0 20))",
89 "FF2FF1212");
90 // disjoint
91 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(3 3,7 3,7 7,3 7,3 3))",
92 "POLYGON((0 20,0 30,10 30,10 20,0 20))",
93 "FF2FF1212");
94
95 // equal non-simple / non-simple hole
96 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(5 5,10 5,5 6,5 5))",
97 "POLYGON((0 0,0 10,10 10,10 5,5 6,5 5,10 5,10 0,0 0))",
98 "2FFF1FFF2");
99
100 // within non-simple / simple
101 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 5,5 6,5 5,10 5,10 0,0 0))",
102 "POLYGON((0 0,5 5,10 5,10 0,0 0))",
103 "212F11FF2");
104 // within non-simple hole / simple
105 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(5 5,10 5,5 6,5 5))",
106 "POLYGON((0 0,5 5,10 5,10 0,0 0))",
107 "212F11FF2");
108
109
110 // not within non-simple / simple
111 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 5,5 6,5 5,10 5,10 0,0 0))",
112 "POLYGON((0 0,0 10,10 10,10 0,0 0))",
113 "2FF11F2F2");
114 // not within non-simple hole / simple
115 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(5 5,10 5,5 6,5 5))",
116 "POLYGON((0 0,0 10,10 10,10 0,0 0))",
117 "2FF11F2F2");
118 // not within simple hole / simple
119 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(5 5,9 5,5 6,5 5))",
120 "POLYGON((0 0,0 10,10 10,9 5,10 0,0 0))",
121 "2121112F2");
122
123 // within non-simple fake hole / simple
124 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 5,4 7,4 3,10 5,10 0,0 0))",
125 "POLYGON((0 0,0 10,10 10,10 0,0 0))",
126 "2FF11F2F2");
127 // within non-simple fake hole / non-simple fake hole
128 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 5,4 7,4 3,10 5,10 0,0 0))",
129 "POLYGON((0 0,0 10,10 10,10 5,5 6,5 4,10 5,10 0,0 0))",
130 "2FF11F212");
131 // within non-simple fake hole / non-simple hole
132 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 5,4 7,4 3,10 5,10 0,0 0))",
133 "POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,5 6,5 4,10 5))",
134 "2FF11F212");
135 // containing non-simple fake hole / non-simple hole
136 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 5,4 7,4 3,10 5,10 0,0 0))",
137 "POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,3 8,3 2,10 5))",
138 "212F11FF2");
139
140 // within non-simple hole / simple
141 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,4 7,4 3,10 5))",
142 "POLYGON((0 0,0 10,10 10,10 0,0 0))",
143 "2FF11F2F2");
144 // within non-simple hole / non-simple fake hole
145 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,4 7,4 3,10 5))",
146 "POLYGON((0 0,0 10,10 10,10 5,5 6,5 4,10 5,10 0,0 0))",
147 "2FF11F212");
148 // containing non-simple hole / non-simple fake hole
149 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,4 7,4 3,10 5))",
150 "POLYGON((0 0,0 10,10 10,10 5,3 8,3 2,10 5,10 0,0 0))",
151 "212F11FF2");
152 // equal non-simple hole / non-simple fake hole
153 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,4 7,4 3,10 5))",
154 "POLYGON((0 0,0 10,10 10,10 5,4 7,4 3,10 5,10 0,0 0))",
155 "2FFF1FFF2");
156 // within non-simple hole / non-simple hole
157 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,4 7,4 3,10 5))",
158 "POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,5 6,5 4,10 5))",
159 "2FF11F212");
160 // containing non-simple hole / non-simple hole
161 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,4 7,4 3,10 5))",
162 "POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,3 8,3 2,10 5))",
163 "212F11FF2");
164 // equal non-simple hole / non-simple hole
165 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,4 7,4 3,10 5))",
166 "POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,4 7,4 3,10 5))",
167 "2FFF1FFF2");
168
169 // intersecting non-simple hole / non-simple hole - touching holes
170 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(5 5,10 5,5 6,5 5))",
171 "POLYGON((0 0,0 10,10 10,10 0,0 0),(0 5,5 4,5 5,0 5))",
172 "21211F2F2");
173 // intersecting non-simple fake hole / non-simple hole - touching holes
174 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 5,5 6,5 5,10 5,10 0,0 0))",
175 "POLYGON((0 0,0 10,10 10,10 0,0 0),(0 5,5 4,5 5,0 5))",
176 "21211F2F2");
177 // intersecting non-simple fake hole / non-simple fake hole - touching holes
178 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 5,5 6,5 5,10 5,10 0,0 0))",
179 "POLYGON((0 0,0 5,5 4,5 5,0 5,0 10,10 10,10 0,0 0))",
180 "21211F2F2");
181
182 // intersecting simple - i/i
183 test_geometry<poly, poly>("POLYGON((0 0,0 10,4 10,6 8,5 5,6 2,4 0,0 0))",
184 "POLYGON((5 5,4 8,6 10,10 10,10 0,6 0,4 2,5 5))",
185 "212101212");
186 // intersecting non-simple hole / non-simple hole - i/i
187 test_geometry<poly, poly>("POLYGON((0 0,0 10,4 10,6 8,5 5,6 2,4 0,0 0),(5 5,2 6,2 4,5 5))",
188 "POLYGON((5 5,4 8,6 10,10 10,10 0,6 0,4 2,5 5),(5 5,8 4,8 6,5 5))",
189 "212101212");
190 // intersecting non-simple hole / simple - i/i
191 test_geometry<poly, poly>("POLYGON((0 0,0 10,4 10,6 8,5 5,6 2,4 0,0 0),(5 5,2 6,2 4,5 5))",
192 "POLYGON((5 5,4 8,6 10,10 10,10 0,6 0,4 2,5 5))",
193 "212101212");
194
195 // no turns - disjoint inside a hole
196 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(1 1,9 1,9 9,1 9,1 1))",
197 "POLYGON((3 3,3 7,7 7,7 3,3 3))",
198 "FF2FF1212");
199 // no turns - within
200 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(1 1,9 1,9 9,1 9,1 1))",
201 "POLYGON((-1 -1,-1 11,11 11,11 -1,-1 -1))",
202 "2FF1FF212");
203 // no-turns - intersects
204 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,8 2,8 8,2 8,2 2))",
205 "POLYGON((1 1,1 9,9 9,9 1,1 1))",
206 "2121F12F2");
207 // no-turns - intersects, hole in a hole
208 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,8 2,8 8,2 8,2 2))",
209 "POLYGON((1 1,1 9,9 9,9 1,1 1),(3 3,7 3,7 7,3 7,3 3))",
210 "2121F1212");
211
212 // no-turns ring - for exteriors
213 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,8 2,8 8,2 8,2 2))",
214 "POLYGON((1 1,1 9,9 9,9 1,1 1),(2 2,8 2,8 8,2 8,2 2))",
215 "212F11FF2");
216 // no-turns ring - for interiors
217 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(3 3,7 3,7 7,3 7,3 3))",
218 "POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,8 2,8 8,2 8,2 2))",
219 "212F11FF2");
220
221 {
222 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
223 "POLYGON((5 5,5 10,6 10,6 5,5 5))",
224 "212F11FF2");
225
226 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
227 "POLYGON((10 0,10 10,20 10,20 0,10 0))",
228 "FF2F11212");
229
230 namespace bgdr = bg::detail::relate;
231 poly p1, p2, p3;
232 bg::read_wkt("POLYGON((0 0,0 10,10 10,10 0,0 0))", p1);
233 bg::read_wkt("POLYGON((10 0,10 10,20 10,20 0,10 0))", p2);
234 bg::read_wkt("POLYGON((5 5,5 10,6 10,6 5,5 5))", p3);
235 BOOST_CHECK(bg::relate(p1, p2, bg::de9im::mask("FT*******")
236 || bg::de9im::mask("F**T*****")
237 || bg::de9im::mask("F***T****"))); // touches()
238 BOOST_CHECK(bg::relate(p1, p3, bg::de9im::mask("T*****FF*"))); // contains()
239 BOOST_CHECK(bg::relate(p2, p3, bg::de9im::mask("FF*FF****"))); // disjoint()
240
241 BOOST_CHECK(bg::relate(p1, p2, bg::de9im::static_mask<'F','T'>()
242 || bg::de9im::static_mask<'F','*','*','T'>()
243 || bg::de9im::static_mask<'F','*','*','*','T'>()));
244 }
245
246 // CCW
247 {
248 typedef bg::model::polygon<P, false> poly;
249 // within non-simple hole / simple
250 test_geometry<poly, poly>("POLYGON((0 0,10 0,10 10,0 10,0 0),(5 5,5 6,10 5,5 5))",
251 "POLYGON((0 0,10 0,10 5,5 5,0 0))",
252 "212F11FF2");
253 }
254 // OPEN
255 {
256 typedef bg::model::polygon<P, true, false> poly;
257 // within non-simple hole / simple
258 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0),(5 5,10 5,5 6))",
259 "POLYGON((0 0,5 5,10 5,10 0))",
260 "212F11FF2");
261 }
262 // CCW, OPEN
263 {
264 typedef bg::model::polygon<P, false, false> poly;
265 // within non-simple hole / simple
266 test_geometry<poly, poly>("POLYGON((0 0,10 0,10 10,0 10),(5 5,5 6,10 5))",
267 "POLYGON((0 0,10 0,10 5,5 5))",
268 "212F11FF2");
269 }
270
271 // https://svn.boost.org/trac/boost/ticket/10912
272 {
273 // external rings touches and G2's hole is inside G1
274 test_geometry<poly, poly>("POLYGON((0 0,0 5,5 5,5 0,0 0))",
275 "POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,4 2,4 4,2 4,2 2),(6 6,8 6,8 8,6 8,6 6))",
276 "21211F212");
277 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
278 "POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,2 4,4 4,4 2,2 2))",
279 "212F1FFF2");
280 // extended
281 // external rings touches and G1's hole is inside G2
282 test_geometry<poly, poly>("POLYGON((0 0,0 9,9 9,9 0,0 0),(2 2,7 2,7 7,2 7,2 2))",
283 "POLYGON((0 0,0 10,10 10,10 0,0 0))",
284 "2FF11F212");
285 // external rings touches, holes doesn't
286 test_geometry<poly, poly>("POLYGON((0 0,0 9,9 9,9 0,0 0),(2 2,7 2,7 7,2 7,2 2))",
287 "POLYGON((0 0,0 10,10 10,10 0,0 0),(1 1,8 1,8 8,1 8,1 1))",
288 "212111212");
289 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,7 2,7 7,2 7,2 2))",
290 "POLYGON((0 0,0 10,10 10,10 0,0 0),(1 1,8 1,8 8,1 8,1 1))",
291 "212F11FF2");
292 // holes touches, external rings doesn't
293 test_geometry<poly, poly>("POLYGON((1 1,1 9,9 9,9 1,1 1),(2 2,8 2,8 8,2 8,2 2))",
294 "POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,7 2,7 7,2 7,2 2))",
295 "2FF11F212");
296 test_geometry<poly, poly>("POLYGON((1 1,1 9,9 9,9 1,1 1),(2 2,7 2,7 7,2 7,2 2))",
297 "POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,8 2,8 8,2 8,2 2))",
298 "212111212");
299 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,8 2,8 8,2 8,2 2))",
300 "POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,7 2,7 7,2 7,2 2))",
301 "2FF11F212");
302
303 test_geometry<poly, poly>("POLYGON((3 3,3 9,9 9,9 3,3 3))",
304 "POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,4 2,4 4,2 4,2 2))",
305 "212101212");
306 }
307
1e59de90 308 if ( BOOST_GEOMETRY_CONDITION((std::is_same<coord_t, double>::value)) )
7c673cae
FG		 309 {
310 // original - assertion for CCW
311 //"POLYGON((-0.593220338983050821113352 -8.05084745762711939676137,1.14285714285714279370154 -4,1.50731707317073171381594 1.10243902439024443751237,1.73758865248226967992196 1.37588652482269591104114,1.21739130434782616418943 -3.82608695652173924628414,2 -2,2 1.68750000000000044408921,2.35384615384615436539661 2.10769230769230775379697,2 2.16666666666666651863693,2 4,1.81967213114754100544701 2.1967213114754100544701,1.5882352941176469673934 2.2352941176470588757752,1.8148148148148146585612 5.4074074074074074403029,-0.538461538461538546940233 4.23076923076923083755219,-1.76510067114094004736558 2.89261744966443012927471,-1.64864864864864868465588 2.7567567567567570208098,-1.83962264150943455298659 2.81132075471698161805989,-1.84337349397590433142113 2.80722891566265086993326,-2.14285714285714279370154 2.85714285714285720629846,-2.11111111111111116045436 2.88888888888888883954564,-2.87234042553191448732264 3.10638297872340407579372,-2.91803278688524558859285 3.4262295081967208965068,-3.1733333333333324510761 3.26666666666666660745477,-2.99999999999999822364316 3.14285714285714234961233,-3.25490196078431326398572 3.21568627450980359938626,-3.47368421052631504153396 3.07894736842105265495206,-7.32000000000000028421709 3.72000000000000019539925,-7.54716981132075481752963 3.62264150943396234794136,-7.75 3.79166666666666651863693,-7.79999999999999982236432 3.79999999999999982236432,-7.59999999999999964472863 3.60000000000000008881784,-8.8556701030927822415606 3.06185567010309300783888,-8.82945736434108674473009 2.8914728682170549589614,-7.73333333333333339254523 2.193939393939393855959,-8 2,-5.94736842105263185942476 -1.42105263157894645686952,-5.32558139534883689947264 -0.488372093023255016142059,-5.85714285714285765038767 1.00000000000000066613381,-4.78723404255319184841255 0.319148936170212838003835,-5.32558139534883689947264 -0.488372093023255016142059,-4.74019607843137258385013 -2.12745098039215774221589,-3.17647058823529437887601 -0.705882352941176627325603,-2.93103448275862055183438 -0.862068965517241436735674,-3 -1,-4.57894736842105309904127 -2.57894736842105265495206,-4.47887323943661996850096 -2.85915492957746497637572,-7.58620689655172419918472 -5.18965517241379359347775,-7.52525252525252508206677 -5.5858585858585865224768,-4.18644067796610119813749 -3.67796610169491522412955,-3.44041450777202051369841 -5.76683937823834202873741,-3.73611111111111116045436 -6.56944444444444464181743,-2.8823529411764705621124 -7.7647058823529411242248,-2.88235294117647100620161 -7.7647058823529411242248,-0.593220338983050821113352 -8.05084745762711939676137),(1.66666666666666696272614 1.59999999999999875655021,1.43749999999999911182158 1.8750000000000002220446,0.0869565217391310429917439 2.26086956521739113057379,0.466666666666667118157363 2.60606060606060552231611,1.04878048780487764801705 2.34146341463414664474385,1.43749999999999911182158 1.8750000000000002220446,1.56756756756756754356275 1.83783783783783771781373,1.66666666666666696272614 1.59999999999999875655021))"
312 //"POLYGON((-2.33333333333333303727386 -8.66666666666666607454772,-2.26315789473684203514381 -8.63157894736842123961651,-2.88235294117647100620161 -7.7647058823529411242248,-2.8823529411764705621124 -7.7647058823529411242248,-4.11949685534591125701809 -7.61006289308176064878353,-4.32530120481927671249878 -8.16867469879518104391991,-2.33333333333333303727386 -8.66666666666666607454772))"
313
314 // assertion failure in 1.57
315 test_geometry<poly, poly>("POLYGON((-0.59322033898305082 -8.0508474576271194,-2.882352941176471 -7.7647058823529411,-2.8823529411764706 -7.7647058823529411,-3.7361111111111112 -6.5694444444444446,-3.4404145077720205 -5.766839378238342,-4.1864406779661012 -3.6779661016949152,-7.5252525252525251 -5.5858585858585865,-7.5862068965517242 -5.1896551724137936,-4.47887323943662 -2.859154929577465,-4.5789473684210531 -2.5789473684210527,-3 -1,-2.9310344827586206 -0.86206896551724144,-3.1764705882352944 -0.70588235294117663,-4.7401960784313726 -2.1274509803921577,-5.3255813953488369 -0.48837209302325502,-4.7872340425531918 0.31914893617021284,-5.8571428571428577 1.0000000000000007,-5.3255813953488369 -0.48837209302325502,-5.9473684210526319 -1.4210526315789465,-8 2,-7.7333333333333334 2.1939393939393939,-8.8294573643410867 2.891472868217055,-8.8556701030927822 3.061855670103093,-7.5999999999999996 3.6000000000000001,-7.7999999999999998 3.7999999999999998,-7.75 3.7916666666666665,-7.5471698113207548 3.6226415094339623,-7.3200000000000003 3.7200000000000002,-3.473684210526315 3.0789473684210527,-3.2549019607843133 3.2156862745098036,-2.9999999999999982 3.1428571428571423,-3.1733333333333325 3.2666666666666666,-2.9180327868852456 3.4262295081967209,-2.8723404255319145 3.1063829787234041,-2.1111111111111112 2.8888888888888888,-2.1428571428571428 2.8571428571428572,-1.8433734939759043 2.8072289156626509,-1.8396226415094346 2.8113207547169816,-1.6486486486486487 2.756756756756757,-1.76510067114094 2.8926174496644301,-0.53846153846153855 4.2307692307692308,1.8148148148148147 5.4074074074074074,1.588235294117647 2.2352941176470589,1.819672131147541 2.1967213114754101,2 4,2 2.1666666666666665,2.3538461538461544 2.1076923076923078,2 1.6875000000000004,2 -2,1.2173913043478262 -3.8260869565217392,1.7375886524822697 1.3758865248226959,1.5073170731707317 1.1024390243902444,1.1428571428571428 -4,-0.59322033898305082 -8.0508474576271194),(1.666666666666667 1.5999999999999988,1.5675675675675675 1.8378378378378377,1.4374999999999991 1.8750000000000002,1.0487804878048776 2.3414634146341466,0.46666666666666712 2.6060606060606055,0.086956521739131043 2.2608695652173911,1.4374999999999991 1.8750000000000002,1.666666666666667 1.5999999999999988))",
316 "POLYGON((-2.333333333333333 -8.6666666666666661,-4.3253012048192767 -8.168674698795181,-4.1194968553459113 -7.6100628930817606,-2.8823529411764706 -7.7647058823529411,-2.882352941176471 -7.7647058823529411,-2.263157894736842 -8.6315789473684212,-2.333333333333333 -8.6666666666666661))",
317 "*********");
318 // simpler case
319 test_geometry<poly, poly>("POLYGON((0 0, -0.59322033898305082 -8.0508474576271194, -2.8823529411764710 -7.7647058823529411, -3.7361111111111112 -6.5694444444444446, 0 0))",
320 "POLYGON((-10 -10, -4.1194968553459113 -7.6100628930817606, -2.8823529411764706 -7.7647058823529411, -2.2631578947368420 -8.6315789473684212, -10 -10))",
321 "FF2F01212");
322 // sanity check
323 test_geometry<poly, poly>("POLYGON((0 0, -0.59322033898305082 -8.0508474576271194, -2.8823529411764710 -7.7647058823529411, -3.7361111111111112 -6.5694444444444446, 0 0))",
324 "POLYGON((-10 -10, -4.1194968553459113 -7.6100628930817606, -2.8823529411764710 -7.7647058823529411, -2.2631578947368420 -8.6315789473684212, -10 -10))",
325 "FF2F01212");
326
327 // compatibility check
328 /*test_geometry<poly, poly>("POLYGON((0 0,0 1,1 1,1 0,0 0))",
329 "POLYGON((1 1,1 2,2 2,2 1,1 1))",
330 "****0****");
331 test_geometry<P, P>("POINT(0.9999999999999998 0.9999999999999998)",
332 "POINT(1 1)",
333 "0********");
334 test_geometry<P, P>("POINT(0.9999999999999995 0.9999999999999995)",
335 "POINT(1 1)",
336 "F********");
337 test_geometry<poly, poly>("POLYGON((0 0,0 1,0.9999999999999998 0.9999999999999998,1 0,0 0))",
338 "POLYGON((1 1,1 2,2 2,2 1,1 1))",
339 "****0****");
340 test_geometry<poly, poly>("POLYGON((0 0,0 1,0.9999999999999995 0.9999999999999995,1 0,0 0))",
341 "POLYGON((1 1,1 2,2 2,2 1,1 1))",
342 "****F****");*/
343 }
344
345 // mysql 21872795 (overlaps)
346 test_geometry<poly, poly>("POLYGON((2 2,2 8,8 8,8 2,2 2))",
347 "POLYGON((0 0,0 10,10 10,10 0,0 0),(8 8,4 6,4 4,8 8))",
348 "21210F212");
349 test_geometry<poly, poly>("POLYGON((2 2,2 8,8 8,8 2,2 2))",
350 "POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,4 4,4 6,2 2))",
351 "21210F212");
352 // mysql 21873343 (touches)
353 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0), (0 8, 8 5, 8 8, 0 8))",
354 "POLYGON((0 8,-8 5,-8 8,0 8))",
355 "FF2F01212");
356 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0), (0 6, 6 3, 6 6, 0 6))",
357 "POLYGON((0 6,-6 3,-6 6,0 6))",
358 "FF2F01212");
359 // similar but touching from the inside of the hole
360 test_geometry<poly, poly>("POLYGON((0 0,0 10,10 10,10 0,0 0), (0 8, 8 5, 8 8, 0 8))",
361 "POLYGON((0 8,7 7, 7 6,0 8))",
362 "FF2F01212");
363}
364
365template <typename P>
366void test_polygon_multi_polygon()
367{
368 typedef bg::model::polygon<P> poly;
369 typedef bg::model::ring<P> ring;
370 typedef bg::model::multi_polygon<poly> mpoly;
371
372 test_geometry<poly, mpoly>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
373 "MULTIPOLYGON(((5 5,5 10,6 10,6 5,5 5)),((0 20,0 30,10 30,10 20,0 20)))",
374 "212F11212");
375 test_geometry<ring, mpoly>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
376 "MULTIPOLYGON(((5 5,5 10,6 10,6 5,5 5)),((0 20,0 30,10 30,10 20,0 20)))",
377 "212F11212");
378
379 test_geometry<mpoly, poly>("MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)),((0 0,0 -10,-10 -10,-10 0,0 0)))",
380 "POLYGON((0 0,0 5,5 5,5 0,0 0))",
381 "212F11FF2");
382 test_geometry<poly, mpoly>("POLYGON((0 0,0 5,5 5,5 0,0 0))",
383 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)),((0 0,0 -10,-10 -10,-10 0,0 0)))",
384 "2FF11F212");
385 test_geometry<poly, mpoly>("POLYGON((0 0,0 5,5 5,5 0,0 0))",
386 "MULTIPOLYGON(((0 0,0 -10,-10 -10,-10 0,0 0)),((0 0,0 10,10 10,10 0,0 0)))",
387 "2FF11F212");
388 test_geometry<poly, mpoly>("POLYGON((0 0,0 10,10 10,10 0,0 0))",
389 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)),((0 0,0 -10,-10 -10,-10 0,0 0)))",
390 "2FFF1F212");
391}
392
393template <typename P>
394void test_multi_polygon_multi_polygon()
395{
396 typedef bg::model::polygon<P> poly;
397 typedef bg::model::multi_polygon<poly> mpoly;
398
399 test_geometry<mpoly, mpoly>("MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)))",
400 "MULTIPOLYGON(((5 5,5 10,6 10,6 5,5 5)),((0 20,0 30,10 30,10 20,0 20)))",
401 "212F11212");
402 test_geometry<mpoly, mpoly>("MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)),((0 20,0 30,10 30,10 20,0 20)))",
403 "MULTIPOLYGON(((5 5,5 10,6 10,6 5,5 5)))",
404 "212F11FF2");
405
406 test_geometry<mpoly, mpoly>("MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)))",
407 "MULTIPOLYGON(((5 5,5 6,6 6,6 5,5 5)),((0 20,0 30,10 30,10 20,0 20)))",
408 "212FF1212");
409 test_geometry<mpoly, mpoly>("MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)),((0 20,0 30,10 30,10 20,0 20)))",
410 "MULTIPOLYGON(((5 5,5 6,6 6,6 5,5 5)))",
411 "212FF1FF2");
412}
413
414template <typename P>
415void test_all()
416{
417 test_polygon_polygon<P>();
418 test_polygon_multi_polygon<P>();
419 test_multi_polygon_multi_polygon<P>();
420}
421
422int test_main( int , char* [] )
423{
424 test_all<bg::model::d2::point_xy<int> >();
425 test_all<bg::model::d2::point_xy<double> >();
426
7c673cae
FG		 427 return 0;
428}
Ceph