[ceph.git] / ceph / src / boost / libs / geometry / test / algorithms / relational_operations / disjoint / disjoint_sph.cpp

// Boost.Geometry

// Copyright (c) 2016 Oracle and/or its affiliates.
// Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle

// Use, modification and distribution is subject to 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)

#include "test_disjoint.hpp"

#include <algorithms/overlay/overlay_cases.hpp>
#include <algorithms/overlay/multi_overlay_cases.hpp>

#include <boost/geometry/geometries/geometries.hpp>

#include <boost/geometry/strategies/strategies.hpp>


template <typename P>
void test_polygon_polygon()
{
    typedef bg::model::polygon<P> poly;
    typedef bg::model::ring<P> ring;

    test_geometry<ring, ring>(case_1[0], case_1[1],
                              false);
    test_geometry<ring, poly>(case_1[0], case_1[1],
                              false);

    test_geometry<poly, poly>(case_1[0], case_1[1],
                              false);
    test_geometry<poly, poly>(case_2[0], case_2[1],
                              false);
    test_geometry<poly, poly>(case_3_sph[0], case_3_sph[1],
                              false);
    test_geometry<poly, poly>(case_3_2_sph[0], case_3_2_sph[1],
                              false);
    test_geometry<poly, poly>(case_4[0], case_4[1],
                              false);
    test_geometry<poly, poly>(case_5[0], case_5[1],
                              false);
    test_geometry<poly, poly>(case_6_sph[0], case_6_sph[1],
                              false);

    test_geometry<poly, poly>(case_7[0], case_7[1],
                              false);
    test_geometry<poly, poly>(case_8_sph[0], case_8_sph[1],
                              false);
    test_geometry<poly, poly>(case_9_sph[0], case_9_sph[1],
                              false);
    test_geometry<poly, poly>(case_10_sph[0], case_10_sph[1],
                              false);
    test_geometry<poly, poly>(case_11_sph[0], case_11_sph[1],
                              false);
    test_geometry<poly, poly>(case_12[0], case_12[1],
                              false);

    test_geometry<poly, poly>(case_13_sph[0], case_13_sph[1],
                              false);
    test_geometry<poly, poly>(case_14_sph[0], case_14_sph[1],
                              false);
    test_geometry<poly, poly>(case_15_sph[0], case_15_sph[1],
                              false);
    test_geometry<poly, poly>(case_16_sph[0], case_16_sph[1],
                              false);
    test_geometry<poly, poly>(case_17_sph[0], case_17_sph[1],
                              false);
    test_geometry<poly, poly>(case_17_sph[1], case_17_sph[0],
                              false);
    test_geometry<poly, poly>(case_18_sph[0], case_18_sph[1],
                              false);
    test_geometry<poly, poly>(case_18_sph[1], case_18_sph[0],
                              false);
}

template <typename P>
void test_polygon_multi_polygon()
{
    typedef bg::model::polygon<P> poly;
    typedef bg::model::ring<P> ring;
    typedef bg::model::multi_polygon<poly> mpoly;

    test_geometry<ring, mpoly>(case_1[0], case_multi_2[0],
                               false);
    test_geometry<poly, mpoly>(case_2[0], case_multi_2[0],
                               false);
}

template <typename P>
void test_multi_polygon_multi_polygon()
{
    typedef bg::model::polygon<P> poly;
    typedef bg::model::multi_polygon<poly> mpoly;

    test_geometry<mpoly, mpoly>(case_multi_2[0], case_multi_2[1],
                                false);
}

template <typename P>
void test_linestring_polygon()
{
    typedef bg::model::linestring<P> ls;
    typedef bg::model::polygon<P> poly;
    typedef bg::model::polygon<P> ring;

    test_geometry<ls, poly>("LINESTRING(11 0,11 10)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", true);
    test_geometry<ls, ring>("LINESTRING(11 0,11 10)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", true);
    test_geometry<ls, poly>("LINESTRING(0 0,10 10)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", false);
    test_geometry<ls, poly>("LINESTRING(5 0,5 5,10 5)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", false);
    test_geometry<ls, poly>("LINESTRING(5 1,5 5,9 5)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", false);
    test_geometry<ls, poly>("LINESTRING(11 1,11 5)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", true);

    test_geometry<ls, poly>("LINESTRING(9 1,10 5,9 9)",
                            "POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5))",
                            false);

    test_geometry<ls, poly>("LINESTRING(9 1,10 5,9 9,1 9,1 1,9 1)",
                            "POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5))",
                            false);

    test_geometry<ls, poly>("LINESTRING(0 0,10 0,10 10,0 10,0 0)",
                            "POLYGON((0 0,0 10,10 10,10 0,0 0))",
                            false);
}

