1 // Boost.Geometry
2 // Unit Test
3
4 // Copyright (c) 2018 Oracle and/or its affiliates.
5
6 // Contributed and/or modified by Vissarion Fysikopoulos, 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 #define BOOST_GEOMETRY_NORMALIZE_LATITUDE
13
14 #include <sstream>
15
16 #include "test_formula.hpp"
17 #include "direct_meridian_cases.hpp"
18
19 #include <boost/geometry/srs/srs.hpp>
20 #include <boost/geometry/formulas/vincenty_direct.hpp>
21 #include <boost/geometry/formulas/meridian_direct.hpp>
22
23 template <typename Result, typename Expected, typename Reference>
check_meridian_direct(Result & result,Expected const & expected,Reference & reference,double reference_error)24 void check_meridian_direct(Result& result,
25 Expected const& expected,
26 Reference& reference,
27 double reference_error)
28 {
29 boost::geometry::math::normalize_spheroidal_coordinates
30 <
31 boost::geometry::radian,
32 double
33 >(result.lon2, result.lat2);
34
35 boost::geometry::math::normalize_spheroidal_coordinates
36 <
37 boost::geometry::radian,
38 double
39 >(reference.lon2, reference.lat2);
40
41 std::stringstream ss;
42 ss << "(" << result.lon2 * bg::math::r2d<double>()
43 << " " << result.lat2 * bg::math::r2d<double>() << ")";
44
45 check_one("lon:" + ss.str(), result.lon2, expected.lon, reference.lon2,
46 reference_error);
47 check_one("lat:" + ss.str(), result.lat2, expected.lat, reference.lat2,
48 reference_error);
49 check_one("rev_az:" + ss.str(), result.reverse_azimuth,
50 result.reverse_azimuth, reference.reverse_azimuth, reference_error);
51 check_one("red len:" + ss.str(), result.reduced_length, result.reduced_length,
52 reference.reduced_length, 0.01);
53 check_one("geo scale:" + ss.str(), result.geodesic_scale, result.geodesic_scale,
54 reference.geodesic_scale, 0.01);
55
56 }
57
test_all(expected_results const & results)58 void test_all(expected_results const& results)
59 {
60 double const d2r = bg::math::d2r<double>();
61
62 double lon1_rad = results.p1.lon * d2r;
63 double lat1_rad = results.p1.lat * d2r;
64 coordinates expected_point;
65 expected_point.lon = results.p2.lon * d2r;
66 expected_point.lat = results.p2.lat * d2r;
67 double distance = results.distance;
68 bool direction = results.direction;
69
70 // WGS84
71 bg::srs::spheroid<double> spheroid(6378137.0, 6356752.3142451793);
72
73 bg::formula::result_direct<double> vincenty_result;
74 bg::formula::result_direct<double> meridian_result;
75
76 typedef bg::formula::vincenty_direct<double, true, true, true, true> vi_t;
77 double vincenty_azimuth = direction ? 0.0 : bg::math::pi<double>();
78 vincenty_result = vi_t::apply(lon1_rad, lat1_rad, distance, vincenty_azimuth, spheroid);
79
80 {
81 typedef bg::formula::meridian_direct<double, true, true, true, true, 1> eli;
82 meridian_result = eli::apply(lon1_rad, lat1_rad, distance, direction, spheroid);
83 check_meridian_direct(meridian_result, expected_point, vincenty_result, 0.001);
84 }
85 {
86 typedef bg::formula::meridian_direct<double, true, true, true, true, 2> eli;
87 meridian_result = eli::apply(lon1_rad, lat1_rad, distance, direction, spheroid);
88 check_meridian_direct(meridian_result, expected_point, vincenty_result, 0.00001);
89 }
90 {
91 typedef bg::formula::meridian_direct<double, true, true, true, true, 3> eli;
92 meridian_result = eli::apply(lon1_rad, lat1_rad, distance, direction, spheroid);
93 check_meridian_direct(meridian_result, expected_point, vincenty_result, 0.00000001);
94 }
95 {
96 typedef bg::formula::meridian_direct<double, true, true, true, true, 4> eli;
97 meridian_result = eli::apply(lon1_rad, lat1_rad, distance, direction, spheroid);
98 check_meridian_direct(meridian_result, expected_point, vincenty_result, 0.00000001);
99 }
100 {
101 typedef bg::formula::meridian_direct<double, true, true, true, true, 5> eli;
102 meridian_result = eli::apply(lon1_rad, lat1_rad, distance, direction, spheroid);
103 check_meridian_direct(meridian_result, expected_point, vincenty_result, 0.00000000001);
104 }
105 }
106
test_main(int,char * [])107 int test_main(int, char*[])
108 {
109 for (size_t i = 0; i < expected_size; ++i)
110 {
111 test_all(expected[i]);
112 }
113
114 return 0;
115 }
116