1 // Boost.Geometry (aka GGL, Generic Geometry Library)
2
3 // Copyright (c) 2007-2015 Barend Gehrels, Amsterdam, the Netherlands.
4
5 // This file was modified by Oracle on 2013, 2014, 2015.
6 // Modifications copyright (c) 2013-2015 Oracle and/or its affiliates.
7
8 // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
9
10 // Use, modification and distribution is subject to the Boost Software License,
11 // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
12 // http://www.boost.org/LICENSE_1_0.txt)
13
14 #include "test_relate.hpp"
15
16 //TEST
17 //#include <to_svg.hpp>
18
19 template <typename P>
test_linestring_linestring()20 void test_linestring_linestring()
21 {
22 typedef bg::model::linestring<P> ls;
23
24 test_geometry<ls, ls>("LINESTRING(0 0, 2 2, 3 2)", "LINESTRING(0 0, 2 2, 3 2)", "1FFF0FFF2");
25 test_geometry<ls, ls>("LINESTRING(0 0,3 2)", "LINESTRING(0 0, 2 2, 3 2)", "FF1F0F1F2");
26 test_geometry<ls, ls>("LINESTRING(1 0,2 2,2 3)", "LINESTRING(0 0, 2 2, 3 2)", "0F1FF0102");
27
28 test_geometry<ls, ls>("LINESTRING(0 0,5 0)", "LINESTRING(0 0,1 1,2 0,2 -1)", "0F1F00102");
29
30 test_geometry<ls, ls>("LINESTRING(0 0, 1 1, 2 2, 3 2)", "LINESTRING(0 0, 2 2, 3 2)", "1FFF0FFF2");
31 test_geometry<ls, ls>("LINESTRING(3 2, 2 2, 1 1, 0 0)", "LINESTRING(0 0, 2 2, 3 2)", "1FFF0FFF2");
32 test_geometry<ls, ls>("LINESTRING(0 0, 1 1, 2 2, 3 2)", "LINESTRING(3 2, 2 2, 0 0)", "1FFF0FFF2");
33 test_geometry<ls, ls>("LINESTRING(3 2, 2 2, 1 1, 0 0)", "LINESTRING(3 2, 2 2, 0 0)", "1FFF0FFF2");
34
35 test_geometry<ls, ls>("LINESTRING(3 1, 2 2, 1 1, 0 0)", "LINESTRING(0 0, 2 2, 3 2)", "1F1F00102");
36 test_geometry<ls, ls>("LINESTRING(3 3, 2 2, 1 1, 0 0)", "LINESTRING(0 0, 2 2, 3 2)", "1F1F00102");
37
38 test_geometry<ls, ls>("LINESTRING(0 0, 1 1, 2 2, 2 3)", "LINESTRING(0 0, 2 2, 2 3)", "1FFF0FFF2");
39 test_geometry<ls, ls>("LINESTRING(2 3, 2 2, 1 1, 0 0)", "LINESTRING(0 0, 2 2, 2 3)", "1FFF0FFF2");
40 test_geometry<ls, ls>("LINESTRING(0 0, 1 1, 2 2, 2 3)", "LINESTRING(2 3, 2 2, 0 0)", "1FFF0FFF2");
41 test_geometry<ls, ls>("LINESTRING(2 3, 2 2, 1 1, 0 0)", "LINESTRING(2 3, 2 2, 0 0)", "1FFF0FFF2");
42
43 test_geometry<ls, ls>("LINESTRING(1 1, 2 2, 3 2)", "LINESTRING(0 0, 2 2, 4 2)", "1FF0FF102");
44
45 test_geometry<ls, ls>("LINESTRING(3 2, 2 2, 1 1)", "LINESTRING(0 0, 2 2, 4 2)", "1FF0FF102");
46 test_geometry<ls, ls>("LINESTRING(1 1, 2 2, 3 2)", "LINESTRING(4 2, 2 2, 0 0)", "1FF0FF102");
47 test_geometry<ls, ls>("LINESTRING(3 2, 2 2, 1 1)", "LINESTRING(4 2, 2 2, 0 0)", "1FF0FF102");
48
49 // test_geometry<ls, ls>("LINESTRING(1 1, 2 2, 2 2)", "LINESTRING(0 0, 2 2, 4 2)", true);
50
51 // test_geometry<ls, ls>("LINESTRING(1 1, 2 2, 3 3)", "LINESTRING(0 0, 2 2, 4 2)", false);
52 // test_geometry<ls, ls>("LINESTRING(1 1, 2 2, 3 2, 3 3)", "LINESTRING(0 0, 2 2, 4 2)", false);
53 // test_geometry<ls, ls>("LINESTRING(1 1, 2 2, 3 1)", "LINESTRING(0 0, 2 2, 4 2)", false);
54 // test_geometry<ls, ls>("LINESTRING(1 1, 2 2, 3 2, 3 1)", "LINESTRING(0 0, 2 2, 4 2)", false);
55
56 // test_geometry<ls, ls>("LINESTRING(0 1, 1 1, 2 2, 3 2)", "LINESTRING(0 0, 2 2, 4 2)", false);
57 // test_geometry<ls, ls>("LINESTRING(0 1, 0 0, 2 2, 3 2)", "LINESTRING(0 0, 2 2, 4 2)", false);
58 // test_geometry<ls, ls>("LINESTRING(1 0, 1 1, 2 2, 3 2)", "LINESTRING(0 0, 2 2, 4 2)", false);
59 // test_geometry<ls, ls>("LINESTRING(1 0, 0 0, 2 2, 3 2)", "LINESTRING(0 0, 2 2, 4 2)", false);
60
61 // test_geometry<ls, ls>("LINESTRING(0 0)", "LINESTRING(0 0)", false);
62 // test_geometry<ls, ls>("LINESTRING(1 1)", "LINESTRING(0 0, 2 2)", true);
63 // test_geometry<ls, ls>("LINESTRING(0 0)", "LINESTRING(0 0, 2 2)", false);
64 // test_geometry<ls, ls>("LINESTRING(0 0, 1 1)", "LINESTRING(0 0)", false);
65
66 // test_geometry<ls, ls>("LINESTRING(0 0,5 0,3 0,6 0)", "LINESTRING(0 0,6 0)", true);
67 // test_geometry<ls, ls>("LINESTRING(0 0,2 2,3 3,1 1)", "LINESTRING(0 0,3 3,6 3)", true);
68
69 // SPIKES!
70 test_geometry<ls, ls>("LINESTRING(0 0,2 2,3 3,1 1,5 3)", "LINESTRING(0 0,3 3,6 3)", "1F100F102");
71 test_geometry<ls, ls>("LINESTRING(5 3,1 1,3 3,2 2,0 0)", "LINESTRING(0 0,3 3,6 3)", "1F100F102");
72 test_geometry<ls, ls>("LINESTRING(0 0,2 2,3 3,1 1,5 3)", "LINESTRING(6 3,3 3,0 0)", "1F100F102");
73 test_geometry<ls, ls>("LINESTRING(5 3,1 1,3 3,2 2,0 0)", "LINESTRING(6 3,3 3,0 0)", "1F100F102");
74
75 test_geometry<ls, ls>("LINESTRING(6 3,3 3,0 0)", "LINESTRING(0 0,2 2,3 3,1 1,5 3)", "101F001F2");
76
77 test_geometry<ls, ls>("LINESTRING(0 0,10 0)", "LINESTRING(1 0,9 0,2 0)", "101FF0FF2");
78 test_geometry<ls, ls>("LINESTRING(0 0,2 2,3 3,1 1)", "LINESTRING(0 0,3 3,6 3)", "1FF00F102");
79 // TODO: REWRITE MATRICES
80 // BEGIN
81 /*test_geometry<ls, ls>("LINESTRING(0 0,2 2,3 3,1 1)", "LINESTRING(0 0,4 4,6 3)", "1FF00F102");
82
83 test_geometry<ls, ls>("LINESTRING(0 0,2 0,1 0)", "LINESTRING(0 1,0 0,2 0)", "1FF00F102");
84 test_geometry<ls, ls>("LINESTRING(2 0,0 0,1 0)", "LINESTRING(0 1,0 0,2 0)", "1FF00F102");
85
86 test_geometry<ls, ls>("LINESTRING(0 0,3 3,1 1)", "LINESTRING(3 0,3 3,3 1)", "0F1FF0102");
87 test_geometry<ls, ls>("LINESTRING(0 0,3 3,1 1)", "LINESTRING(2 0,2 3,2 1)", "0F1FF0102");
88 test_geometry<ls, ls>("LINESTRING(0 0,3 3,1 1)", "LINESTRING(2 0,2 2,2 1)", "0F1FF0102");
89
90 test_geometry<ls, ls>("LINESTRING(0 0,2 2,3 3,4 4)", "LINESTRING(0 0,1 1,4 4)", "1FFF0FFF2");*/
91 // END
92
93 test_geometry<ls, ls>("LINESTRING(0 0,2 2,3 3,4 4)", "LINESTRING(0 0,1 1,4 4)", "1FFF0FFF2");
94
95 // loop i/i i/i u/u u/u
96 test_geometry<ls, ls>("LINESTRING(0 0,10 0)",
97 "LINESTRING(1 1,1 0,6 0,6 1,4 1,4 0,9 0,9 1)", "1F1FF0102");
98
99 // self-intersecting and self-touching equal
100 test_geometry<ls, ls>("LINESTRING(0 5,5 5,10 5,10 10,5 10,5 5,5 0)",
101 "LINESTRING(0 5,5 5,5 10,10 10,10 5,5 5,5 0)", "1FFF0FFF2");
102 // self-intersecting loop and self-touching equal
103 test_geometry<ls, ls>("LINESTRING(0 5,5 5,10 5,10 10,5 10,5 5,10 5,10 10,5 10,5 5,5 0)",
104 "LINESTRING(0 5,5 5,5 10,10 10,10 5,5 5,5 0)", "1FFF0FFF2");
105
106 test_geometry<ls, ls>("LINESTRING(0 0,1 1)", "LINESTRING(0 1,1 0)", "0F1FF0102");
107 test_geometry<ls, ls>("LINESTRING(0 0,1 1)", "LINESTRING(1 1,2 0)", "FF1F00102");
108 test_geometry<ls, ls>("LINESTRING(0 0,1 1)", "LINESTRING(2 0,1 1)", "FF1F00102");
109
110 test_geometry<ls, ls>("LINESTRING(0 0,1 0,2 1,3 5,4 0)", "LINESTRING(1 0,2 1,3 5)", "101FF0FF2");
111 test_geometry<ls, ls>("LINESTRING(0 0,1 0,2 1,3 5,4 0)", "LINESTRING(3 5,2 1,1 0)", "101FF0FF2");
112 test_geometry<ls, ls>("LINESTRING(1 0,2 1,3 5)", "LINESTRING(4 0,3 5,2 1,1 0,0 0)", "1FF0FF102");
113 test_geometry<ls, ls>("LINESTRING(3 5,2 1,1 0)", "LINESTRING(4 0,3 5,2 1,1 0,0 0)", "1FF0FF102");
114
115 test_geometry<ls, ls>("LINESTRING(0 0,10 0)", "LINESTRING(-1 -1,1 0,10 0,20 -1)", "1F10F0102");
116 test_geometry<ls, ls>("LINESTRING(0 0,10 0)", "LINESTRING(20 -1,10 0,1 0,-1 -1)", "1F10F0102");
117
118 test_geometry<ls, ls>("LINESTRING(-1 1,0 0,1 0,5 0,5 5,10 5,15 0,31 0)",
119 "LINESTRING(-1 -1,0 0,1 0,2 0,3 1,4 0,30 0)",
120 "101FF0102");
121 test_geometry<ls, ls>("LINESTRING(-1 1,0 0,1 0,5 0,5 5,10 5,15 0,31 0)",
122 "LINESTRING(30 0,4 0,3 1,2 0,1 0,0 0,-1 -1)",
123 "101FF0102");
124 test_geometry<ls, ls>("LINESTRING(31 0,15 0,10 5,5 5,5 0,1 0,0 0,-1 1)",
125 "LINESTRING(-1 -1,0 0,1 0,2 0,3 1,4 0,30 0)",
126 "101FF0102");
127 test_geometry<ls, ls>("LINESTRING(31 0,15 0,10 5,5 5,5 0,1 0,0 0,-1 1)",
128 "LINESTRING(30 0,4 0,3 1,2 0,1 0,0 0,-1 -1)",
129 "101FF0102");
130
131 // self-IP
132 test_geometry<ls, ls>("LINESTRING(1 0,9 0)",
133 "LINESTRING(0 0,10 0,10 10,5 0,0 10)",
134 "1FF0FF102");
135 test_geometry<ls, ls>("LINESTRING(1 0,5 0,9 0)",
136 "LINESTRING(0 0,10 0,10 10,5 0,0 10)",
137 "1FF0FF102");
138 test_geometry<ls, ls>("LINESTRING(1 0,9 0)",
139 "LINESTRING(0 0,10 0,10 10,5 10,5 -1)",
140 "1FF0FF102");
141 test_geometry<ls, ls>("LINESTRING(1 0,9 0)",
142 "LINESTRING(0 0,10 0,5 0,5 5)",
143 "1FF0FF102");
144 test_geometry<ls, ls>("LINESTRING(1 0,7 0)", "LINESTRING(0 0,10 0,10 10,4 -1)",
145 "1FF0FF102");
146 test_geometry<ls, ls>("LINESTRING(1 0,5 0,7 0)", "LINESTRING(0 0,10 0,10 10,4 -1)",
147 "1FF0FF102");
148 test_geometry<ls, ls>("LINESTRING(1 0,7 0,8 1)", "LINESTRING(0 0,10 0,10 10,4 -1)",
149 "1F10F0102");
150 test_geometry<ls, ls>("LINESTRING(1 0,5 0,7 0,8 1)", "LINESTRING(0 0,10 0,10 10,4 -1)",
151 "1F10F0102");
152
153 // self-IP going out and in on the same point
154 test_geometry<ls, ls>("LINESTRING(2 0,5 0,5 5,6 5,5 0,8 0)", "LINESTRING(1 0,9 0)",
155 "1F10FF102");
156
157 // duplicated points
158 test_geometry<ls, ls>("LINESTRING(1 1, 2 2, 2 2)", "LINESTRING(0 0, 2 2, 4 2)", "1FF0FF102");
159 test_geometry<ls, ls>("LINESTRING(1 1, 1 1, 2 2)", "LINESTRING(0 0, 2 2, 4 2)", "1FF0FF102");
160
161 // linear ring
162 test_geometry<ls, ls>("LINESTRING(0 0,10 0)", "LINESTRING(5 0,9 0,5 5,1 0,5 0)", "1F1FF01F2");
163 test_geometry<ls, ls>("LINESTRING(0 0,5 0,10 0)", "LINESTRING(5 0,9 0,5 5,1 0,5 0)", "1F1FF01F2");
164 test_geometry<ls, ls>("LINESTRING(0 0,5 0,10 0)", "LINESTRING(5 0,10 0,5 5,1 0,5 0)", "1F10F01F2");
165
166 test_geometry<ls, ls>("LINESTRING(0 0,5 0)", "LINESTRING(5 0,10 0,5 5,0 0,5 0)", "1FF0FF1F2");
167 test_geometry<ls, ls>("LINESTRING(0 0,5 0)", "LINESTRING(5 0,10 0,5 5,5 0)", "FF10F01F2");
168
169 test_geometry<ls, ls>("LINESTRING(1 0,1 6)", "LINESTRING(0 0,5 0,5 5,0 5)", "0F10F0102");
170
171 // point-size Linestring
172 test_geometry<ls, ls>("LINESTRING(1 0,1 0)", "LINESTRING(0 0,5 0)", "0FFFFF102");
173 test_geometry<ls, ls>("LINESTRING(1 0,1 0)", "LINESTRING(1 0,5 0)", "F0FFFF102");
174 test_geometry<ls, ls>("LINESTRING(1 0,1 0)", "LINESTRING(0 0,1 0)", "F0FFFF102");
175 test_geometry<ls, ls>("LINESTRING(1 0,1 0)", "LINESTRING(1 0,1 0)", "0FFFFFFF2");
176 test_geometry<ls, ls>("LINESTRING(1 0,1 0)", "LINESTRING(0 0,0 0)", "FF0FFF0F2");
177
178 //to_svg<ls, ls>("LINESTRING(0 0,5 0)", "LINESTRING(5 0,10 0,5 5,5 0)", "test_relate_00.svg");
179
180 // INVALID LINESTRINGS
181 // 1-point LS (a Point) NOT disjoint
182 //test_geometry<ls, ls>("LINESTRING(1 0)", "LINESTRING(0 0,5 0)", "0FFFFF102");
183 //test_geometry<ls, ls>("LINESTRING(0 0,5 0)", "LINESTRING(1 0)", "0F1FF0FF2");
184 //test_geometry<ls, ls>("LINESTRING(0 0,5 0)", "LINESTRING(1 0,1 0,1 0)", "0F1FF0FF2");
185 // Point/Point
186 //test_geometry<ls, ls>("LINESTRING(0 0)", "LINESTRING(0 0)", "0FFFFFFF2");
187
188 if ( BOOST_GEOMETRY_CONDITION(
189 boost::is_floating_point<typename bg::coordinate_type<ls>::type>::value ) )
190 {
191 // https://svn.boost.org/trac/boost/ticket/10904
192 // very small segments
193 test_geometry<ls, ls>("LINESTRING(5.6956521739130430148634331999347 -0.60869565217391330413931882503675,5.5 -0.50000000000000066613381477509392)",
194 "LINESTRING(5.5 -0.50000000000000066613381477509392,5.5 -0.5)",
195 "FF1F00102");
196 test_geometry<ls, ls>("LINESTRING(-3.2333333333333333925452279800083 5.5999999999999978683717927196994,-3.2333333333333333925452279800083 5.5999999999999996447286321199499)",
197 "LINESTRING(-3.2333333333333325043668082798831 5.5999999999999996447286321199499,-3.2333333333333333925452279800083 5.5999999999999996447286321199499)",
198 "FF1F00102", "F0FFFF102"); // on some platforms the first Linestring may be detected as degenerated to Point
199 }
200
201 if ( BOOST_GEOMETRY_CONDITION((
202 boost::is_same<typename bg::coordinate_type<ls>::type, double>::value )) )
203 {
204 // detected as collinear
205 test_geometry<ls, ls>("LINESTRING(1 -10, 3.069359e+307 3.069359e+307)",
206 "LINESTRING(1 6, 1 0)",
207 "*********"); // TODO: be more specific with the result
208 }
209
210 // OTHER MASKS
211 {
212 namespace bgdr = bg::detail::relate;
213 ls ls1, ls2, ls3, ls4;
214 bg::read_wkt("LINESTRING(0 0,2 0)", ls1);
215 bg::read_wkt("LINESTRING(2 0,4 0)", ls2);
216 bg::read_wkt("LINESTRING(1 0,1 1)", ls3);
217 bg::read_wkt("LINESTRING(1 0,4 0)", ls4);
218 BOOST_CHECK(bg::relate(ls1, ls2, bg::de9im::mask("FT*******")
219 || bg::de9im::mask("F**T*****")
220 || bg::de9im::mask("F***T****")));
221 BOOST_CHECK(bg::relate(ls1, ls3, bg::de9im::mask("FT*******")
222 || bg::de9im::mask("F**T*****")
223 || bg::de9im::mask("F***T****")));
224 BOOST_CHECK(bg::relate(ls3, ls1, bg::de9im::mask("FT*******")
225 || bg::de9im::mask("F**T*****")
226 || bg::de9im::mask("F***T****")));
227 BOOST_CHECK(bg::relate(ls2, ls4, bg::de9im::mask("T*F**F***"))); // within
228
229 BOOST_CHECK(bg::relate(ls1, ls2, bg::de9im::static_mask<'F','T'>()
230 || bg::de9im::static_mask<'F','*','*','T'>()
231 || bg::de9im::static_mask<'F','*','*','*','T'>()));
232 BOOST_CHECK(bg::relate(ls1, ls2, bg::de9im::mask("FT")
233 || bg::de9im::mask("F**T")
234 || bg::de9im::mask("F***T**************")));
235
236 BOOST_CHECK_THROW(bg::relate(ls1, ls2, bg::de9im::mask("A")), bg::invalid_input_exception);
237 }
238
239 // spike - boundary and interior on the same point
240 test_geometry<ls, ls>("LINESTRING(3 7, 8 8, 2 6)", "LINESTRING(5 7, 10 7, 0 7)", "0010F0102");
241
242 // 22.01.2015
243 test_geometry<ls, ls>("LINESTRING(5 5,10 10)", "LINESTRING(6 6,3 3)", "1010F0102");
244 test_geometry<ls, ls>("LINESTRING(5 5,2 8)", "LINESTRING(4 6,7 3)", "1010F0102");
245 }
246
247 template <typename P>
test_linestring_multi_linestring()248 void test_linestring_multi_linestring()
249 {
250 typedef bg::model::linestring<P> ls;
251 typedef bg::model::multi_linestring<ls> mls;
252
253 // LS disjoint
254 test_geometry<ls, mls>("LINESTRING(0 0,10 0)", "MULTILINESTRING((1 0,2 0),(1 1,2 1))", "101FF0102");
255 // linear ring disjoint
256 test_geometry<ls, mls>("LINESTRING(0 0,10 0)", "MULTILINESTRING((1 0,2 0),(1 1,2 1,2 2,1 1))", "101FF01F2");
257 // 2xLS forming non-simple linear ring disjoint
258 test_geometry<ls, mls>("LINESTRING(0 0,10 0)", "MULTILINESTRING((1 0,2 0),(1 1,2 1,2 2),(1 1,2 2))", "101FF01F2");
259
260 test_geometry<ls, mls>("LINESTRING(0 0,10 0)",
261 "MULTILINESTRING((1 0,9 0),(9 0,2 0))",
262 "101FF0FF2");
263
264 // rings
265 test_geometry<ls, mls>("LINESTRING(0 0,5 0,5 5,0 5,0 0)",
266 "MULTILINESTRING((5 5,0 5,0 0),(0 0,5 0,5 5))",
267 "1FFFFFFF2");
268 test_geometry<ls, mls>("LINESTRING(0 0,5 0,5 5,0 5,0 0)",
269 "MULTILINESTRING((5 5,5 0,0 0),(0 0,0 5,5 5))",
270 "1FFFFFFF2");
271 // overlapping rings
272 test_geometry<ls, mls>("LINESTRING(0 0,5 0,5 5,0 5,0 0)",
273 "MULTILINESTRING((5 5,0 5,0 0),(0 0,5 0,5 5,0 5))",
274 "10FFFFFF2");
275 test_geometry<ls, mls>("LINESTRING(0 0,5 0,5 5,0 5,0 0)",
276 "MULTILINESTRING((5 5,5 0,0 0),(0 0,0 5,5 5,5 0))",
277 "10FFFFFF2");
278
279 // INVALID LINESTRINGS
280 // 1-point LS (a Point) disjoint
281 //test_geometry<ls, mls>("LINESTRING(0 0,10 0)", "MULTILINESTRING((1 0,2 0),(1 1))", "101FF00F2");
282 //test_geometry<ls, mls>("LINESTRING(0 0,10 0)", "MULTILINESTRING((1 0,2 0),(1 1,1 1))", "101FF00F2");
283 //test_geometry<ls, mls>("LINESTRING(0 0,10 0)", "MULTILINESTRING((1 0,2 0),(1 1,1 1,1 1))", "101FF00F2");
284 // 1-point LS (a Point) NOT disjoint
285 //test_geometry<ls, mls>("LINESTRING(0 0,10 0)", "MULTILINESTRING((1 0,9 0),(2 0))", "101FF0FF2");
286 //test_geometry<ls, mls>("LINESTRING(0 0,10 0)", "MULTILINESTRING((1 0,9 0),(2 0,2 0))", "101FF0FF2");
287 //test_geometry<ls, mls>("LINESTRING(0 0,10 0)", "MULTILINESTRING((1 0,9 0),(2 0,2 0,2 0))", "101FF0FF2");
288
289 // point-like
290 test_geometry<ls, mls>("LINESTRING(0 0, 5 0)", // |--------------|
291 "MULTILINESTRING((0 0, 1 0),(2 0, 2 0))", // |------| *
292 "101F00FF2");
293 test_geometry<ls, mls>("LINESTRING(0 0, 5 0)", // |--------------|
294 "MULTILINESTRING((0 0, 1 0),(1 0, 1 0))", // |------*
295 "101F00FF2");
296 test_geometry<ls, mls>("LINESTRING(0 0, 5 0)", // |--------------|
297 "MULTILINESTRING((5 0, 1 0),(1 0, 1 0))", // *-------|
298 "101F00FF2");
299 test_geometry<ls, mls>("LINESTRING(0 0, 5 0)", // |--------------|
300 "MULTILINESTRING((0 0, 1 0),(5 0, 5 0))", // |------| *
301 "10100FFF2");
302 test_geometry<ls, mls>("LINESTRING(0 0, 5 0)", // |--------------|
303 "MULTILINESTRING((0 0, 1 0),(0 0, 0 0))", // *------|
304 "101000FF2");
305 test_geometry<ls, mls>("LINESTRING(0 0, 5 0)", // |--------------|
306 "MULTILINESTRING((4 0, 5 0),(5 0, 5 0))", // |------*
307 "101000FF2");
308 test_geometry<ls, mls>("LINESTRING(0 0, 5 0)", // |--------------|
309 "MULTILINESTRING((1 0, 2 0),(0 0, 0 0))", // * |------|
310 "1010F0FF2");
311
312 test_geometry<ls, mls>("LINESTRING(0 0, 5 0)", // |--------------|
313 "MULTILINESTRING((2 0, 2 0),(2 0, 2 2))", // *
314 "001FF0102"); // |
315
316 // for consistency
317 test_geometry<ls, mls>("LINESTRING(0 0, 5 0)", // |--------------|
318 "MULTILINESTRING((0 0, 5 0),(0 0, 2 0))", // |--------------|
319 "10F00FFF2"); // |------|
320
321 test_geometry<ls, mls>("LINESTRING(0 0, 5 0)", // |--------------|
322 "MULTILINESTRING((0 0, 5 0),(3 0, 5 0))", // |--------------|
323 "10F00FFF2"); // |------|
324
325 test_geometry<ls, mls>("LINESTRING(0 0, 5 0)", // |--------------|
326 "MULTILINESTRING((0 0, 5 0),(0 0, 6 0))", // |--------------|
327 "1FF00F102"); // |----------------|
328
329 test_geometry<ls, mls>("LINESTRING(0 0, 5 0)", // |--------------|
330 "MULTILINESTRING((0 0, 5 0),(-1 0, 5 0))", // |--------------|
331 "1FF00F102"); // |----------------|
332
333 test_geometry<ls, mls>("LINESTRING(0 0, 5 0)", // |--------------|
334 "MULTILINESTRING((0 0, 5 0),(-1 0, 6 0))", // |--------------|
335 "1FF00F102"); // |------------------|
336
337 test_geometry<ls, mls>("LINESTRING(0 0, 5 0)", // |--------------|
338 "MULTILINESTRING((0 0, 5 0),(-1 0, 2 0))", // |--------------|
339 "10F00F102"); // |-------|
340
341 test_geometry<ls, mls>("LINESTRING(0 0, 5 0)", // |--------------|
342 "MULTILINESTRING((0 0, 5 0),(2 0, 6 0))", // |--------------|
343 "10F00F102"); // |-------|
344
345 test_geometry<ls, mls>("LINESTRING(0 0, 5 0)", // |--------------|
346 "MULTILINESTRING((0 0, 5 0),(2 0, 2 2))", // |--------------|
347 "10FF0F102"); // |
348 // |
349
350 if ( BOOST_GEOMETRY_CONDITION(boost::is_floating_point<typename bg::coordinate_type<ls>::type>::value) )
351 {
352 // related to https://svn.boost.org/trac/boost/ticket/10904
353 test_geometry<ls, mls>("LINESTRING(-2305843009213693956 4611686018427387906, -33 -92, 78 83)",
354 "MULTILINESTRING((20 100, 31 -97, -46 57, -20 -4),(-71 -4))",
355 "*********"); // TODO: be more specific with the result
356 }
357
358 // 22.01.2015
359 // inspired by L/A and A/A
360 test_geometry<ls, mls>("LINESTRING(1 1,2 2)", "MULTILINESTRING((0 0,1 1),(1 1,3 3))", "1FF0FF102");
361
362 // 25.01.2015
363 test_geometry<ls, mls>("LINESTRING(1 1, 5 5, 4 4)",
364 "MULTILINESTRING((2 5, 7 5, 8 3, 6 3, 4 0),(0 0,10 10))",
365 "1FF0FF102");
366 test_geometry<ls, mls>("LINESTRING(4 1, 4 5, 4 4)",
367 "MULTILINESTRING((2 5, 7 5, 8 3, 6 3, 4 0),(4 0, 4 8, 0 4))",
368 "1FF0FF102");
369 test_geometry<ls, mls>("LINESTRING(1 1,5 5,4 4)",
370 "MULTILINESTRING((5 0,5 5,5 10),(0 0,10 10))",
371 "1FF0FF102");
372 test_geometry<ls, mls>("LINESTRING(1 1,5 5,1 0)",
373 "MULTILINESTRING((5 0,5 5,5 10),(0 0,5 5,1 0))",
374 "1FF00F102");
375 test_geometry<ls, mls>("LINESTRING(5 5,4 4)",
376 "MULTILINESTRING((5 0,5 5,5 10),(0 0,10 10))",
377 "1FF0FF102");
378 test_geometry<ls, mls>("LINESTRING(6 6,5 5,4 4)",
379 "MULTILINESTRING((5 0,5 5,5 10),(0 0,10 10))",
380 "1FF0FF102");
381 test_geometry<ls, mls>("LINESTRING(5 5,4 4)",
382 "MULTILINESTRING((5 0,5 5,5 10))",
383 "FF10F0102");
384 }
385
386 template <typename P>
test_multi_linestring_multi_linestring()387 void test_multi_linestring_multi_linestring()
388 {
389 typedef bg::model::linestring<P> ls;
390 typedef bg::model::multi_linestring<ls> mls;
391
392 test_geometry<mls, mls>("MULTILINESTRING((0 0,0 0,18 0,18 0,19 0,19 0,19 0,30 0,30 0))",
393 "MULTILINESTRING((0 10,5 0,20 0,20 0,30 0))",
394 "1F1F00102");
395 test_geometry<mls, mls>("MULTILINESTRING((0 0,0 0,18 0,18 0,19 0,19 0,19 0,30 0,30 0))",
396 //"MULTILINESTRING((0 10,5 0,20 0,20 0,30 0),(1 10,1 10,1 0,1 0,1 -10),(2 0,2 0),(3 0,3 0,3 0),(0 0,0 0,0 10,0 10),(30 0,30 0,31 0,31 0))",
397 "MULTILINESTRING((0 10,5 0,20 0,20 0,30 0),(1 10,1 10,1 0,1 0,1 -10),(0 0,0 0,0 10,0 10),(30 0,30 0,31 0,31 0))",
398 "1F100F102");
399 test_geometry<mls, mls>("MULTILINESTRING((0 0,0 0,18 0,18 0,19 0,19 0,19 0,30 0,30 0))",
400 "MULTILINESTRING((0 10,5 0,20 0,20 0,30 0),(0 0,0 0,0 10,0 10))",
401 "1F1F0F1F2");
402
403 // point-like
404 test_geometry<mls, mls>("MULTILINESTRING((0 0, 0 0),(1 1, 1 1))",
405 "MULTILINESTRING((0 0, 0 0))",
406 "0F0FFFFF2");
407 test_geometry<mls, mls>("MULTILINESTRING((0 0, 0 0),(1 1, 1 1))",
408 "MULTILINESTRING((0 0, 0 0),(1 1, 1 1))",
409 "0FFFFFFF2");
410 test_geometry<mls, mls>("MULTILINESTRING((0 0, 0 0),(1 1, 1 1))",
411 "MULTILINESTRING((2 2, 2 2),(3 3, 3 3))",
412 "FF0FFF0F2");
413
414 test_geometry<mls, mls>("MULTILINESTRING((0 5,10 5,10 10,5 10),(5 10,5 0,5 2),(5 2,5 5,0 5))",
415 "MULTILINESTRING((5 5,0 5),(5 5,5 0),(10 10,10 5,5 5,5 10,10 10))",
416 "10FFFFFF2");
417
418 if ( BOOST_GEOMETRY_CONDITION((
419 boost::is_same<typename bg::coordinate_type<P>::type, double>::value )) )
420 {
421 // assertion failure in 1.57
422 test_geometry<mls, mls>("MULTILINESTRING((-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))",
423 "MULTILINESTRING((-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))",
424 "*********");
425 }
426 }
427
428 template <typename P>
test_all()429 void test_all()
430 {
431 test_linestring_linestring<P>();
432 test_linestring_multi_linestring<P>();
433 test_multi_linestring_multi_linestring<P>();
434 }
435
test_main(int,char * [])436 int test_main( int , char* [] )
437 {
438 test_all<bg::model::d2::point_xy<int> >();
439 test_all<bg::model::d2::point_xy<double> >();
440
441 #if defined(HAVE_TTMATH)
442 test_all<bg::model::d2::point_xy<ttmath_big> >();
443 #endif
444
445 return 0;
446 }
447