template <typename P>
void test_linestring_multi_polygon()
{
    typedef bg::model::linestring<P> ls;
    typedef bg::model::polygon<P> poly;
    typedef bg::model::multi_polygon<poly> mpoly;

    test_geometry<ls, mpoly>("LINESTRING(10 1,10 5,10 9)",
                             "MULTIPOLYGON(((0 20,0 30,10 30,10 20,0 20)),((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5)))",
                             false);
}

template <typename P>
void test_multi_linestring_polygon()
{
    typedef bg::model::linestring<P> ls;
    typedef bg::model::polygon<P> poly;
    typedef bg::model::ring<P> ring;
    typedef bg::model::multi_linestring<ls> mls;

    test_geometry<mls, poly>("MULTILINESTRING((11 11, 20 20),(5 7, 4 1))",
                             "POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,4 2,4 4,2 4,2 2))",
                             false);

    test_geometry<mls, ring>("MULTILINESTRING((6 6,15 15),(0 0, 7 7))",
                             "POLYGON((5 5,5 15,15 15,15 5,5 5))",
                             false);

    test_geometry<mls, poly>("MULTILINESTRING((3 10.031432746397092, 1 5, 1 10.013467818052765, 3 4, 7 8, 6 10.035925377760330, 10 2))",
                             "POLYGON((0 0,0 10,10 10,10 0,0 0))",
                             false);
}

template <typename P>
void test_multi_linestring_multi_polygon()
{
    typedef bg::model::linestring<P> ls;
    typedef bg::model::polygon<P> poly;
    typedef bg::model::multi_linestring<ls> mls;
    typedef bg::model::multi_polygon<poly> mpoly;

    test_geometry<mls, mpoly>("MULTILINESTRING((0 0,10 0,10 10,0 10,0 0),(2 2,5 5,2 8,2 2))",
                              "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0),(2 2,5 5,2 8,2 2)))",
                              false);

    test_geometry<mls, mpoly>("MULTILINESTRING((0 0,10 0,10 10),(10 10,0 10,0 0),(20 20,50 50,20 80,20 20))",
                              "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)))",
                              false);

    test_geometry<mls, mpoly>("MULTILINESTRING((5 -2,4 -2,5 0),(5 -2,6 -2,5 0))",
                              "MULTIPOLYGON(((5 0,0 5,10 5,5 0)),((5 0,10 -5,0 -5,5 0)))",
                              false);
}

template <typename P>
void test_linestring_linestring()
{
    typedef bg::model::linestring<P> ls;

    test_geometry<ls, ls>("LINESTRING(0 0, 2 2, 3 2)", "LINESTRING(0 0, 2 2, 3 2)", false);

    test_geometry<ls, ls>("LINESTRING(1 0,2 2,2 3)", "LINESTRING(0 0, 2 2, 3 2)", false);
}

template <typename P>
void test_linestring_multi_linestring()
{
    typedef bg::model::linestring<P> ls;
    typedef bg::model::multi_linestring<ls> mls;

    test_geometry<ls, mls>("LINESTRING(0 0,10 0)",
                           "MULTILINESTRING((1 0,2 0),(1 1,2 1))",
                           false);

    test_geometry<ls, mls>("LINESTRING(0 0,5 0,5 5,0 5,0 0)",
                           "MULTILINESTRING((5 5,0 5,0 0),(0 0,5 0,5 5))",
                           false);
}

template <typename P>
void test_multi_linestring_multi_linestring()
{
    typedef bg::model::linestring<P> ls;
    typedef bg::model::multi_linestring<ls> mls;

    test_geometry<mls, mls>("MULTILINESTRING((0 0,0 0,18 0,18 0,19 0,19 0,19 0,30 0,30 0))",
                            "MULTILINESTRING((0 10,5 0,20 0,20 0,30 0))",
                            false);
}

