update ceph source to reef 18.1.2

[ceph.git] / ceph / src / boost / boost / geometry / strategies / cartesian / azimuth.hpp

diff --git a/ceph/src/boost/boost/geometry/strategies/cartesian/azimuth.hpp b/ceph/src/boost/boost/geometry/strategies/cartesian/azimuth.hpp
index 2ff16e94584a03eff4e0f035101a86d43bec129b..9035fdddcbe7e1d019916d77b9414a10eb42832f 100644 (file)
--- a/ceph/src/boost/boost/geometry/strategies/cartesian/azimuth.hpp
+++ b/ceph/src/boost/boost/geometry/strategies/cartesian/azimuth.hpp
@@ -1,6 +1,6 @@
 // Boost.Geometry (aka GGL, Generic Geometry Library)
 
-// Copyright (c) 2016-2018 Oracle and/or its affiliates.
+// Copyright (c) 2016-2021 Oracle and/or its affiliates.
 // Contributed and/or modified by Vissarion Fisikopoulos, on behalf of Oracle
 // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
 
@@ -11,32 +11,81 @@
 #ifndef BOOST_GEOMETRY_STRATEGIES_CARTESIAN_AZIMUTH_HPP
 #define BOOST_GEOMETRY_STRATEGIES_CARTESIAN_AZIMUTH_HPP
 
+#include <cmath>
+
 #include <boost/geometry/core/tags.hpp>
+#include <boost/geometry/core/coordinate_promotion.hpp>
 
 #include <boost/geometry/strategies/azimuth.hpp>
 
+#include <boost/geometry/util/select_most_precise.hpp>
+
 namespace boost { namespace geometry
 {
 
 namespace strategy { namespace azimuth
 {
 
-template
-<
-    typename CalculationType = void
->
+template <typename CalculationType = void>
 class cartesian
-{};
+{
+public:
+    template <typename T1, typename T2>
+    struct result_type
+        : geometry::select_most_precise
+              <
+                  // NOTE: this promotes any integer type to double
+                  typename geometry::promote_floating_point<T1, double>::type,
+                  typename geometry::promote_floating_point<T2, double>::type,
+                  CalculationType
+              >
+    {};
+
+    template <typename T1, typename T2, typename Result>
+    static inline void apply(T1 const& x1, T1 const& y1,
+                             T2 const& x2, T2 const& y2,
+                             Result& a1, Result& a2)
+    {
+        compute(x1, y1, x2, y2, a1, a2);
+    }
+    template <typename T1, typename T2, typename Result>
+    static inline void apply(T1 const& x1, T1 const& y1,
+                             T2 const& x2, T2 const& y2,
+                             Result& a1)
+    {
+        compute(x1, y1, x2, y2, a1, a1);
+    }
+    template <typename T1, typename T2, typename Result>
+    static inline void apply_reverse(T1 const& x1, T1 const& y1,
+                                     T2 const& x2, T2 const& y2,
+                                     Result& a2)
+    {
+        compute(x1, y1, x2, y2, a2, a2);
+    }
+
+private:
+    template <typename T1, typename T2, typename Result>
+    static inline void compute(T1 const& x1, T1 const& y1,
+                               T2 const& x2, T2 const& y2,
+                               Result& a1, Result& a2)
+    {
+        typedef typename result_type<T1, T2>::type calc_t;
+
+        // NOTE: azimuth 0 is at Y axis, increasing right
+        // as in spherical/geographic where 0 is at North axis
+        a1 = a2 = atan2(calc_t(x2) - calc_t(x1), calc_t(y2) - calc_t(y1));
+    }
+};
 
 #ifndef DOXYGEN_NO_STRATEGY_SPECIALIZATIONS
 
 namespace services
 {
 
-template <typename CalculationType>
-struct default_strategy<cartesian_tag, CalculationType>
+template <>
+struct default_strategy<cartesian_tag>
 {
-    typedef strategy::azimuth::cartesian<CalculationType> type;
+    typedef strategy::azimuth::cartesian<> type;
 };
 
 }