// Boost.Geometry (aka GGL, Generic Geometry Library)
|
|
// Copyright (c) 2016-2020 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
|
|
// Use, modification and distribution is subject to the Boost Software License,
|
// Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
|
// http://www.boost.org/LICENSE_1_0.txt)
|
|
#ifndef BOOST_GEOMETRY_STRATEGIES_SPHERICAL_AZIMUTH_HPP
|
#define BOOST_GEOMETRY_STRATEGIES_SPHERICAL_AZIMUTH_HPP
|
|
|
#include <type_traits>
|
|
#include <boost/geometry/strategies/azimuth.hpp>
|
#include <boost/geometry/formulas/spherical.hpp>
|
|
|
namespace boost { namespace geometry
|
{
|
|
namespace strategy { namespace azimuth
|
{
|
|
template
|
<
|
typename CalculationType = void
|
>
|
class spherical
|
{
|
public :
|
|
inline spherical()
|
{}
|
|
template <typename T>
|
inline void apply(T const& lon1_rad, T const& lat1_rad,
|
T const& lon2_rad, T const& lat2_rad,
|
T& a1, T& a2) const
|
{
|
compute<true, true>(lon1_rad, lat1_rad,
|
lon2_rad, lat2_rad,
|
a1, a2);
|
}
|
template <typename T>
|
inline void apply(T const& lon1_rad, T const& lat1_rad,
|
T const& lon2_rad, T const& lat2_rad,
|
T& a1) const
|
{
|
compute<true, false>(lon1_rad, lat1_rad,
|
lon2_rad, lat2_rad,
|
a1, a1);
|
}
|
template <typename T>
|
inline void apply_reverse(T const& lon1_rad, T const& lat1_rad,
|
T const& lon2_rad, T const& lat2_rad,
|
T& a2) const
|
{
|
compute<false, true>(lon1_rad, lat1_rad,
|
lon2_rad, lat2_rad,
|
a2, a2);
|
}
|
|
private :
|
|
template
|
<
|
bool EnableAzimuth,
|
bool EnableReverseAzimuth,
|
typename T
|
>
|
inline void compute(T const& lon1_rad, T const& lat1_rad,
|
T const& lon2_rad, T const& lat2_rad,
|
T& a1, T& a2) const
|
{
|
typedef std::conditional_t
|
<
|
std::is_void<CalculationType>::value, T, CalculationType
|
> calc_t;
|
|
geometry::formula::result_spherical<calc_t>
|
result = geometry::formula::spherical_azimuth
|
<
|
calc_t,
|
EnableReverseAzimuth
|
>(calc_t(lon1_rad), calc_t(lat1_rad),
|
calc_t(lon2_rad), calc_t(lat2_rad));
|
|
if (EnableAzimuth)
|
{
|
a1 = result.azimuth;
|
}
|
if (EnableReverseAzimuth)
|
{
|
a2 = result.reverse_azimuth;
|
}
|
}
|
};
|
|
#ifndef DOXYGEN_NO_STRATEGY_SPECIALIZATIONS
|
|
namespace services
|
{
|
|
template <typename CalculationType>
|
struct default_strategy<spherical_equatorial_tag, CalculationType>
|
{
|
typedef strategy::azimuth::spherical<CalculationType> type;
|
};
|
|
/*
|
template <typename CalculationType>
|
struct default_strategy<spherical_polar_tag, CalculationType>
|
{
|
typedef strategy::azimuth::spherical<CalculationType> type;
|
};
|
*/
|
}
|
|
#endif // DOXYGEN_NO_STRATEGY_SPECIALIZATIONS
|
|
}} // namespace strategy::azimuth
|
|
|
}} // namespace boost::geometry
|
|
#endif // BOOST_GEOMETRY_STRATEGIES_SPHERICAL_AZIMUTH_HPP
|