template <typename P>
void test_point_polygon()
{
    typedef bg::model::polygon<P> poly;
    
    // https://svn.boost.org/trac/boost/ticket/9162
    test_geometry<P, poly>("POINT(0 90)",
                           "POLYGON((0 80,-90 80, -180 80, 90 80, 0 80))",
                           false);
    test_geometry<P, poly>("POINT(-120 21)",
                           "POLYGON((30 0,30 30,90 30, 90 0, 30 0))",
                           true);
    // extended
    test_geometry<P, poly>("POINT(0 -90)",
                           "POLYGON((0 -80,90 -80, -180 -80, -90 -80, 0 -80))",
                           false);
    test_geometry<P, poly>("POINT(0 89)",
                           "POLYGON((0 80,-90 80, -180 80, 90 80, 0 80))",
                           false);
    test_geometry<P, poly>("POINT(-180 89)",
                           "POLYGON((0 80,-90 80, -180 80, 90 80, 0 80))",
                           false);
}


template <typename P>
void test_all()
{
    test_polygon_polygon<P>();
    test_polygon_multi_polygon<P>();
    test_multi_polygon_multi_polygon<P>();

    test_linestring_polygon<P>();
    test_linestring_multi_polygon<P>();
    test_multi_linestring_polygon<P>();
    test_multi_linestring_multi_polygon<P>();

    test_linestring_linestring<P>();
    test_linestring_multi_linestring<P>();
    test_multi_linestring_multi_linestring<P>();

    test_point_polygon<P>();
}


