Geometry azimuth does not compile for geographic coordinates

[Ken Barker <ken.barker@…>]

Calling azimuth to find the relative bearing of one point to another in geographic coordinates does not compile, e.g.:

    typedef boost::geometry::cs::geographic<boost::geometry::radian> Wgs84Coords;
    typedef boost::geometry::model::point<double, 2, Wgs84Coords> GeographicPoint;
    typedef boost::geometry::srs::spheroid<double> SpheroidType;

    double azimuth(GeographicPoint boost_a, GeographicPoint boost_b)
    {
      return boost::geometry::detail::azimuth<double>(boost_a, boost_b);
    }

or

    double wgs84_azimuth(GeographicPoint boost_a, GeographicPoint boost_b)
    {
      SpheroidType spheriod(wgs84::a, wgs84::b);
      return boost::geometry::detail::azimuth<double>(boost_a, boost_b, spheriod);
    }

The problem seems to be caused on line 52 of azimuth.hpp in geometry/algorithms/detail, it is:

    return geometry::detail::vincenty_inverse<ReturnType, false, true>

when it should be:

    return geometry::detail::vincenty_inverse<ReturnType, false, true>::apply

The amended line compiles and gives reasonable results.

[awulkiew]

Thanks!

Fix: ​https://github.com/boostorg/geometry/commit/2442cdd1e3e618c243ded917a8050df10e684d64