libs/geometry/test/geometries/infinite_line.cpp

   1 // Boost.Geometry
   2 // Unit Test
   3
   4 // Copyright (c) 2019 Barend Gehrels, Amsterdam, the Netherlands.
   5
   6 // This file was modified by Oracle on 2019.
   7 // Modifications copyright (c) 2019, 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/algorithms/detail/make/make.hpp>
  17 #include <boost/geometry/geometries/infinite_line.hpp>
  18 #include <boost/geometry/geometries/point.hpp>
  19 #include <boost/geometry/util/math.hpp>
  20
  21 namespace
  22 {
  23     // Boost.Test does not support BOOST_CHECK_CLOSE for integral types
  24     template <typename T>
  25     bool is_small(T const& value)
  26     {
  27         static long double const epsilon = 1.0e-5;
  28         return bg::math::abs(value) < epsilon;
  29     }
  30 }
  31
  32 template <typename T, typename C>
  33 void verify_point_on_line(bg::model::infinite_line<T> const& line,
  34                           C const& x, C const& y)
  35 {
  36     BOOST_CHECK_MESSAGE(is_small(line.a * x + line.b * y + line.c),
  37                         "Point is not located on the line");
  38 }
  39
  40 template <typename T>
  41 void test_make()
  42 {
  43     typedef bg::model::infinite_line<T> line_type;
  44
  45     // Horizontal through origin
  46     line_type line = bg::detail::make::make_infinite_line<T>(0, 0, 10, 0);
  47     verify_point_on_line(line, 0, 0);
  48     verify_point_on_line(line, 10, 0);
  49
  50     // Horizontal line above origin
  51     line = bg::detail::make::make_infinite_line<T>(0, 5, 10, 5);
  52     verify_point_on_line(line, 0, 5);
  53     verify_point_on_line(line, 10, 5);
  54
  55     // Vertical through origin
  56     line = bg::detail::make::make_infinite_line<T>(0, 0, 0, 10);
  57     verify_point_on_line(line, 0, 0);
  58     verify_point_on_line(line, 0, 10);
  59
  60     // Vertical line left from origin
  61     line = bg::detail::make::make_infinite_line<T>(5, 0, 5, 10);
  62     verify_point_on_line(line, 5, 0);
  63     verify_point_on_line(line, 5, 10);
  64
  65     // Diagonal through origin
  66     line = bg::detail::make::make_infinite_line<T>(0, 0, 8, 10);
  67     verify_point_on_line(line, 0, 0);
  68     verify_point_on_line(line, 8, 10);
  69
  70     // Diagonal not through origin
  71     line = bg::detail::make::make_infinite_line<T>(5, 2, -8, 10);
  72     verify_point_on_line(line, 5, 2);
  73     verify_point_on_line(line, -8, 10);
  74 }
  75
  76
  77 template <typename T>
  78 void test_all()
  79 {
  80     test_make<T>();
  81 }
  82
  83 int test_main(int, char* [])
  84 {
  85     test_all<double>();
  86     test_all<long double>();
  87     test_all<float>();
  88     test_all<int>();
  89     return 0;
  90 }