int test_main( int , char* [] )
{
    test_all<bg::model::point<double, 2, bg::cs::spherical_equatorial<bg::degree> > >();

#if defined(HAVE_TTMATH)
    test_cs<bg::model::point<ttmath_big, 2, bg::cs::spherical_equatorial<bg::degree> > >();
#endif

    return 0;
}
Commit	Line	Data
7c673cae FG	1	// Boost.Geometry
	2
	3	// Copyright (c) 2016 Oracle and/or its affiliates.
	4	// Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
	5
	6	// Use, modification and distribution is subject to the Boost Software License,
	7	// Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
	8	// http://www.boost.org/LICENSE_1_0.txt)
	9
	10	#include "test_disjoint.hpp"
	11
	12	#include <algorithms/overlay/overlay_cases.hpp>
	13	#include <algorithms/overlay/multi_overlay_cases.hpp>
	14
	15	#include <boost/geometry/geometries/geometries.hpp>
	16
	17	#include <boost/geometry/strategies/strategies.hpp>
	18
	19
	20	template <typename P>
	21	void test_polygon_polygon()
	22	{
	23	typedef bg::model::polygon<P> poly;
	24	typedef bg::model::ring<P> ring;
	25
	26	test_geometry<ring, ring>(case_1[0], case_1[1],
	27	false);
	28	test_geometry<ring, poly>(case_1[0], case_1[1],
	29	false);
	30
	31	test_geometry<poly, poly>(case_1[0], case_1[1],
	32	false);
	33	test_geometry<poly, poly>(case_2[0], case_2[1],
	34	false);
	35	test_geometry<poly, poly>(case_3_sph[0], case_3_sph[1],
	36	false);
	37	test_geometry<poly, poly>(case_3_2_sph[0], case_3_2_sph[1],
	38	false);
	39	test_geometry<poly, poly>(case_4[0], case_4[1],
	40	false);
	41	test_geometry<poly, poly>(case_5[0], case_5[1],
	42	false);
	43	test_geometry<poly, poly>(case_6_sph[0], case_6_sph[1],
	44	false);
	45
	46	test_geometry<poly, poly>(case_7[0], case_7[1],
	47	false);
	48	test_geometry<poly, poly>(case_8_sph[0], case_8_sph[1],
	49	false);
	50	test_geometry<poly, poly>(case_9_sph[0], case_9_sph[1],
	51	false);
	52	test_geometry<poly, poly>(case_10_sph[0], case_10_sph[1],
	53	false);
	54	test_geometry<poly, poly>(case_11_sph[0], case_11_sph[1],
	55	false);
	56	test_geometry<poly, poly>(case_12[0], case_12[1],
	57	false);
	58
	59	test_geometry<poly, poly>(case_13_sph[0], case_13_sph[1],
	60	false);
	61	test_geometry<poly, poly>(case_14_sph[0], case_14_sph[1],
	62	false);
	63	test_geometry<poly, poly>(case_15_sph[0], case_15_sph[1],
	64	false);
65	test_geometry<poly, poly>(case_16_sph[0], case_16_sph[1],
66	false);
67	test_geometry<poly, poly>(case_17_sph[0], case_17_sph[1],
68	false);
69	test_geometry<poly, poly>(case_17_sph[1], case_17_sph[0],
70	false);
71	test_geometry<poly, poly>(case_18_sph[0], case_18_sph[1],
72	false);
73	test_geometry<poly, poly>(case_18_sph[1], case_18_sph[0],
74	false);
75	}
76
77	template <typename P>
78	void test_polygon_multi_polygon()
79	{
80	typedef bg::model::polygon<P> poly;
81	typedef bg::model::ring<P> ring;
82	typedef bg::model::multi_polygon<poly> mpoly;
83
84	test_geometry<ring, mpoly>(case_1[0], case_multi_2[0],
85	false);
86	test_geometry<poly, mpoly>(case_2[0], case_multi_2[0],
87	false);
88	}
89
90	template <typename P>
91	void test_multi_polygon_multi_polygon()
92	{
93	typedef bg::model::polygon<P> poly;
94	typedef bg::model::multi_polygon<poly> mpoly;
95
96	test_geometry<mpoly, mpoly>(case_multi_2[0], case_multi_2[1],
97	false);
98	}
99
100	template <typename P>
101	void test_linestring_polygon()
102	{
103	typedef bg::model::linestring<P> ls;
104	typedef bg::model::polygon<P> poly;
105	typedef bg::model::polygon<P> ring;
106
107	test_geometry<ls, poly>("LINESTRING(11 0,11 10)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", true);
108	test_geometry<ls, ring>("LINESTRING(11 0,11 10)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", true);
109	test_geometry<ls, poly>("LINESTRING(0 0,10 10)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", false);
110	test_geometry<ls, poly>("LINESTRING(5 0,5 5,10 5)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", false);
111	test_geometry<ls, poly>("LINESTRING(5 1,5 5,9 5)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", false);
112	test_geometry<ls, poly>("LINESTRING(11 1,11 5)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", true);
113
114	test_geometry<ls, poly>("LINESTRING(9 1,10 5,9 9)",
115	"POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5))",
116	false);
117
118	test_geometry<ls, poly>("LINESTRING(9 1,10 5,9 9,1 9,1 1,9 1)",
119	"POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5))",
120	false);
121
122	test_geometry<ls, poly>("LINESTRING(0 0,10 0,10 10,0 10,0 0)",
123	"POLYGON((0 0,0 10,10 10,10 0,0 0))",
124	false);
125	}
126
127	template <typename P>
128	void test_linestring_multi_polygon()
129	{
130	typedef bg::model::linestring<P> ls;
131	typedef bg::model::polygon<P> poly;
132	typedef bg::model::multi_polygon<poly> mpoly;
133
134	test_geometry<ls, mpoly>("LINESTRING(10 1,10 5,10 9)",
135	"MULTIPOLYGON(((0 20,0 30,10 30,10 20,0 20)),((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5)))",
136	false);
137	}
138
139	template <typename P>
140	void test_multi_linestring_polygon()
141	{
142	typedef bg::model::linestring<P> ls;
143	typedef bg::model::polygon<P> poly;
144	typedef bg::model::ring<P> ring;
145	typedef bg::model::multi_linestring<ls> mls;
146
147	test_geometry<mls, poly>("MULTILINESTRING((11 11, 20 20),(5 7, 4 1))",
148	"POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,4 2,4 4,2 4,2 2))",
149	false);
150
151	test_geometry<mls, ring>("MULTILINESTRING((6 6,15 15),(0 0, 7 7))",
152	"POLYGON((5 5,5 15,15 15,15 5,5 5))",
153	false);
154
155	test_geometry<mls, poly>("MULTILINESTRING((3 10.031432746397092, 1 5, 1 10.013467818052765, 3 4, 7 8, 6 10.035925377760330, 10 2))",
156	"POLYGON((0 0,0 10,10 10,10 0,0 0))",
157	false);
158	}
159
160	template <typename P>
161	void test_multi_linestring_multi_polygon()
162	{
163	typedef bg::model::linestring<P> ls;
164	typedef bg::model::polygon<P> poly;
165	typedef bg::model::multi_linestring<ls> mls;
166	typedef bg::model::multi_polygon<poly> mpoly;
167
168	test_geometry<mls, mpoly>("MULTILINESTRING((0 0,10 0,10 10,0 10,0 0),(2 2,5 5,2 8,2 2))",
169	"MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0),(2 2,5 5,2 8,2 2)))",
170	false);
171
172	test_geometry<mls, mpoly>("MULTILINESTRING((0 0,10 0,10 10),(10 10,0 10,0 0),(20 20,50 50,20 80,20 20))",
173	"MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)))",
174	false);
175
176	test_geometry<mls, mpoly>("MULTILINESTRING((5 -2,4 -2,5 0),(5 -2,6 -2,5 0))",
177	"MULTIPOLYGON(((5 0,0 5,10 5,5 0)),((5 0,10 -5,0 -5,5 0)))",
178	false);
179	}
180
181	template <typename P>
182	void test_linestring_linestring()
183	{
184	typedef bg::model::linestring<P> ls;
185
186	test_geometry<ls, ls>("LINESTRING(0 0, 2 2, 3 2)", "LINESTRING(0 0, 2 2, 3 2)", false);
187
188	test_geometry<ls, ls>("LINESTRING(1 0,2 2,2 3)", "LINESTRING(0 0, 2 2, 3 2)", false);
189	}
190
191	template <typename P>
192	void test_linestring_multi_linestring()
193	{
194	typedef bg::model::linestring<P> ls;
195	typedef bg::model::multi_linestring<ls> mls;
196
197	test_geometry<ls, mls>("LINESTRING(0 0,10 0)",
198	"MULTILINESTRING((1 0,2 0),(1 1,2 1))",
199	false);
200
201	test_geometry<ls, mls>("LINESTRING(0 0,5 0,5 5,0 5,0 0)",
202	"MULTILINESTRING((5 5,0 5,0 0),(0 0,5 0,5 5))",
203	false);
204	}
205
206	template <typename P>
207	void test_multi_linestring_multi_linestring()
208	{
209	typedef bg::model::linestring<P> ls;
210	typedef bg::model::multi_linestring<ls> mls;
211
212	test_geometry<mls, mls>("MULTILINESTRING((0 0,0 0,18 0,18 0,19 0,19 0,19 0,30 0,30 0))",
213	"MULTILINESTRING((0 10,5 0,20 0,20 0,30 0))",
214	false);
215	}
216
217	template <typename P>
218	void test_point_polygon()
219	{
220	typedef bg::model::polygon<P> poly;
221
222	// https://svn.boost.org/trac/boost/ticket/9162
223	test_geometry<P, poly>("POINT(0 90)",
224	"POLYGON((0 80,-90 80, -180 80, 90 80, 0 80))",
225	false);
226	test_geometry<P, poly>("POINT(-120 21)",
227	"POLYGON((30 0,30 30,90 30, 90 0, 30 0))",
228	true);
229	// extended
230	test_geometry<P, poly>("POINT(0 -90)",
231	"POLYGON((0 -80,90 -80, -180 -80, -90 -80, 0 -80))",
232	false);
233	test_geometry<P, poly>("POINT(0 89)",
234	"POLYGON((0 80,-90 80, -180 80, 90 80, 0 80))",
235	false);
236	test_geometry<P, poly>("POINT(-180 89)",
237	"POLYGON((0 80,-90 80, -180 80, 90 80, 0 80))",
238	false);
239	}
240
241
242	template <typename P>
243	void test_all()
244	{
245	test_polygon_polygon<P>();
246	test_polygon_multi_polygon<P>();
247	test_multi_polygon_multi_polygon<P>();
248
249	test_linestring_polygon<P>();
250	test_linestring_multi_polygon<P>();
251	test_multi_linestring_polygon<P>();
252	test_multi_linestring_multi_polygon<P>();
253
254	test_linestring_linestring<P>();
255	test_linestring_multi_linestring<P>();
256	test_multi_linestring_multi_linestring<P>();
257
258	test_point_polygon<P>();
259	}
260
261
262	int test_main( int , char* [] )
263	{
264	test_all<bg::model::point<double, 2, bg::cs::spherical_equatorial<bg::degree> > >();
265
266	#if defined(HAVE_TTMATH)
267	test_cs<bg::model::point<ttmath_big, 2, bg::cs::spherical_equatorial<bg::degree> > >();
268	#endif
269
270	return 0;
271	}