Chris@102: // Boost.Geometry (aka GGL, Generic Geometry Library)
Chris@102: 
Chris@102: // Copyright (c) 2013 Barend Gehrels, Amsterdam, the Netherlands.
Chris@102: // Copyright (c) 2013 Bruno Lalande, Paris, France.
Chris@102: // Copyright (c) 2013 Mateusz Loskot, London, UK.
Chris@102: // Copyright (c) 2013 Adam Wulkiewicz, Lodz, Poland.
Chris@102: 
Chris@102: // Use, modification and distribution is subject to the Boost Software License,
Chris@102: // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
Chris@102: // http://www.boost.org/LICENSE_1_0.txt)
Chris@102: 
Chris@102: #ifndef BOOST_GEOMETRY_ALGORITHMS_DETAIL_RECALCULATE_HPP
Chris@102: #define BOOST_GEOMETRY_ALGORITHMS_DETAIL_RECALCULATE_HPP
Chris@102: 
Chris@102: 
Chris@102: #include <cstddef>
Chris@102: 
Chris@102: #include <boost/concept/requires.hpp>
Chris@102: #include <boost/concept_check.hpp>
Chris@102: #include <boost/mpl/assert.hpp>
Chris@102: #include <boost/mpl/if.hpp>
Chris@102: #include <boost/numeric/conversion/bounds.hpp>
Chris@102: #include <boost/numeric/conversion/cast.hpp>
Chris@102: #include <boost/type_traits.hpp>
Chris@102: 
Chris@102: #include <boost/geometry/arithmetic/arithmetic.hpp>
Chris@102: #include <boost/geometry/algorithms/append.hpp>
Chris@102: #include <boost/geometry/algorithms/clear.hpp>
Chris@102: #include <boost/geometry/core/access.hpp>
Chris@102: #include <boost/geometry/core/interior_rings.hpp>
Chris@102: #include <boost/geometry/core/exterior_ring.hpp>
Chris@102: #include <boost/geometry/core/tags.hpp>
Chris@102: 
Chris@102: #include <boost/geometry/geometries/concepts/check.hpp>
Chris@102: 
Chris@102: 
Chris@102: namespace boost { namespace geometry
Chris@102: {
Chris@102: 
Chris@102: #ifndef DOXYGEN_NO_DETAIL
Chris@102: namespace detail { namespace recalculate
Chris@102: {
Chris@102: 
Chris@102: template <std::size_t Dimension>
Chris@102: struct recalculate_point
Chris@102: {
Chris@102:     template <typename Point1, typename Point2, typename Strategy>
Chris@102:     static inline void apply(Point1& point1, Point2 const& point2, Strategy const& strategy)
Chris@102:     {
Chris@102:         std::size_t const dim = Dimension - 1;
Chris@102:         geometry::set<dim>(point1, strategy.template apply<dim>(geometry::get<dim>(point2)));
Chris@102:         recalculate_point<dim>::apply(point1, point2, strategy);
Chris@102:     }
Chris@102: };
Chris@102: 
Chris@102: template <>
Chris@102: struct recalculate_point<0>
Chris@102: {
Chris@102:     template <typename Point1, typename Point2, typename Strategy>
Chris@102:     static inline void apply(Point1&, Point2 const&, Strategy const&)
Chris@102:     {
Chris@102:     }
Chris@102: };
Chris@102: 
Chris@102: 
Chris@102: template <std::size_t Dimension>
Chris@102: struct recalculate_indexed
Chris@102: {
Chris@102:     template <typename Geometry1, typename Geometry2, typename Strategy>
Chris@102:     static inline void apply(Geometry1& geometry1, Geometry2 const& geometry2, Strategy const& strategy)
Chris@102:     {
Chris@102:         // Do it for both indices in one dimension
Chris@102:         static std::size_t const dim = Dimension - 1;
Chris@102:         geometry::set<0, dim>(geometry1, strategy.template apply<dim>(geometry::get<0, dim>(geometry2)));
Chris@102:         geometry::set<1, dim>(geometry1, strategy.template apply<dim>(geometry::get<1, dim>(geometry2)));
Chris@102:         recalculate_indexed<dim>::apply(geometry1, geometry2, strategy);
Chris@102:     }
Chris@102: };
Chris@102: 
Chris@102: template <>
Chris@102: struct recalculate_indexed<0>
Chris@102: {
Chris@102: 
Chris@102:     template <typename Geometry1, typename Geometry2, typename Strategy>
Chris@102:     static inline void apply(Geometry1& , Geometry2 const& , Strategy const& )
Chris@102:     {
Chris@102:     }
Chris@102: };
Chris@102: 
Chris@102: struct range_to_range
Chris@102: {
Chris@102:     template
Chris@102:     <
Chris@102:         typename Range1,
Chris@102:         typename Range2,
Chris@102:         typename Strategy
Chris@102:     >
Chris@102:     static inline void apply(Range1& destination, Range2 const& source,
Chris@102:             Strategy const& strategy)
Chris@102:     {
Chris@102:         typedef typename geometry::point_type<Range2>::type point_type;
Chris@102:         typedef recalculate_point<geometry::dimension<point_type>::value> per_point;
Chris@102:         geometry::clear(destination);
Chris@102: 
Chris@102:         for (typename boost::range_iterator<Range2 const>::type it
Chris@102:                 = boost::begin(source);
Chris@102:             it != boost::end(source);
Chris@102:             ++it)
Chris@102:         {
Chris@102:             point_type p;
Chris@102:             per_point::apply(p, *it, strategy);
Chris@102:             geometry::append(destination, p);
Chris@102:         }
Chris@102:     }
Chris@102: };
Chris@102: 
Chris@102: struct polygon_to_polygon
Chris@102: {
Chris@102: private:
Chris@102:     template
Chris@102:     <
Chris@102:         typename IteratorIn,
Chris@102:         typename IteratorOut,
Chris@102:         typename Strategy
Chris@102:     >
Chris@102:     static inline void iterate(IteratorIn begin, IteratorIn end,
Chris@102:                     IteratorOut it_out,
Chris@102:                     Strategy const& strategy)
Chris@102:     {
Chris@102:         for (IteratorIn it_in = begin; it_in != end;  ++it_in, ++it_out)
Chris@102:         {
Chris@102:             range_to_range::apply(*it_out, *it_in, strategy);
Chris@102:         }
Chris@102:     }
Chris@102: 
Chris@102:     template
Chris@102:     <
Chris@102:         typename InteriorRingsOut,
Chris@102:         typename InteriorRingsIn,
Chris@102:         typename Strategy
Chris@102:     >
Chris@102:     static inline void apply_interior_rings(
Chris@102:                     InteriorRingsOut& interior_rings_out,
Chris@102:                     InteriorRingsIn const& interior_rings_in,
Chris@102:                     Strategy const& strategy)
Chris@102:     {
Chris@102:         traits::resize<InteriorRingsOut>::apply(interior_rings_out,
Chris@102:             boost::size(interior_rings_in));
Chris@102: 
Chris@102:         iterate(
Chris@102:             boost::begin(interior_rings_in), boost::end(interior_rings_in),
Chris@102:             boost::begin(interior_rings_out),
Chris@102:             strategy);
Chris@102:     }
Chris@102: 
Chris@102: public:
Chris@102:     template
Chris@102:     <
Chris@102:         typename Polygon1,
Chris@102:         typename Polygon2,
Chris@102:         typename Strategy
Chris@102:     >
Chris@102:     static inline void apply(Polygon1& destination, Polygon2 const& source,
Chris@102:             Strategy const& strategy)
Chris@102:     {
Chris@102:         range_to_range::apply(geometry::exterior_ring(destination),
Chris@102:             geometry::exterior_ring(source), strategy);
Chris@102: 
Chris@102:         apply_interior_rings(geometry::interior_rings(destination),
Chris@102:             geometry::interior_rings(source), strategy);
Chris@102:     }
Chris@102: };
Chris@102: 
Chris@102: }} // namespace detail::recalculate
Chris@102: #endif // DOXYGEN_NO_DETAIL
Chris@102: 
Chris@102: #ifndef DOXYGEN_NO_DISPATCH
Chris@102: namespace dispatch
Chris@102: {
Chris@102: 
Chris@102: template
Chris@102: <
Chris@102:     typename Geometry1,
Chris@102:     typename Geometry2,
Chris@102:     typename Tag1 = typename geometry::tag<Geometry1>::type,
Chris@102:     typename Tag2 = typename geometry::tag<Geometry2>::type
Chris@102: >
Chris@102: struct recalculate : not_implemented<Tag1, Tag2>
Chris@102: {};
Chris@102: 
Chris@102: template <typename Point1, typename Point2>
Chris@102: struct recalculate<Point1, Point2, point_tag, point_tag>
Chris@102:     : detail::recalculate::recalculate_point<geometry::dimension<Point1>::value>
Chris@102: {};
Chris@102: 
Chris@102: template <typename Box1, typename Box2>
Chris@102: struct recalculate<Box1, Box2, box_tag, box_tag>
Chris@102:     : detail::recalculate::recalculate_indexed<geometry::dimension<Box1>::value>
Chris@102: {};
Chris@102: 
Chris@102: template <typename Segment1, typename Segment2>
Chris@102: struct recalculate<Segment1, Segment2, segment_tag, segment_tag>
Chris@102:     : detail::recalculate::recalculate_indexed<geometry::dimension<Segment1>::value>
Chris@102: {};
Chris@102: 
Chris@102: template <typename Polygon1, typename Polygon2>
Chris@102: struct recalculate<Polygon1, Polygon2, polygon_tag, polygon_tag>
Chris@102:     : detail::recalculate::polygon_to_polygon
Chris@102: {};
Chris@102: 
Chris@102: } // namespace dispatch
Chris@102: #endif // DOXYGEN_NO_DISPATCH
Chris@102: 
Chris@102: 
Chris@102: 
Chris@102: template <typename Geometry1, typename Geometry2, typename Strategy>
Chris@102: inline void recalculate(Geometry1& geometry1, Geometry2 const& geometry2, Strategy const& strategy)
Chris@102: {
Chris@102:     concept::check<Geometry1>();
Chris@102:     concept::check<Geometry2 const>();
Chris@102: 
Chris@102:     // static assert dimensions (/types) are the same
Chris@102: 
Chris@102:     dispatch::recalculate<Geometry1, Geometry2>::apply(geometry1, geometry2, strategy);
Chris@102: }
Chris@102: 
Chris@102: 
Chris@102: }} // namespace boost::geometry
Chris@102: 
Chris@102: 
Chris@102: #endif // BOOST_GEOMETRY_ALGORITHMS_DETAIL_RECALCULATE_HPP