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

   1 // Boost.Geometry (aka GGL, Generic Geometry Library)
   2 // Unit Test
   3
   4 // Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands.
   5
   6 // This file was modified by Oracle on 2014, 2015.
   7 // Modifications copyright (c) 2014-2015, Oracle and/or its affiliates.
   8
   9 // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
  10
  11 // Use, modification and distribution is subject to the Boost Software License,
  12 // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
  13 // http://www.boost.org/LICENSE_1_0.txt)
  14
  15
  16 #include <geometry_test_common.hpp>
  17
  18
  19 #include <boost/geometry/algorithms/assign.hpp>
  20
  21 #include <boost/geometry/core/coordinate_type.hpp>
  22 #include <boost/geometry/core/cs.hpp>
  23
  24 #include <boost/geometry/geometries/point.hpp>
  25 #include <boost/geometry/geometries/segment.hpp>
  26
  27 #include <boost/geometry/strategies/spherical/side_by_cross_track.hpp>
  28 #include <boost/geometry/strategies/spherical/ssf.hpp>
  29 #include <boost/geometry/strategy/cartesian/side_robust.hpp>
  30 #include <boost/geometry/strategy/cartesian/side_by_triangle.hpp>
  31
  32 #include <boost/geometry/strategies/geographic/mapping_ssf.hpp>
  33 #include <boost/geometry/strategies/geographic/side_andoyer.hpp>
  34 #include <boost/geometry/strategies/geographic/side_thomas.hpp>
  35 #include <boost/geometry/strategies/geographic/side_vincenty.hpp>
  36
  37 #include <boost/geometry/util/math.hpp>
  38
  39
  40 namespace boost { namespace geometry {
  41
  42 template <typename Vector, typename Point1, typename Point2>
  43 static inline Vector create_vector(Point1 const& p1, Point2 const& p2)
  44 {
  45     Vector v;
  46     convert(p1, v);
  47     subtract_point(v, p2);
  48     return v;
  49 }
  50
  51 }}
  52
  53 inline char side_char(int side)
  54 {
  55     return side == 1 ? 'L'
  56         : side == -1 ? 'R'
  57         : '-'
  58         ;
  59 }
  60
  61 template <typename Point>
  62 void test_side1(std::string const& /*case_id*/, Point const& p1, Point const& p2, Point const& p3,
  63                    int expected, int expected_cartesian)
  64 {
  65     namespace bgss = bg::strategy::side;
  66
  67     // std::cout << case_id << ": ";
  68     //int s = bgss::side_via_plane<>::apply(p1, p2, p3);
  69     int side_ssf = bgss::spherical_side_formula<>::apply(p1, p2, p3);
  70     //int side2 = bgss::side_via_plane<>::apply(p1, p2, p3);
  71     int side_ct = bgss::side_by_cross_track<>::apply(p1, p2, p3);
  72
  73     // non-official
  74     typedef bg::srs::spheroid<double> spheroid;
  75     spheroid const sph(1.0, 1.0);
  76     int side_mssf1 = bgss::mapping_spherical_side_formula<spheroid>(sph).apply(p1, p2, p3);
  77     int side_mssf2 = bgss::mapping_spherical_side_formula<spheroid, bgss::mapping_reduced>(sph).apply(p1, p2, p3);
  78     int side_mssf3 = bgss::mapping_spherical_side_formula<spheroid, bgss::mapping_geocentric>(sph).apply(p1, p2, p3);
  79     int side_andoyer = bgss::andoyer<spheroid>(sph).apply(p1, p2, p3);
  80     int side_thomas = bgss::thomas<spheroid>(sph).apply(p1, p2, p3);
  81     int side_vincenty = bgss::vincenty<spheroid>(sph).apply(p1, p2, p3);
  82
  83     // cartesian
  84     typedef bg::strategy::side::services::default_strategy<bg::cartesian_tag>::type cartesian_strategy;
  85     int side_cart = cartesian_strategy::apply(p1, p2, p3);
  86
  87     BOOST_CHECK_EQUAL(side_ssf, expected);
  88     BOOST_CHECK_EQUAL(side_ct, expected);
  89     BOOST_CHECK_EQUAL(side_mssf1, expected);
  90     BOOST_CHECK_EQUAL(side_mssf2, expected);
  91     BOOST_CHECK_EQUAL(side_mssf3, expected);
  92     BOOST_CHECK_EQUAL(side_andoyer, expected);
  93     BOOST_CHECK_EQUAL(side_thomas, expected);
  94     BOOST_CHECK_EQUAL(side_vincenty, expected);
  95     BOOST_CHECK_EQUAL(side_cart, expected_cartesian);
  96
  97     /*
  98     std::cout
  99         << "exp: " << side_char(expected)
 100         << " ssf: " << side_char(side1)
 101         << " pln: " << side_char(side2)
 102         << " ct: " << side_char(side3)
 103         //<< " def: " << side_char(side4)
 104         << " cart: " << side_char(side5)
 105         << std::endl;
 106     */
 107 }
 108
 109 template <typename Point>
 110 void test_side(std::string const& case_id, Point const& p1, Point const& p2, Point const& p3,
 111                    int expected, int expected_cartesian = -999)
 112 {
 113     if (expected_cartesian == -999)
 114     {
 115         expected_cartesian = expected;
 116     }
 117     test_side1(case_id, p1, p2, p3, expected, expected_cartesian);
 118     test_side1(case_id, p2, p1, p3, -expected, -expected_cartesian);
 119 }
 120
 121
 122 template <typename Point>
 123 void test_all()
 124 {
 125     typedef typename bg::coordinate_type<Point>::type CT;
 126
 127     Point amsterdam(bg::math::rounding_cast<CT>(5.9), bg::math::rounding_cast<CT>(52.4));
 128     Point barcelona(2.0, 41.0);
 129     Point paris(2.0, 48.0);
 130     Point milan(7.0, 45.0);
 131
 132     //goto wrong;
 133
 134     test_side<Point>("bp-m", barcelona, paris, milan, -1);
 135     test_side<Point>("bm-p", barcelona, milan, paris, 1);
 136     test_side<Point>("mp-b", milan, paris, barcelona, 1);
 137
 138     test_side<Point>("am-p", amsterdam, milan, paris, -1);
 139     test_side<Point>("pm-a", paris, milan, amsterdam, 1);
 140
 141     // http://www.gcmap.com/mapui?P=30N+10E-50N+50E,39N+30E
 142     Point gcmap_p1(10.0, 30.0);
 143     Point gcmap_p2(50.0, 50.0);
 144     test_side<Point>("blog1", gcmap_p1, gcmap_p2, Point(30.0, 41.0), -1, 1);
 145     test_side<Point>("blog1", gcmap_p1, gcmap_p2, Point(30.0, 42.0), -1, 1);
 146     test_side<Point>("blog1", gcmap_p1, gcmap_p2, Point(30.0, 43.0), -1, 1);
 147     test_side<Point>("blog1", gcmap_p1, gcmap_p2, Point(30.0, 44.0), 1);
 148
 149     // http://www.gcmap.com/mapui?P=50N+80E-60N+50W,65N+30E
 150     Point gcmap_np1(80.0, 50.0);
 151     Point gcmap_np2(-50.0, 60.0);
 152     // http://www.gcmap.com/mapui?P=50N+140E-60N+10E,65N+30E
 153     //Point gcmap_np1(140.0, 50.0);
 154     //Point gcmap_np2(10.0, 60.0);
 155     //test_side<Point>(gcmap_np1, gcmap_np2, gcmap_np, 1);
 156     test_side<Point>("40", gcmap_np1, gcmap_np2, Point(30.0, 60.0), 1, -1);
 157     test_side<Point>("45", gcmap_np1, gcmap_np2, Point(30.0, 65.0), 1, -1);
 158     test_side<Point>("70", gcmap_np1, gcmap_np2, Point(30.0, 70.0), 1, -1);
 159     test_side<Point>("75", gcmap_np1, gcmap_np2, Point(30.0, 75.0), -1);
 160 }
 161
 162 int test_main(int, char* [])
 163 {
 164     test_all<bg::model::point<int, 2, bg::cs::spherical<bg::degree> > >();
 165     test_all<bg::model::point<double, 2, bg::cs::spherical_equatorial<bg::degree> > >();
 166
 167     return 0;
 168 }