1 files changed, 129 insertions, 0 deletions
diff --git a/boost/geometry/strategies/cartesian/line_interpolate.hpp b/boost/geometry/strategies/cartesian/line_interpolate.hpp
new file mode 100644
index 0000000000..ffad476ee4
--- /dev/null
+++ b/boost/geometry/strategies/cartesian/line_interpolate.hpp
@@ -0,0 +1,129 @@
+// Boost.Geometry
+
+// Copyright (c) 2018, Oracle and/or its affiliates.
+
+// Contributed and/or modified by Vissarion Fysikopoulos, on behalf of Oracle
+
+// Licensed under the Boost Software License version 1.0.
+// http://www.boost.org/users/license.html
+
+#ifndef BOOST_GEOMETRY_STRATEGIES_CARTESIAN_LINE_INTERPOLATE_HPP
+#define BOOST_GEOMETRY_STRATEGIES_CARTESIAN_LINE_INTERPOLATE_HPP
+
+#include <boost/geometry/core/assert.hpp>
+#include <boost/geometry/core/coordinate_dimension.hpp>
+#include <boost/geometry/core/coordinate_type.hpp>
+#include <boost/geometry/strategies/line_interpolate.hpp>
+#include <boost/geometry/strategies/cartesian/distance_pythagoras.hpp>
+#include <boost/geometry/util/select_calculation_type.hpp>
+
+
+namespace boost { namespace geometry
+{
+
+namespace strategy { namespace line_interpolate
+{
+
+
+/*!
+\brief Interpolate point on a cartesian segment.
+\ingroup strategies
+\tparam CalculationType \tparam_calculation
+\tparam DistanceStrategy The underlying point-point distance strategy
+
+\qbk{
+[heading See also]
+\* [link geometry.reference.algorithms.line_interpolate.line_interpolate_4_with_strategy line_interpolate (with strategy)]
+}
+
+*/
+template
+<
+    typename CalculationType = void,
+    typename DistanceStrategy = distance::pythagoras<CalculationType>
+>
+class cartesian
+{
+public:
+
+    // point-point strategy getters
+    struct distance_pp_strategy
+    {
+        typedef DistanceStrategy type;
+    };
+
+    inline typename distance_pp_strategy::type get_distance_pp_strategy() const
+    {
+        typedef typename distance_pp_strategy::type distance_type;
+        return distance_type();
+    }
+
+    template <typename Point, typename Fraction, typename Distance>
+    inline void apply(Point const& p0,
+                      Point const& p1,
+                      Fraction const& fraction,
+                      Point & p,
+                      Distance const&) const
+    {
+        typedef typename select_calculation_type_alt
+            <
+                CalculationType,
+                Point
+            >::type calc_t;
+
+        typedef model::point
+            <
+                calc_t,
+                geometry::dimension<Point>::value,
+                cs::cartesian
+            > calc_point_t;
+
+        calc_point_t cp0, cp1;
+        geometry::detail::conversion::convert_point_to_point(p0, cp0);
+        geometry::detail::conversion::convert_point_to_point(p1, cp1);
+
+        //segment convex combination: p0*fraction + p1*(1-fraction)
+        Fraction const one_minus_fraction = 1-fraction;
+        for_each_coordinate(cp1, detail::value_operation
+                                 <
+                                    Fraction,
+                                    std::multiplies
+                                 >(fraction));
+        for_each_coordinate(cp0, detail::value_operation
+                                 <
+                                    Fraction,
+                                    std::multiplies
+                                 >(one_minus_fraction));
+        for_each_coordinate(cp1, detail::point_operation
+                                 <
+                                    calc_point_t,
+                                    std::plus
+                                 >(cp0));
+
+        assert_dimension_equal<calc_point_t, Point>();
+        geometry::detail::conversion::convert_point_to_point(cp1, p);
+    }
+};
+
+
+#ifndef DOXYGEN_NO_STRATEGY_SPECIALIZATIONS
+namespace services
+{
+
+template <>
+struct default_strategy<cartesian_tag>
+{
+    typedef strategy::line_interpolate::cartesian<> type;
+};
+
+
+} // namespace services
+#endif // DOXYGEN_NO_STRATEGY_SPECIALIZATIONS
+
+
+}} // namespace strategy::line_interpolate
+
+
+}} // namespace boost::geometry
+
+#endif // BOOST_GEOMETRY_STRATEGIES_CARTESIAN_LINE_INTERPOLATE_HPP