1 | // Boost.Geometry |
2 | // Unit Test |
3 | |
4 | // Copyright (c) 2022 Oracle and/or its affiliates. |
5 | |
6 | // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle |
7 | |
8 | // Use, modification and distribution is subject to the Boost Software License, |
9 | // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at |
10 | // http://www.boost.org/LICENSE_1_0.txt) |
11 | |
12 | |
13 | #include <limits> |
14 | #include <type_traits> |
15 | |
16 | #include <geometry_test_common.hpp> |
17 | |
18 | #include <boost/geometry/util/condition.hpp> |
19 | #include <boost/geometry/util/math.hpp> |
20 | #include <boost/geometry/util/normalize_spheroidal_box_coordinates.hpp> |
21 | #include <boost/geometry/util/normalize_spheroidal_coordinates.hpp> |
22 | |
23 | |
24 | namespace bgm = bg::math; |
25 | |
26 | template <typename T, std::enable_if_t<std::numeric_limits<T>::has_quiet_NaN, int> = 0> |
27 | void test_equals(T const& v, T const& r) |
28 | { |
29 | bool const check = bg::math::equals(v, r) || (std::isnan(v) && std::isnan(r)); |
30 | BOOST_CHECK_MESSAGE(check, v << " different than " << r); |
31 | } |
32 | |
33 | template <typename T, std::enable_if_t<! std::numeric_limits<T>::has_quiet_NaN, int> = 0> |
34 | void test_equals(T const& v, T const& r) |
35 | { |
36 | bool const check = bg::math::equals(v, r); |
37 | BOOST_CHECK_MESSAGE(check, v << " different than " << r); |
38 | } |
39 | |
40 | template <typename T> |
41 | void test_one(T lon, T lat, T const& res_lon, T const& res_lat) |
42 | { |
43 | bg::math::normalize_spheroidal_coordinates<bg::degree>(lon, lat); |
44 | test_equals(lon, res_lon); |
45 | test_equals(lat, res_lat); |
46 | } |
47 | |
48 | template <typename T> |
49 | void test_box(T lon1, T lat1, T lon2, T lat2, |
50 | T const& res_lon1, T const& res_lat1, T const& res_lon2, T const& res_lat2) |
51 | { |
52 | bg::math::normalize_spheroidal_box_coordinates<bg::degree>(lon1, lat1, lon2, lat2); |
53 | test_equals(lon1, res_lon1); |
54 | test_equals(lat1, res_lat1); |
55 | test_equals(lon2, res_lon2); |
56 | test_equals(lat2, res_lat2); |
57 | } |
58 | |
59 | template <typename T> |
60 | void test_all() |
61 | { |
62 | test_one<T>(270, 0, -90, 0); |
63 | |
64 | test_one<T>(-180, 10, 180, 10); |
65 | test_one<T>(-180, -90, 0, -90); |
66 | |
67 | if (BOOST_GEOMETRY_CONDITION(std::numeric_limits<T>::has_quiet_NaN)) |
68 | { |
69 | T nan = std::numeric_limits<T>::quiet_NaN(); |
70 | test_one<T>(nan, 10, nan, 10); |
71 | test_one<T>(10, nan, 10, nan); |
72 | test_one<T>(nan, nan, nan, nan); |
73 | |
74 | test_box<T>(nan, 10, 20, 20, nan, 10, 20, 20); |
75 | test_box<T>(10, nan, 20, 20, 10, nan, 20, 20); |
76 | test_box<T>(10, 10, nan, 20, 10, 10, nan, 20); |
77 | test_box<T>(10, 10, 20, nan, 10, 10, 20, nan); |
78 | |
79 | test_box<T>(nan, 10, 270, 20, nan, 10, -90, 20); |
80 | test_box<T>(270, 10, nan, 20, -90, 10, nan, 20); |
81 | } |
82 | } |
83 | |
84 | int test_main(int, char* []) |
85 | { |
86 | test_all<int>(); |
87 | test_all<float>(); |
88 | test_all<double>(); |
89 | |
90 | return 0; |
91 | } |
92 | |