FM / formulas

/usr/include/boost/geometry/formulas

Name	Size	Mode	Actions
andoyer_inverse.hpp	9909	0644	edit dl rm
area_formulas.hpp	21897	0644	edit dl rm
authalic_radius_sqr.hpp	2654	0644	edit dl rm
differential_quantities.hpp	10976	0644	edit dl rm
eccentricity_sqr.hpp	1959	0644	edit dl rm
flattening.hpp	1847	0644	edit dl rm
geographic.hpp	14251	0644	edit dl rm
gnomonic_intersection.hpp	4905	0644	edit dl rm
gnomonic_spheroid.hpp	3989	0644	edit dl rm
interpolate_point_spherical.hpp	3586	0644	edit dl rm
karney_direct.hpp	10063	0644	edit dl rm
karney_inverse.hpp	35130	0644	edit dl rm
mean_radius.hpp	1894	0644	edit dl rm
meridian_direct.hpp	5139	0644	edit dl rm
meridian_inverse.hpp	4782	0644	edit dl rm
meridian_segment.hpp	2009	0644	edit dl rm
quarter_meridian.hpp	3174	0644	edit dl rm
result_direct.hpp	879	0644	edit dl rm
result_inverse.hpp	903	0644	edit dl rm
sjoberg_intersection.hpp	43099	0644	edit dl rm
spherical.hpp	8788	0644	edit dl rm
thomas_direct.hpp	9256	0644	edit dl rm
thomas_inverse.hpp	7768	0644	edit dl rm
unit_spheroid.hpp	1199	0644	edit dl rm
vertex_latitude.hpp	4489	0644	edit dl rm
vertex_longitude.hpp	10613	0644	edit dl rm
vincenty_direct.hpp	7105	0644	edit dl rm
vincenty_inverse.hpp	8254	0644	edit dl rm

Edit: /usr/include/boost/geometry/formulas/meridian_segment.hpp (2009B)

// Boost.Geometry

// Contributed and/or modified by Vissarion Fysikopoulos, 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_FORMULAS_MERIDIAN_SEGMENT_HPP
#define BOOST_GEOMETRY_FORMULAS_MERIDIAN_SEGMENT_HPP

#include <boost/math/constants/constants.hpp>

#include <boost/geometry/core/radius.hpp>

#include <boost/geometry/util/condition.hpp>
#include <boost/geometry/util/math.hpp>
#include <boost/geometry/util/normalize_spheroidal_coordinates.hpp>

namespace boost { namespace geometry { namespace formula
{

/*!
\brief Test if a segment is meridian or not.
*/

class meridian_segment
{

public :

enum SegmentType {NonMeridian, MeridianCrossingPole, MeridianNotCrossingPole};

template <typename T>
    static inline SegmentType is_meridian(T lon1, T lat1, T lon2, T lat2)
    {
        SegmentType res = NonMeridian;
        T diff = geometry::math::longitude_distance_signed<geometry::radian>(lon1, lon2);

if ( meridian_not_crossing_pole(lat1, lat2, diff) )
        {
            res = MeridianNotCrossingPole;
        }
        else if ( meridian_crossing_pole(diff) )
        {
            res = MeridianCrossingPole;
        }
        return res;
    }

template <typename T>
    static bool meridian_not_crossing_pole(T lat1, T lat2, T diff)
    {
        T half_pi = math::half_pi<T>();
        return math::equals(diff, T(0)) ||
               (math::equals(lat2, half_pi) && math::equals(lat1, -half_pi));
    }

template <typename T>
    static bool meridian_crossing_pole(T diff)
    {
        return math::equals(math::abs(diff), math::pi<T>());
    }

};

}}} // namespace boost::geometry::formula

#endif //BOOST_GEOMETRY_FORMULAS_MERIDIAN_SEGMENT_HPP

// Boost.Geometry // Copyright (c) 2017-2018 Oracle and/or its affiliates. // Contributed and/or modified by Vissarion Fysikopoulos, 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_FORMULAS_MERIDIAN_SEGMENT_HPP #define BOOST_GEOMETRY_FORMULAS_MERIDIAN_SEGMENT_HPP #include #include #include #include #include namespace boost { namespace geometry { namespace formula { /*! \brief Test if a segment is meridian or not. */ class meridian_segment { public : enum SegmentType {NonMeridian, MeridianCrossingPole, MeridianNotCrossingPole}; template static inline SegmentType is_meridian(T lon1, T lat1, T lon2, T lat2) { SegmentType res = NonMeridian; T diff = geometry::math::longitude_distance_signed(lon1, lon2); if ( meridian_not_crossing_pole(lat1, lat2, diff) ) { res = MeridianNotCrossingPole; } else if ( meridian_crossing_pole(diff) ) { res = MeridianCrossingPole; } return res; } template static bool meridian_not_crossing_pole(T lat1, T lat2, T diff) { T half_pi = math::half_pi(); return math::equals(diff, T(0)) || (math::equals(lat2, half_pi) && math::equals(lat1, -half_pi)); } template static bool meridian_crossing_pole(T diff) { return math::equals(math::abs(diff), math::pi()); } }; }}} // namespace boost::geometry::formula #endif //BOOST_GEOMETRY_FORMULAS_MERIDIAN_SEGMENT_HPP