]>
Commit | Line | Data |
---|---|---|
11fdf7f2 TL |
1 | // Boost.Geometry (aka GGL, Generic Geometry Library) |
2 | // This file is manually converted from PROJ4 | |
3 | ||
4 | // Copyright (c) 2008-2012 Barend Gehrels, Amsterdam, the Netherlands. | |
5 | ||
92f5a8d4 TL |
6 | // This file was modified by Oracle on 2018. |
7 | // Modifications copyright (c) 2018, Oracle and/or its affiliates. | |
8 | // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle | |
9 | ||
11fdf7f2 TL |
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 | // This file is converted from PROJ4, http://trac.osgeo.org/proj | |
15 | // PROJ4 is originally written by Gerald Evenden (then of the USGS) | |
16 | // PROJ4 is maintained by Frank Warmerdam | |
17 | // PROJ4 is converted to Geometry Library by Barend Gehrels (Geodan, Amsterdam) | |
18 | ||
19 | // Original copyright notice: | |
20 | ||
21 | // Permission is hereby granted, free of charge, to any person obtaining a | |
22 | // copy of this software and associated documentation files (the "Software"), | |
23 | // to deal in the Software without restriction, including without limitation | |
24 | // the rights to use, copy, modify, merge, publish, distribute, sublicense, | |
25 | // and/or sell copies of the Software, and to permit persons to whom the | |
26 | // Software is furnished to do so, subject to the following conditions: | |
27 | ||
28 | // The above copyright notice and this permission notice shall be included | |
29 | // in all copies or substantial portions of the Software. | |
30 | ||
31 | // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS | |
32 | // OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, | |
33 | // FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL | |
34 | // THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER | |
35 | // LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING | |
36 | // FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER | |
37 | // DEALINGS IN THE SOFTWARE. | |
38 | ||
39 | #ifndef BOOST_GEOMETRY_PROJECTIONS_IMPL_AASINCOS_HPP | |
40 | #define BOOST_GEOMETRY_PROJECTIONS_IMPL_AASINCOS_HPP | |
41 | ||
42 | ||
43 | #include <cmath> | |
44 | ||
92f5a8d4 TL |
45 | #include <boost/geometry/srs/projections/exception.hpp> |
46 | #include <boost/geometry/srs/projections/impl/pj_strerrno.hpp> | |
11fdf7f2 TL |
47 | #include <boost/geometry/util/math.hpp> |
48 | ||
49 | ||
50 | namespace boost { namespace geometry { namespace projections | |
51 | { | |
52 | ||
53 | namespace detail | |
54 | { | |
55 | ||
56 | namespace aasincos | |
57 | { | |
58 | template <typename T> | |
59 | inline T ONE_TOL() { return 1.00000000000001; } | |
60 | //template <typename T> | |
61 | //inline T TOL() { return 0.000000001; } | |
62 | template <typename T> | |
63 | inline T ATOL() { return 1e-50; } | |
64 | } | |
65 | ||
66 | template <typename T> | |
67 | inline T aasin(T const& v) | |
68 | { | |
69 | T av = 0; | |
70 | ||
71 | if ((av = geometry::math::abs(v)) >= 1.0) | |
72 | { | |
73 | if (av > aasincos::ONE_TOL<T>()) | |
74 | { | |
92f5a8d4 | 75 | BOOST_THROW_EXCEPTION( projection_exception(error_acos_asin_arg_too_large) ); |
11fdf7f2 TL |
76 | } |
77 | return (v < 0.0 ? -geometry::math::half_pi<T>() : geometry::math::half_pi<T>()); | |
78 | } | |
79 | ||
80 | return asin(v); | |
81 | } | |
82 | ||
83 | template <typename T> | |
84 | inline T aacos(T const& v) | |
85 | { | |
86 | T av = 0; | |
87 | ||
88 | if ((av = geometry::math::abs(v)) >= 1.0) | |
89 | { | |
90 | if (av > aasincos::ONE_TOL<T>()) | |
91 | { | |
92f5a8d4 | 92 | BOOST_THROW_EXCEPTION( projection_exception(error_acos_asin_arg_too_large) ); |
11fdf7f2 TL |
93 | } |
94 | return (v < 0.0 ? geometry::math::pi<T>() : 0.0); | |
95 | } | |
96 | ||
97 | return acos(v); | |
98 | } | |
99 | ||
100 | template <typename T> | |
101 | inline T asqrt(T const& v) | |
102 | { | |
103 | return ((v <= 0) ? 0 : sqrt(v)); | |
104 | } | |
105 | ||
106 | template <typename T> | |
107 | inline T aatan2(T const& n, T const& d) | |
108 | { | |
109 | return ((geometry::math::abs(n) < aasincos::ATOL<T>() | |
110 | && geometry::math::abs(d) < aasincos::ATOL<T>()) ? 0.0 : atan2(n, d)); | |
111 | } | |
112 | ||
113 | ||
114 | } // namespace detail | |
115 | ||
116 | ||
117 | }}} // namespace boost::geometry::projections | |
118 | ||
119 | ||
120 | #endif // BOOST_GEOMETRY_PROJECTIONS_IMPL_AASINCOS_HPP |