ceph/src/boost/libs/geometry/test/strategies/winding.cpp

   1 // Boost.Geometry (aka GGL, Generic Geometry Library)
   2 // Unit Test
   3
   4 // Copyright (c) 2010-2012 Barend Gehrels, Amsterdam, the Netherlands.
   5
   6 // This file was modified by Oracle on 2014, 2016.
   7 // Modifications copyright (c) 2014-2016 Oracle and/or its affiliates.
   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
  15 #include <strategies/test_within.hpp>
  16
  17
  18 template <typename Point>
  19 void test_cartesian()
  20 {
  21     typedef bg::model::polygon<Point> polygon;
  22
  23     std::string const box = "POLYGON((0 0,0 2,2 2,2 0,0 0))";
  24     std::string const triangle = "POLYGON((0 0,0 4,6 0,0 0))";
  25     std::string const with_hole = "POLYGON((0 0,0 3,3 3,3 0,0 0),(1 1,2 1,2 2,1 2,1 1))";
  26
  27     bg::strategy::within::winding<Point> s;
  28
  29
  30     test_geometry<Point, polygon>("b1", "POINT(1 1)", box, s, true);
  31     test_geometry<Point, polygon>("b2", "POINT(3 3)", box, s, false);
  32
  33     // Test ALL corners (officialy false but some strategies might answer true)
  34     test_geometry<Point, polygon>("b3a", "POINT(0 0)", box, s, false);
  35     test_geometry<Point, polygon>("b3b", "POINT(0 2)", box, s, false);
  36     test_geometry<Point, polygon>("b3c", "POINT(2 2)", box, s, false);
  37     test_geometry<Point, polygon>("b3d", "POINT(2 0)", box, s, false);
  38
  39     // Test ALL sides (officialy false but some strategies might answer true)
  40     test_geometry<Point, polygon>("b4a", "POINT(0 1)", box, s, false);
  41     test_geometry<Point, polygon>("b4b", "POINT(1 2)", box, s, false);
  42     test_geometry<Point, polygon>("b4c", "POINT(2 1)", box, s, false);
  43     test_geometry<Point, polygon>("b4d", "POINT(1 0)", box, s, false);
  44
  45
  46     test_geometry<Point, polygon>("t1", "POINT(1 1)", triangle, s, true);
  47     test_geometry<Point, polygon>("t2", "POINT(3 3)", triangle, s, false);
  48
  49     test_geometry<Point, polygon>("t3a", "POINT(0 0)", triangle, s, false);
  50     test_geometry<Point, polygon>("t3b", "POINT(0 4)", triangle, s, false);
  51     test_geometry<Point, polygon>("t3c", "POINT(5 0)", triangle, s, false);
  52
  53     test_geometry<Point, polygon>("t4a", "POINT(0 2)", triangle, s, false);
  54     test_geometry<Point, polygon>("t4b", "POINT(3 2)", triangle, s, false);
  55     test_geometry<Point, polygon>("t4c", "POINT(2 0)", triangle, s, false);
  56
  57
  58     test_geometry<Point, polygon>("h1", "POINT(0.5 0.5)", with_hole, s, true);
  59     test_geometry<Point, polygon>("h2a", "POINT(1.5 1.5)", with_hole, s, false);
  60     test_geometry<Point, polygon>("h2b", "POINT(5 5)", with_hole, s, false);
  61
  62     test_geometry<Point, polygon>("h3a", "POINT(1 1)", with_hole, s, false);
  63     test_geometry<Point, polygon>("h3b", "POINT(2 2)", with_hole, s, false);
  64     test_geometry<Point, polygon>("h3c", "POINT(0 0)", with_hole, s, false);
  65
  66     test_geometry<Point, polygon>("h4a", "POINT(1 1.5)", with_hole, s, false);
  67     test_geometry<Point, polygon>("h4b", "POINT(1.5 2)", with_hole, s, false);
  68
  69     // Lying ON (one of the sides of) interior ring
  70     test_geometry<Point, polygon>("#77-1", "POINT(6 3.5)",
  71         "POLYGON((5 3,5 4,4 4,4 5,3 5,3 6,5 6,5 5,7 5,7 6,8 6,8 5,9 5,9 2,8 2,8 1,7 1,7 2,5 2,5 3),(6 3,8 3,8 4,6 4,6 3))",
  72         s, false);
  73 }
  74
  75 template <typename T>
  76 void test_spherical()
  77 {
  78     typedef bg::model::point<T, 2, bg::cs::spherical_equatorial<bg::degree> > point;
  79     typedef bg::model::polygon<point> polygon;
  80
  81     bg::strategy::within::winding<point> s;
  82
  83
  84     // Ticket #9354
  85     test_geometry<point, polygon>(
  86         "#9354",
  87         "POINT(-78.1239 25.9556)",
  88         "POLYGON((-97.08466667 25.95683333, -97.13683333 25.954, -97.1 26, -97.08466667 25.95683333))",
  89         s,
  90         false);
  91
  92     test_geometry<point, polygon>(
  93         "sph1N",
  94         "POINT(0 10.001)",
  95         "POLYGON((-10 10, 10 10, 10 -10, -10 -10, -10 10))",
  96         s,
  97         bg::strategy::side::spherical_side_formula<>::apply(
  98             point(-10, 10),
  99             point(10, 10),
 100             point(0, (T)10.001)) == -1 // right side
 101         /*true*/);
 102     test_geometry<point, polygon>(
 103         "sph1S",
 104         "POINT(0 -10.001)",
 105         "POLYGON((-10 10, 10 10, 10 -10, -10 -10, -10 10))",
 106         s,
 107         bg::strategy::side::spherical_side_formula<>::apply(
 108               point(10, -10),
 109               point(-10, -10),
 110               point(0, (T)-10.001)) == -1 // right side
 111       /*true*/);
 112
 113     test_geometry<point, polygon>(
 114         "sph2S",
 115         "POINT(0 10.001)",
 116         "POLYGON((-10 20, 10 20, 10 10, -10 10, -10 20))",
 117         s,
 118         bg::strategy::side::spherical_side_formula<>::apply(
 119               point(10, 10),
 120               point(-10, 10),
 121               point(0, (T)10.001)) == -1 // right side
 122         /*false*/);
 123
 124     test_geometry<point, polygon>(
 125         "sph3N",
 126         "POINT(0 10)",
 127         "POLYGON((-10 10, 10 10, 10 -10, -10 -10, -10 10))",
 128         s,
 129         bg::strategy::side::spherical_side_formula<>::apply(
 130             point(-10, 10),
 131             point(10, 10),
 132             point(0, (T)10.001)) == -1 // right side
 133         /*true*/);
 134     test_geometry<point, polygon>(
 135         "sph3S",
 136         "POINT(0 -10)",
 137         "POLYGON((-10 10, 10 10, 10 -10, -10 -10, -10 10))",
 138         s,
 139         bg::strategy::side::spherical_side_formula<>::apply(
 140               point(10, -10),
 141               point(-10, -10),
 142               point(0, (T)-10.001)) == -1 // right side
 143       /*true*/);
 144
 145     test_geometry<point, polygon>(
 146         "sphEq1",
 147         "POINT(179 10)",
 148         "POLYGON((170 10, -170 10, -170 0, 170 0, 170 10))",
 149         s,
 150         true,
 151         false);
 152     test_geometry<point, polygon>(
 153         "sphEq2",
 154         "POINT(179 10)",
 155         "POLYGON((170 20, -170 20, -170 10, 170 10, 170 20))",
 156         s,
 157         false,
 158         false);
 159     test_geometry<point, polygon>(
 160         "sphEq3",
 161         "POINT(-179 10)",
 162         "POLYGON((170 10, -170 10, -170 0, 170 0, 170 10))",
 163         s,
 164         true,
 165         false);
 166     test_geometry<point, polygon>(
 167         "sphEq4",
 168         "POINT(-179 10)",
 169         "POLYGON((170 20, -170 20, -170 10, 170 10, 170 20))",
 170         s,
 171         false,
 172         false);
 173
 174     test_geometry<point, polygon>(
 175         "sphEq5",
 176         "POINT(169 10)",
 177         "POLYGON((170 20, -170 20, -170 10, 170 10, 170 20))",
 178         s,
 179         false,
 180         false);
 181     test_geometry<point, polygon>(
 182         "sphEq6",
 183         "POINT(-169 10)",
 184         "POLYGON((170 20, -170 20, -170 10, 170 10, 170 20))",
 185         s,
 186         false,
 187         false);
 188     test_geometry<point, polygon>(
 189         "sphEq7",
 190         "POINT(169 10)",
 191         "POLYGON((170 10, -170 10, -170 0, 170 0, 170 10))",
 192         s,
 193         false,
 194         false);
 195     test_geometry<point, polygon>(
 196         "sphEq8",
 197         "POINT(-169 10)",
 198         "POLYGON((170 10, -170 10, -170 0, 170 0, 170 10))",
 199         s,
 200         false,
 201         false);
 202 }
 203
 204 int test_main(int, char* [])
 205 {
 206     test_cartesian<bg::model::point<float, 2, bg::cs::cartesian> >();
 207     test_cartesian<bg::model::point<double, 2, bg::cs::cartesian> >();
 208
 209     test_spherical<float>();
 210     test_spherical<double>();
 211
 212 #if defined(HAVE_TTMATH)
 213     test_cartesian<bg::model::point<ttmath_big, 2, bg::cs::cartesian> >();
 214     test_spherical<ttmath_big>();
 215 #endif
 216
 217     return 0;
 218 }