Chris@102: // Boost.Geometry (aka GGL, Generic Geometry Library)
Chris@102: 
Chris@102: // Copyright (c) 2014-2015 Barend Gehrels, Amsterdam, the Netherlands.
Chris@102: // Copyright (c) 2014-2015 Bruno Lalande, Paris, France.
Chris@102: // Copyright (c) 2014-2015 Mateusz Loskot, London, UK.
Chris@102: // Copyright (c) 2014-2015 Adam Wulkiewicz, Lodz, Poland.
Chris@102: 
Chris@102: // This file was modified by Oracle on 2015.
Chris@102: // Modifications copyright (c) 2015, Oracle and/or its affiliates.
Chris@102: 
Chris@102: // Contributed and/or modified by Menelaos Karavelas, on behalf of Oracle
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_POLICIES_ROBUSTNESS_GET_RESCALE_POLICY_HPP
Chris@102: #define BOOST_GEOMETRY_POLICIES_ROBUSTNESS_GET_RESCALE_POLICY_HPP
Chris@102: 
Chris@102: 
Chris@102: #include <cstddef>
Chris@102: 
Chris@102: #include <boost/type_traits.hpp>
Chris@102: #include <boost/mpl/assert.hpp>
Chris@102: 
Chris@102: #include <boost/geometry/core/tag_cast.hpp>
Chris@102: 
Chris@102: #include <boost/geometry/algorithms/envelope.hpp>
Chris@102: #include <boost/geometry/algorithms/expand.hpp>
Chris@102: #include <boost/geometry/algorithms/detail/recalculate.hpp>
Chris@102: #include <boost/geometry/algorithms/detail/get_max_size.hpp>
Chris@102: #include <boost/geometry/policies/robustness/robust_type.hpp>
Chris@102: 
Chris@102: #include <boost/geometry/geometries/point.hpp>
Chris@102: #include <boost/geometry/geometries/box.hpp>
Chris@102: 
Chris@102: #include <boost/geometry/policies/robustness/no_rescale_policy.hpp>
Chris@102: #include <boost/geometry/policies/robustness/rescale_policy.hpp>
Chris@102: 
Chris@102: #include <boost/geometry/util/promote_floating_point.hpp>
Chris@102: 
Chris@102: namespace boost { namespace geometry
Chris@102: {
Chris@102: 
Chris@102: #ifndef DOXYGEN_NO_DETAIL
Chris@102: namespace detail { namespace get_rescale_policy
Chris@102: {
Chris@102: 
Chris@102: template
Chris@102: <
Chris@102:     typename Box,
Chris@102:     typename Point,
Chris@102:     typename RobustPoint,
Chris@102:     typename Factor
Chris@102: >
Chris@102: inline void scale_box_to_integer_range(Box const& box,
Chris@102:                                        Point& min_point,
Chris@102:                                        RobustPoint& min_robust_point,
Chris@102:                                        Factor& factor)
Chris@102: {
Chris@102:     // Scale box to integer-range
Chris@102:     typedef typename promote_floating_point
Chris@102:         <
Chris@102:             typename geometry::coordinate_type<Point>::type
Chris@102:         >::type num_type;
Chris@102:     num_type const diff = boost::numeric_cast<num_type>(detail::get_max_size(box));
Chris@102:     num_type const range = 10000000.0; // Define a large range to get precise integer coordinates
Chris@102:     num_type const half = 0.5;
Chris@102:     factor = math::equals(diff, num_type()) || diff >= range ? 1
Chris@102:         : boost::numeric_cast<num_type>(
Chris@102:             boost::numeric_cast<boost::long_long_type>(half + range / diff));
Chris@102: 
Chris@102:     BOOST_ASSERT(factor >= 1);
Chris@102: 
Chris@102:     // Assign input/output minimal points
Chris@102:     detail::assign_point_from_index<0>(box, min_point);
Chris@102:     num_type const two = 2;
Chris@102:     boost::long_long_type const min_coordinate
Chris@102:         = boost::numeric_cast<boost::long_long_type>(-range / two);
Chris@102:     assign_values(min_robust_point, min_coordinate, min_coordinate);
Chris@102: }
Chris@102: 
Chris@102: template <typename Point, typename RobustPoint, typename Geometry, typename Factor>
Chris@102: static inline void init_rescale_policy(Geometry const& geometry,
Chris@102:         Point& min_point,
Chris@102:         RobustPoint& min_robust_point,
Chris@102:         Factor& factor)
Chris@102: {
Chris@102:     // Get bounding boxes
Chris@102:     model::box<Point> env = geometry::return_envelope<model::box<Point> >(geometry);
Chris@102: 
Chris@102:     scale_box_to_integer_range(env, min_point, min_robust_point, factor);
Chris@102: }
Chris@102: 
Chris@102: template <typename Point, typename RobustPoint, typename Geometry1, typename Geometry2, typename Factor>
Chris@102: static inline void init_rescale_policy(Geometry1 const& geometry1,
Chris@102:         Geometry2 const& geometry2,
Chris@102:         Point& min_point,
Chris@102:         RobustPoint& min_robust_point,
Chris@102:         Factor& factor)
Chris@102: {
Chris@102:     // Get bounding boxes
Chris@102:     model::box<Point> env = geometry::return_envelope<model::box<Point> >(geometry1);
Chris@102:     model::box<Point> env2 = geometry::return_envelope<model::box<Point> >(geometry2);
Chris@102:     geometry::expand(env, env2);
Chris@102: 
Chris@102:     scale_box_to_integer_range(env, min_point, min_robust_point, factor);
Chris@102: }
Chris@102: 
Chris@102: 
Chris@102: template
Chris@102: <
Chris@102:     typename Point,
Chris@102:     bool IsFloatingPoint
Chris@102: >
Chris@102: struct rescale_policy_type
Chris@102: {
Chris@102:     typedef no_rescale_policy type;
Chris@102: };
Chris@102: 
Chris@102: // We rescale only all FP types
Chris@102: template
Chris@102: <
Chris@102:     typename Point
Chris@102: >
Chris@102: struct rescale_policy_type<Point, true>
Chris@102: {
Chris@102:     typedef typename geometry::coordinate_type<Point>::type coordinate_type;
Chris@102:     typedef model::point
Chris@102:     <
Chris@102:         typename detail::robust_type<coordinate_type>::type,
Chris@102:         geometry::dimension<Point>::value,
Chris@102:         typename geometry::coordinate_system<Point>::type
Chris@102:     > robust_point_type;
Chris@102:     typedef typename promote_floating_point<coordinate_type>::type factor_type;
Chris@102:     typedef detail::robust_policy<Point, robust_point_type, factor_type> type;
Chris@102: };
Chris@102: 
Chris@102: template <typename Policy>
Chris@102: struct get_rescale_policy
Chris@102: {
Chris@102:     template <typename Geometry>
Chris@102:     static inline Policy apply(Geometry const& geometry)
Chris@102:     {
Chris@102:         typedef typename point_type<Geometry>::type point_type;
Chris@102:         typedef typename geometry::coordinate_type<Geometry>::type coordinate_type;
Chris@102:         typedef typename promote_floating_point<coordinate_type>::type factor_type;
Chris@102:         typedef model::point
Chris@102:         <
Chris@102:             typename detail::robust_type<coordinate_type>::type,
Chris@102:             geometry::dimension<point_type>::value,
Chris@102:             typename geometry::coordinate_system<point_type>::type
Chris@102:         > robust_point_type;
Chris@102: 
Chris@102:         point_type min_point;
Chris@102:         robust_point_type min_robust_point;
Chris@102:         factor_type factor;
Chris@102:         init_rescale_policy(geometry, min_point, min_robust_point, factor);
Chris@102: 
Chris@102:         return Policy(min_point, min_robust_point, factor);
Chris@102:     }
Chris@102: 
Chris@102:     template <typename Geometry1, typename Geometry2>
Chris@102:     static inline Policy apply(Geometry1 const& geometry1, Geometry2 const& geometry2)
Chris@102:     {
Chris@102:         typedef typename point_type<Geometry1>::type point_type;
Chris@102:         typedef typename geometry::coordinate_type<Geometry1>::type coordinate_type;
Chris@102:         typedef typename promote_floating_point<coordinate_type>::type factor_type;
Chris@102:         typedef model::point
Chris@102:         <
Chris@102:             typename detail::robust_type<coordinate_type>::type,
Chris@102:             geometry::dimension<point_type>::value,
Chris@102:             typename geometry::coordinate_system<point_type>::type
Chris@102:         > robust_point_type;
Chris@102: 
Chris@102:         point_type min_point;
Chris@102:         robust_point_type min_robust_point;
Chris@102:         factor_type factor;
Chris@102:         init_rescale_policy(geometry1, geometry2, min_point, min_robust_point, factor);
Chris@102: 
Chris@102:         return Policy(min_point, min_robust_point, factor);
Chris@102:     }
Chris@102: };
Chris@102: 
Chris@102: // Specialization for no-rescaling
Chris@102: template <>
Chris@102: struct get_rescale_policy<no_rescale_policy>
Chris@102: {
Chris@102:     template <typename Geometry>
Chris@102:     static inline no_rescale_policy apply(Geometry const& )
Chris@102:     {
Chris@102:         return no_rescale_policy();
Chris@102:     }
Chris@102: 
Chris@102:     template <typename Geometry1, typename Geometry2>
Chris@102:     static inline no_rescale_policy apply(Geometry1 const& , Geometry2 const& )
Chris@102:     {
Chris@102:         return no_rescale_policy();
Chris@102:     }
Chris@102: };
Chris@102: 
Chris@102: 
Chris@102: }} // namespace detail::get_rescale_policy
Chris@102: #endif // DOXYGEN_NO_DETAIL
Chris@102: 
Chris@102: template<typename Point>
Chris@102: struct rescale_policy_type
Chris@102:     : public detail::get_rescale_policy::rescale_policy_type
Chris@102:     <
Chris@102:         Point,
Chris@102: #if defined(BOOST_GEOMETRY_NO_ROBUSTNESS)
Chris@102:         false
Chris@102: #else
Chris@102:         boost::is_floating_point
Chris@102:         <
Chris@102:             typename geometry::coordinate_type<Point>::type
Chris@102:         >::type::value
Chris@102: #endif
Chris@102:     >
Chris@102: {
Chris@102:     static const bool is_point
Chris@102:         = boost::is_same
Chris@102:             <
Chris@102:                 typename geometry::tag<Point>::type,
Chris@102:                 geometry::point_tag
Chris@102:             >::type::value;
Chris@102: 
Chris@102:     BOOST_MPL_ASSERT_MSG((is_point),
Chris@102:                          INVALID_INPUT_GEOMETRY,
Chris@102:                          (typename geometry::tag<Point>::type));
Chris@102: };
Chris@102: 
Chris@102: 
Chris@102: template
Chris@102: <
Chris@102:     typename Geometry1,
Chris@102:     typename Geometry2,
Chris@102:     typename Tag1 = typename tag_cast
Chris@102:     <
Chris@102:         typename tag<Geometry1>::type,
Chris@102:         box_tag,
Chris@102:         pointlike_tag,
Chris@102:         linear_tag,
Chris@102:         areal_tag
Chris@102:     >::type,
Chris@102:     typename Tag2 = typename tag_cast
Chris@102:     <
Chris@102:         typename tag<Geometry2>::type,
Chris@102:         box_tag,
Chris@102:         pointlike_tag,
Chris@102:         linear_tag,
Chris@102:         areal_tag
Chris@102:     >::type
Chris@102: >
Chris@102: struct rescale_overlay_policy_type
Chris@102:     // Default: no rescaling
Chris@102:     : public detail::get_rescale_policy::rescale_policy_type
Chris@102:         <
Chris@102:             typename geometry::point_type<Geometry1>::type,
Chris@102:             false
Chris@102:         >
Chris@102: {};
Chris@102: 
Chris@102: // Areal/areal: get rescale policy based on coordinate type
Chris@102: template
Chris@102: <
Chris@102:     typename Geometry1,
Chris@102:     typename Geometry2
Chris@102: >
Chris@102: struct rescale_overlay_policy_type<Geometry1, Geometry2, areal_tag, areal_tag>
Chris@102:     : public rescale_policy_type
Chris@102:         <
Chris@102:             typename geometry::point_type<Geometry1>::type
Chris@102:         >
Chris@102: {};
Chris@102: 
Chris@102: 
Chris@102: template <typename Policy, typename Geometry>
Chris@102: inline Policy get_rescale_policy(Geometry const& geometry)
Chris@102: {
Chris@102:     return detail::get_rescale_policy::get_rescale_policy<Policy>::apply(geometry);
Chris@102: }
Chris@102: 
Chris@102: template <typename Policy, typename Geometry1, typename Geometry2>
Chris@102: inline Policy get_rescale_policy(Geometry1 const& geometry1, Geometry2 const& geometry2)
Chris@102: {
Chris@102:     return detail::get_rescale_policy::get_rescale_policy<Policy>::apply(geometry1, geometry2);
Chris@102: }
Chris@102: 
Chris@102: 
Chris@102: }} // namespace boost::geometry
Chris@102: 
Chris@102: 
Chris@102: #endif // BOOST_GEOMETRY_POLICIES_ROBUSTNESS_GET_RESCALE_POLICY_HPP