ceph/src/boost/libs/geometry/test/arithmetic/infinite_line_functions.cpp

   1 // Boost.Geometry (aka GGL, Generic Geometry Library)
   2 // Unit Test
   3
   4 // Copyright (c) 2018-2019 Barend Gehrels, Amsterdam, the Netherlands.
   5
   6 // This file was modified by Oracle on 2019, 2020.
   7 // Modifications copyright (c) 2019, 2020, 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 #include <geometry_test_common.hpp>
  15
  16 #include <boost/geometry/arithmetic/infinite_line_functions.hpp>
  17 #include <boost/geometry/geometries/infinite_line.hpp>
  18 #include <boost/geometry/algorithms/detail/make/make.hpp>
  19 #include <boost/geometry/geometries/point.hpp>
  20 #include <boost/geometry/io/wkt/wkt.hpp>
  21
  22 namespace
  23 {
  24     // Boost.Test does not support BOOST_CHECK_CLOSE for integral types
  25     template <typename T>
  26     bool is_small(T const& value)
  27     {
  28         static long double const epsilon = 1.0e-5;
  29         return bg::math::abs(value) < epsilon;
  30     }
  31 }
  32
  33 template <typename T, typename C>
  34 void verify_point_on_line(bg::model::infinite_line<T> const& line,
  35                           C const& x, C const& y)
  36 {
  37     BOOST_CHECK_MESSAGE(is_small(line.a * x + line.b * y + line.c),
  38                         "Point is not located on the line");
  39 }
  40
  41 template <typename T>
  42 void test_side_value()
  43 {
  44     typedef bg::model::infinite_line<T> line_type;
  45
  46     // Horizontal line going right
  47     line_type line = bg::detail::make::make_infinite_line<T>(0, 0, 10, 0);
  48
  49     // Point above (= on left side)
  50     T d = bg::arithmetic::side_value(line, 5, 5);
  51     BOOST_CHECK_MESSAGE(d > 0, "point not on left side");
  52
  53     // Point below (= on right side)
  54     d = bg::arithmetic::side_value(line, 5, -5);
  55     BOOST_CHECK_MESSAGE(d < 0, "point not on right side");
  56
  57     // Diagonal not through origin, from right (down) to left (up)
  58     line = bg::detail::make::make_infinite_line<T>(5, 2, -7, 10);
  59     d = bg::arithmetic::side_value(line, 5, 2);
  60     BOOST_CHECK_MESSAGE(is_small(d), "point not on line");
  61     d = bg::arithmetic::side_value(line, -7, 10);
  62     BOOST_CHECK_MESSAGE(is_small(d), "point not on line");
  63
  64     // vector is (-12, 8), move (-3,2) on the line from (5,2)
  65     d = bg::arithmetic::side_value(line, 2, 4);
  66     BOOST_CHECK_MESSAGE(is_small(d), "point not on line");
  67
  68     // Go perpendicular (2,3) from (2,4) up, so right of the line (negative)
  69     d = bg::arithmetic::side_value(line, 4, 7);
  70     BOOST_CHECK_MESSAGE(d < 0, "point not on right side");
  71
  72     // Go perpendicular (2,3) from (2,4) down, so left of the line (positive)
  73     d = bg::arithmetic::side_value(line, 0, 1);
  74     BOOST_CHECK_MESSAGE(d > 0, "point not on left side");
  75 }
  76
  77
  78 template <typename T>
  79 void test_get_intersection()
  80 {
  81     typedef bg::model::infinite_line<T> line_type;
  82
  83     // Diagonal lines (first is same as in distance measure,
  84     // second is perpendicular and used there for distance measures)
  85     line_type p = bg::detail::make::make_infinite_line<T>(5, 2, -7, 10);
  86     line_type q = bg::detail::make::make_infinite_line<T>(4, 7, 0, 1);
  87
  88     typedef bg::model::point<T, 2, bg::cs::cartesian> point_type;
  89     point_type ip;
  90     BOOST_CHECK(bg::arithmetic::intersection_point(p, q, ip));
  91
  92     BOOST_CHECK_MESSAGE(is_small(bg::get<0>(ip) - 2), "x-coordinate wrong");
  93     BOOST_CHECK_MESSAGE(is_small(bg::get<1>(ip) - 4), "y-coordinate wrong");
  94
  95     verify_point_on_line(p, bg::get<0>(ip), bg::get<1>(ip));
  96     verify_point_on_line(q, bg::get<0>(ip), bg::get<1>(ip));
  97 }
  98
  99 template <typename T>
 100 void test_degenerate()
 101 {
 102     typedef bg::model::infinite_line<T> line_type;
 103
 104     line_type line = bg::detail::make::make_infinite_line<T>(0, 0, 10, 0);
 105     BOOST_CHECK(! bg::arithmetic::is_degenerate(line));
 106
 107     line = bg::detail::make::make_infinite_line<T>(0, 0, 0, 10);
 108     BOOST_CHECK(! bg::arithmetic::is_degenerate(line));
 109
 110     line = bg::detail::make::make_infinite_line<T>(0, 0, 10, 10);
 111     BOOST_CHECK(! bg::arithmetic::is_degenerate(line));
 112
 113     line = bg::detail::make::make_infinite_line<T>(0, 0, 0, 0);
 114     BOOST_CHECK(bg::arithmetic::is_degenerate(line));
 115 }
 116
 117
 118 template <typename T>
 119 void test_all()
 120 {
 121     test_side_value<T>();
 122     test_get_intersection<T>();
 123     test_degenerate<T>();
 124 }
 125
 126 int test_main(int, char* [])
 127 {
 128     test_all<double>();
 129     test_all<long double>();
 130     test_all<float>();
 131     test_all<int>();
 132     return 0;
 133 }