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/authalic_radius_sqr.hpp (2654B)

// Boost.Geometry

// 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_AUTHALIC_RADIUS_SQR_HPP
#define BOOST_GEOMETRY_FORMULAS_AUTHALIC_RADIUS_SQR_HPP

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

#include <boost/geometry/formulas/eccentricity_sqr.hpp>

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

#include <boost/geometry/algorithms/not_implemented.hpp>

#include <boost/math/special_functions/atanh.hpp>

namespace boost { namespace geometry
{

#ifndef DOXYGEN_NO_DISPATCH
namespace formula_dispatch
{

template <typename ResultType, typename Geometry, typename Tag = typename tag<Geometry>::type>
struct authalic_radius_sqr
    : not_implemented<Tag>
{};

template <typename ResultType, typename Geometry>
struct authalic_radius_sqr<ResultType, Geometry, srs_sphere_tag>
{
    static inline ResultType apply(Geometry const& geometry)
    {
        return math::sqr<ResultType>(get_radius<0>(geometry));
    }
};

template <typename ResultType, typename Geometry>
struct authalic_radius_sqr<ResultType, Geometry, srs_spheroid_tag>
{
    static inline ResultType apply(Geometry const& geometry)
    {
        ResultType const a2 = math::sqr<ResultType>(get_radius<0>(geometry));
        ResultType const e2 = formula::eccentricity_sqr<ResultType>(geometry);

return apply(a2, e2);
    }

static inline ResultType apply(ResultType const& a2, ResultType const& e2)
    {
        ResultType const c0 = 0;

if (math::equals(e2, c0))
        {
            return a2;
        }

ResultType const e = math::sqrt(e2);
        ResultType const c2 = 2;

//ResultType const b2 = math::sqr(get_radius<2>(geometry));
        //return a2 / c2 + b2 * boost::math::atanh(e) / (c2 * e);

ResultType const c1 = 1;
        return (a2 / c2) * ( c1 + (c1 - e2) * boost::math::atanh(e) / e );
    }
};

} // namespace formula_dispatch
#endif // DOXYGEN_NO_DISPATCH

#ifndef DOXYGEN_NO_DETAIL
namespace formula
{

template <typename ResultType, typename Geometry>
inline ResultType authalic_radius_sqr(Geometry const& geometry)
{
    return formula_dispatch::authalic_radius_sqr<ResultType, Geometry>::apply(geometry);
}

} // namespace formula
#endif // DOXYGEN_NO_DETAIL

}} // namespace boost::geometry

#endif // BOOST_GEOMETRY_FORMULAS_AUTHALIC_RADIUS_SQR_HPP

// Boost.Geometry // Copyright (c) 2017 Oracle and/or its affiliates. // 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_AUTHALIC_RADIUS_SQR_HPP #define BOOST_GEOMETRY_FORMULAS_AUTHALIC_RADIUS_SQR_HPP #include #include #include #include #include #include #include namespace boost { namespace geometry { #ifndef DOXYGEN_NO_DISPATCH namespace formula_dispatch { template ::type> struct authalic_radius_sqr : not_implemented {}; template struct authalic_radius_sqr { static inline ResultType apply(Geometry const& geometry) { return math::sqr(get_radius<0>(geometry)); } }; template struct authalic_radius_sqr { static inline ResultType apply(Geometry const& geometry) { ResultType const a2 = math::sqr(get_radius<0>(geometry)); ResultType const e2 = formula::eccentricity_sqr(geometry); return apply(a2, e2); } static inline ResultType apply(ResultType const& a2, ResultType const& e2) { ResultType const c0 = 0; if (math::equals(e2, c0)) { return a2; } ResultType const e = math::sqrt(e2); ResultType const c2 = 2; //ResultType const b2 = math::sqr(get_radius<2>(geometry)); //return a2 / c2 + b2 * boost::math::atanh(e) / (c2 * e); ResultType const c1 = 1; return (a2 / c2) * ( c1 + (c1 - e2) * boost::math::atanh(e) / e ); } }; } // namespace formula_dispatch #endif // DOXYGEN_NO_DISPATCH #ifndef DOXYGEN_NO_DETAIL namespace formula { template inline ResultType authalic_radius_sqr(Geometry const& geometry) { return formula_dispatch::authalic_radius_sqr::apply(geometry); } } // namespace formula #endif // DOXYGEN_NO_DETAIL }} // namespace boost::geometry #endif // BOOST_GEOMETRY_FORMULAS_AUTHALIC_RADIUS_SQR_HPP