1 // Boost.Geometry (aka GGL, Generic Geometry Library)
2 // Unit Test
3
4 // Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands.
5 // Copyright (c) 2008-2012 Bruno Lalande, Paris, France.
6 // Copyright (c) 2009-2012 Mateusz Loskot, London, UK.
7
8 // Parts of Boost.Geometry are redesigned from Geodan's Geographic Library
9 // (geolib/GGL), copyright (c) 1995-2010 Geodan, Amsterdam, the Netherlands.
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 #include <boost/geometry/core/radian_access.hpp>
19
20 #include <boost/geometry/core/cs.hpp>
21 #include <boost/geometry/core/coordinate_type.hpp>
22 #include <boost/geometry/geometries/point.hpp>
23
24
25 template <std::size_t D, typename P, typename T>
test_get(T const & c,T const & e)26 void test_get(T const& c, T const& e)
27 {
28 //std::cout << "get_as_radian " << typeid(P).name() << std::endl;
29
30 typedef typename bg::coordinate_type<P>::type coordinate_type;
31
32 P p;
33 bg::set<D>(p, coordinate_type(c));
34
35 coordinate_type g = bg::get_as_radian<D>(p);
36
37 BOOST_CHECK_CLOSE(double(g), double(e), 0.0001);
38 }
39
40
41 template <std::size_t D, typename P, typename T>
test_set(T const & c,T const & e)42 void test_set(T const& c, T const& e)
43 {
44 //std::cout << "set_from_radian " << typeid(P).name() << std::endl;
45
46 typedef typename bg::coordinate_type<P>::type coordinate_type;
47
48 P p;
49 bg::set_from_radian<D>(p, coordinate_type(c));
50
51 coordinate_type g = bg::get<D>(p);
52
53 BOOST_CHECK_CLOSE(double(g), double(e), 0.0001);
54 }
55
56
57 template <typename T>
test()58 void test()
59 {
60 double d2r = 3.1415926535897932384626433832795 / 180.0;
61
62 // Degree
63 test_get<0, bg::model::point<T, 2, bg::cs::spherical<bg::degree> > >
64 (1.0, 1.0 * d2r);
65 test_get<1, bg::model::point<T, 2, bg::cs::spherical<bg::degree> > >
66 (2.0, 2.0 * d2r);
67
68 test_set<0, bg::model::point<T, 2, bg::cs::spherical<bg::degree> > >
69 (1.0, 1.0 / d2r);
70 test_set<1, bg::model::point<T, 2, bg::cs::spherical<bg::degree> > >
71 (2.0, 2.0 / d2r);
72
73
74 // Radian
75 test_get<0, bg::model::point<T, 2, bg::cs::spherical<bg::radian> > >
76 (1.0, 1.0);
77 test_get<1, bg::model::point<T, 2, bg::cs::spherical<bg::radian> > >
78 (1.0, 1.0);
79
80 test_set<0, bg::model::point<T, 2, bg::cs::spherical<bg::radian> > >
81 (1.0, 1.0);
82 test_set<1, bg::model::point<T, 2, bg::cs::spherical<bg::radian> > >
83 (1.0, 1.0);
84
85
86 // Cartesian (== untouched, no conversion)
87 test_get<0, bg::model::point<T, 2, bg::cs::cartesian> >
88 (1.0, 1.0);
89 test_get<1, bg::model::point<T, 2, bg::cs::cartesian> >
90 (1.0, 1.0);
91
92 // Dimension >=2, should always be untouched, even for degree
93 // (assuming lat/lon + height)
94 test_set<2, bg::model::point<T, 3, bg::cs::spherical<bg::radian> > >
95 (1.0, 1.0);
96 test_set<2, bg::model::point<T, 3, bg::cs::spherical<bg::degree> > >
97 (1.0, 1.0);
98
99 }
100
101
102
103
test_main(int,char * [])104 int test_main(int, char* [])
105 {
106 //test<float>();
107 test<double>();
108
109 return 0;
110 }
111