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

   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 //https://svn.boost.org/trac10/ticket/13057
 192 template <typename P>
 193 void test_linestring_linestring_radians()
 194 {
 195     typedef bg::model::linestring<P> ls;
 196
 197     test_geometry<ls, ls>("LINESTRING(0 -0.31415926535897897853,\
 198                                       0.26179938779914918578 0,\
 199                                      -0.034906585039886556254 0.13962634015954622502,\
 200                                      -0.12217304763960294689 0.12217304763960294689)",\
 201                           "LINESTRING(-0.034906585039886556254 0.13962634015954622502,\
 202                                       -0.26179938779914918578 0)", false);
 203
 204 }
 205
 206 template <typename P>
 207 void test_linestring_multi_linestring()
 208 {
 209     typedef bg::model::linestring<P> ls;
 210     typedef bg::model::multi_linestring<ls> mls;
 211
 212     test_geometry<ls, mls>("LINESTRING(0 0,10 0)",
 213                            "MULTILINESTRING((1 0,2 0),(1 1,2 1))",
 214                            false);
 215
 216     test_geometry<ls, mls>("LINESTRING(0 0,5 0,5 5,0 5,0 0)",
 217                            "MULTILINESTRING((5 5,0 5,0 0),(0 0,5 0,5 5))",
 218                            false);
 219 }
 220
 221 template <typename P>
 222 void test_multi_linestring_multi_linestring()
 223 {
 224     typedef bg::model::linestring<P> ls;
 225     typedef bg::model::multi_linestring<ls> mls;
 226
 227     test_geometry<mls, mls>("MULTILINESTRING((0 0,0 0,18 0,18 0,19 0,19 0,19 0,30 0,30 0))",
 228                             "MULTILINESTRING((0 10,5 0,20 0,20 0,30 0))",
 229                             false);
 230 }
 231
 232 template <typename P>
 233 void test_point_polygon()
 234 {
 235     typedef bg::model::polygon<P> poly;
 236
 237     // https://svn.boost.org/trac/boost/ticket/9162
 238     test_geometry<P, poly>("POINT(0 90)",
 239                            "POLYGON((0 80,-90 80, -180 80, 90 80, 0 80))",
 240                            false);
 241     test_geometry<P, poly>("POINT(-120 21)",
 242                            "POLYGON((30 0,30 30,90 30, 90 0, 30 0))",
 243                            true);
 244     // extended
 245     test_geometry<P, poly>("POINT(0 -90)",
 246                            "POLYGON((0 -80,90 -80, -180 -80, -90 -80, 0 -80))",
 247                            false);
 248     test_geometry<P, poly>("POINT(0 89)",
 249                            "POLYGON((0 80,-90 80, -180 80, 90 80, 0 80))",
 250                            false);
 251     test_geometry<P, poly>("POINT(-180 89)",
 252                            "POLYGON((0 80,-90 80, -180 80, 90 80, 0 80))",
 253                            false);
 254 }
 255
 256
 257 template <typename P>
 258 void test_all()
 259 {
 260     test_polygon_polygon<P>();
 261     test_polygon_multi_polygon<P>();
 262     test_multi_polygon_multi_polygon<P>();
 263
 264     test_linestring_polygon<P>();
 265     test_linestring_multi_polygon<P>();
 266     test_multi_linestring_polygon<P>();
 267     test_multi_linestring_multi_polygon<P>();
 268
 269     test_linestring_linestring<P>();
 270     test_linestring_multi_linestring<P>();
 271     test_multi_linestring_multi_linestring<P>();
 272
 273     test_linestring_linestring_radians<bg::model::point<double, 2,
 274                                      bg::cs::spherical_equatorial<bg::radian> > >();
 275
 276     test_point_polygon<P>();
 277 }
 278
 279
 280 int test_main( int , char* [] )
 281 {
 282     test_all<bg::model::point<double, 2, bg::cs::spherical_equatorial<bg::degree> > >();
 283
 284 #if defined(HAVE_TTMATH)
 285     test_cs<bg::model::point<ttmath_big, 2, bg::cs::spherical_equatorial<bg::degree> > >();
 286 #endif
 287
 288     return 0;
 289 }