Chris@102: // Boost.Geometry (aka GGL, Generic Geometry Library)
Chris@102: 
Chris@102: // Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands.
Chris@102: 
Chris@102: // This file was modified by Oracle on 2013, 2014.
Chris@102: // Modifications copyright (c) 2013-2014 Oracle and/or its affiliates.
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: // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
Chris@102: 
Chris@102: #ifndef BOOST_GEOMETRY_ALGORITHMS_DETAIL_RELATE_RELATE_HPP
Chris@102: #define BOOST_GEOMETRY_ALGORITHMS_DETAIL_RELATE_RELATE_HPP
Chris@102: 
Chris@102: #include <cstddef>
Chris@102: 
Chris@102: #include <boost/concept_check.hpp>
Chris@102: #include <boost/range.hpp>
Chris@102: 
Chris@102: #include <boost/mpl/if.hpp>
Chris@102: #include <boost/mpl/or.hpp>
Chris@102: #include <boost/type_traits/is_base_of.hpp>
Chris@102: 
Chris@102: #include <boost/geometry/algorithms/make.hpp>
Chris@102: #include <boost/geometry/algorithms/not_implemented.hpp>
Chris@102: 
Chris@102: #include <boost/geometry/core/access.hpp>
Chris@102: #include <boost/geometry/core/closure.hpp>
Chris@102: #include <boost/geometry/core/cs.hpp>
Chris@102: #include <boost/geometry/core/exterior_ring.hpp>
Chris@102: #include <boost/geometry/core/interior_rings.hpp>
Chris@102: #include <boost/geometry/core/point_order.hpp>
Chris@102: #include <boost/geometry/core/ring_type.hpp>
Chris@102: #include <boost/geometry/core/interior_rings.hpp>
Chris@102: #include <boost/geometry/core/tags.hpp>
Chris@102: 
Chris@102: #include <boost/geometry/geometries/concepts/check.hpp>
Chris@102: #include <boost/geometry/strategies/concepts/within_concept.hpp>
Chris@102: #include <boost/geometry/strategies/default_strategy.hpp>
Chris@102: #include <boost/geometry/strategies/within.hpp>
Chris@102: #include <boost/geometry/util/math.hpp>
Chris@102: #include <boost/geometry/util/order_as_direction.hpp>
Chris@102: #include <boost/geometry/views/closeable_view.hpp>
Chris@102: #include <boost/geometry/views/reversible_view.hpp>
Chris@102: 
Chris@102: #include <boost/geometry/algorithms/detail/relate/result.hpp>
Chris@102: 
Chris@102: #include <boost/geometry/algorithms/detail/relate/point_point.hpp>
Chris@102: #include <boost/geometry/algorithms/detail/relate/point_geometry.hpp>
Chris@102: #include <boost/geometry/algorithms/detail/relate/linear_linear.hpp>
Chris@102: #include <boost/geometry/algorithms/detail/relate/linear_areal.hpp>
Chris@102: #include <boost/geometry/algorithms/detail/relate/areal_areal.hpp>
Chris@102: 
Chris@102: namespace boost { namespace geometry
Chris@102: {
Chris@102: 
Chris@102: #ifndef DOXYGEN_NO_DETAIL
Chris@102: namespace detail { namespace relate {
Chris@102: 
Chris@102: // Those are used only to allow dispatch::relate to produce compile-time error
Chris@102: 
Chris@102: template <typename Geometry,
Chris@102:           typename Tag = typename geometry::tag<Geometry>::type>
Chris@102: struct is_supported_by_generic
Chris@102: {
Chris@102:     static const bool value
Chris@102:         = boost::is_same<Tag, linestring_tag>::value
Chris@102:        || boost::is_same<Tag, multi_linestring_tag>::value
Chris@102:        || boost::is_same<Tag, ring_tag>::value
Chris@102:        || boost::is_same<Tag, polygon_tag>::value
Chris@102:        || boost::is_same<Tag, multi_polygon_tag>::value;
Chris@102: };
Chris@102: 
Chris@102: template <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: struct is_generic
Chris@102: {
Chris@102:     static const bool value = is_supported_by_generic<Geometry1>::value
Chris@102:                            && is_supported_by_generic<Geometry2>::value;
Chris@102: };
Chris@102: 
Chris@102: 
Chris@102: template <typename Point, typename Geometry, typename Tag>
Chris@102: struct is_generic<Point, Geometry, point_tag, Tag>
Chris@102: {
Chris@102:     static const bool value = is_supported_by_generic<Geometry>::value;
Chris@102: };
Chris@102: 
Chris@102: template <typename Geometry, typename Point, typename Tag>
Chris@102: struct is_generic<Geometry, Point, Tag, point_tag>
Chris@102: {
Chris@102:     static const bool value = is_supported_by_generic<Geometry>::value;
Chris@102: };
Chris@102: 
Chris@102: template <typename Point1, typename Point2>
Chris@102: struct is_generic<Point1, Point2, point_tag, point_tag>
Chris@102: {
Chris@102:     static const bool value = false;
Chris@102: };
Chris@102: 
Chris@102: 
Chris@102: }} // namespace detail::relate
Chris@102: 
Chris@102: #ifndef DOXYGEN_NO_DISPATCH
Chris@102: namespace detail_dispatch { namespace relate {
Chris@102: 
Chris@102: 
Chris@102: template <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:           int TopDim1 = geometry::topological_dimension<Geometry1>::value,
Chris@102:           int TopDim2 = geometry::topological_dimension<Geometry2>::value,
Chris@102:           bool IsGeneric = detail::relate::is_generic<Geometry1, Geometry2>::value
Chris@102: >
Chris@102: struct relate : not_implemented<Tag1, Tag2>
Chris@102: {};
Chris@102: 
Chris@102: 
Chris@102: template <typename Point1, typename Point2>
Chris@102: struct relate<Point1, Point2, point_tag, point_tag, 0, 0, false>
Chris@102:     : detail::relate::point_point<Point1, Point2>
Chris@102: {};
Chris@102: 
Chris@102: template <typename Point, typename MultiPoint>
Chris@102: struct relate<Point, MultiPoint, point_tag, multi_point_tag, 0, 0, false>
Chris@102:     : detail::relate::point_multipoint<Point, MultiPoint>
Chris@102: {};
Chris@102: 
Chris@102: template <typename MultiPoint, typename Point>
Chris@102: struct relate<MultiPoint, Point, multi_point_tag, point_tag, 0, 0, false>
Chris@102:     : detail::relate::multipoint_point<MultiPoint, Point>
Chris@102: {};
Chris@102: 
Chris@102: template <typename MultiPoint1, typename MultiPoint2>
Chris@102: struct relate<MultiPoint1, MultiPoint2, multi_point_tag, multi_point_tag, 0, 0, false>
Chris@102:     : detail::relate::multipoint_multipoint<MultiPoint1, MultiPoint2>
Chris@102: {};
Chris@102: 
Chris@102: //template <typename Point, typename Box, int TopDim2>
Chris@102: //struct relate<Point, Box, point_tag, box_tag, 0, TopDim2, false>
Chris@102: //    : detail::relate::point_box<Point, Box>
Chris@102: //{};
Chris@102: //
Chris@102: //template <typename Box, typename Point, int TopDim1>
Chris@102: //struct relate<Box, Point, box_tag, point_tag, TopDim1, 0, false>
Chris@102: //    : detail::relate::box_point<Box, Point>
Chris@102: //{};
Chris@102: 
Chris@102: 
Chris@102: template <typename Point, typename Geometry, typename Tag2, int TopDim2>
Chris@102: struct relate<Point, Geometry, point_tag, Tag2, 0, TopDim2, true>
Chris@102:     : detail::relate::point_geometry<Point, Geometry>
Chris@102: {};
Chris@102: 
Chris@102: template <typename Geometry, typename Point, typename Tag1, int TopDim1>
Chris@102: struct relate<Geometry, Point, Tag1, point_tag, TopDim1, 0, true>
Chris@102:     : detail::relate::geometry_point<Geometry, Point>
Chris@102: {};
Chris@102: 
Chris@102: 
Chris@102: template <typename Linear1, typename Linear2, typename Tag1, typename Tag2>
Chris@102: struct relate<Linear1, Linear2, Tag1, Tag2, 1, 1, true>
Chris@102:     : detail::relate::linear_linear<Linear1, Linear2>
Chris@102: {};
Chris@102: 
Chris@102: 
Chris@102: template <typename Linear, typename Areal, typename Tag1, typename Tag2>
Chris@102: struct relate<Linear, Areal, Tag1, Tag2, 1, 2, true>
Chris@102:     : detail::relate::linear_areal<Linear, Areal>
Chris@102: {};
Chris@102: 
Chris@102: template <typename Areal, typename Linear, typename Tag1, typename Tag2>
Chris@102: struct relate<Areal, Linear, Tag1, Tag2, 2, 1, true>
Chris@102:     : detail::relate::areal_linear<Areal, Linear>
Chris@102: {};
Chris@102: 
Chris@102: 
Chris@102: template <typename Areal1, typename Areal2, typename Tag1, typename Tag2>
Chris@102: struct relate<Areal1, Areal2, Tag1, Tag2, 2, 2, true>
Chris@102:     : detail::relate::areal_areal<Areal1, Areal2>
Chris@102: {};
Chris@102: 
Chris@102: 
Chris@102: }} // namespace detail_dispatch::relate
Chris@102: #endif // DOXYGEN_NO_DISPATCH
Chris@102: 
Chris@102: namespace detail { namespace relate {
Chris@102: 
Chris@102: template <typename Geometry1, typename Geometry2>
Chris@102: struct interruption_enabled
Chris@102: {
Chris@102:     static const bool value =
Chris@102:         detail_dispatch::relate::relate<Geometry1, Geometry2>::interruption_enabled;
Chris@102: };
Chris@102: 
Chris@102: template <typename Geometry1,
Chris@102:           typename Geometry2,
Chris@102:           typename Result,
Chris@102:           bool IsSequence = boost::mpl::is_sequence<Result>::value>
Chris@102: struct result_handler_type
Chris@102:     : not_implemented<Result>
Chris@102: {};
Chris@102: 
Chris@102: template <typename Geometry1, typename Geometry2>
Chris@102: struct result_handler_type<Geometry1, Geometry2, matrix9, false>
Chris@102: {
Chris@102:     typedef matrix_handler<matrix9> type;
Chris@102: };
Chris@102: 
Chris@102: template <typename Geometry1, typename Geometry2>
Chris@102: struct result_handler_type<Geometry1, Geometry2, mask9, false>
Chris@102: {
Chris@102:     typedef mask_handler
Chris@102:         <
Chris@102:             mask9,
Chris@102:             interruption_enabled
Chris@102:                 <
Chris@102:                     Geometry1,
Chris@102:                     Geometry2
Chris@102:                 >::value
Chris@102:         > type;
Chris@102: };
Chris@102: 
Chris@102: template <typename Geometry1, typename Geometry2, typename Head, typename Tail>
Chris@102: struct result_handler_type<Geometry1, Geometry2, boost::tuples::cons<Head, Tail>, false>
Chris@102: {
Chris@102:     typedef mask_handler
Chris@102:         <
Chris@102:             boost::tuples::cons<Head, Tail>,
Chris@102:             interruption_enabled
Chris@102:                 <
Chris@102:                     Geometry1,
Chris@102:                     Geometry2
Chris@102:                 >::value
Chris@102:         > type;
Chris@102: };
Chris@102: 
Chris@102: template <typename Geometry1, typename Geometry2,
Chris@102:           char II, char IB, char IE,
Chris@102:           char BI, char BB, char BE,
Chris@102:           char EI, char EB, char EE>
Chris@102: struct result_handler_type<Geometry1, Geometry2, static_mask<II, IB, IE, BI, BB, BE, EI, EB, EE>, false>
Chris@102: {
Chris@102:     typedef static_mask_handler
Chris@102:         <
Chris@102:             static_mask<II, IB, IE, BI, BB, BE, EI, EB, EE>,
Chris@102:             interruption_enabled
Chris@102:                 <
Chris@102:                     Geometry1,
Chris@102:                     Geometry2
Chris@102:                 >::value
Chris@102:         > type;
Chris@102: };
Chris@102: 
Chris@102: template <typename Geometry1, typename Geometry2, typename StaticSequence>
Chris@102: struct result_handler_type<Geometry1, Geometry2, StaticSequence, true>
Chris@102: {
Chris@102:     typedef static_mask_handler
Chris@102:         <
Chris@102:             StaticSequence,
Chris@102:             interruption_enabled
Chris@102:                 <
Chris@102:                     Geometry1,
Chris@102:                     Geometry2
Chris@102:                 >::value
Chris@102:         > type;
Chris@102: };
Chris@102: 
Chris@102: template <typename MatrixOrMask, typename Geometry1, typename Geometry2>
Chris@102: inline
Chris@102: typename result_handler_type
Chris@102:     <
Chris@102:         Geometry1,
Chris@102:         Geometry2,
Chris@102:         MatrixOrMask
Chris@102:     >::type::result_type
Chris@102: relate(Geometry1 const& geometry1,
Chris@102:        Geometry2 const& geometry2,
Chris@102:        MatrixOrMask const& matrix_or_mask = MatrixOrMask())
Chris@102: {
Chris@102:     typedef typename result_handler_type
Chris@102:         <
Chris@102:             Geometry1,
Chris@102:             Geometry2,
Chris@102:             MatrixOrMask
Chris@102:         >::type handler_type;
Chris@102: 
Chris@102:     handler_type handler(matrix_or_mask);
Chris@102:     detail_dispatch::relate::relate<Geometry1, Geometry2>::apply(geometry1, geometry2, handler);
Chris@102:     return handler.result();
Chris@102: }
Chris@102: 
Chris@102: struct implemented_tag {};
Chris@102: 
Chris@102: template <template <typename, typename> class StaticMaskTrait,
Chris@102:           typename Geometry1,
Chris@102:           typename Geometry2>
Chris@102: struct relate_base
Chris@102:     : boost::mpl::if_
Chris@102:         <
Chris@102:             boost::mpl::or_
Chris@102:                 <
Chris@102:                     boost::is_base_of
Chris@102:                         <
Chris@102:                             nyi::not_implemented_tag,
Chris@102:                             StaticMaskTrait<Geometry1, Geometry2>
Chris@102:                         >,
Chris@102:                     boost::is_base_of
Chris@102:                         <
Chris@102:                             nyi::not_implemented_tag,
Chris@102:                             detail_dispatch::relate::relate<Geometry1, Geometry2>
Chris@102:                         >
Chris@102:                 >,
Chris@102:             not_implemented
Chris@102:                 <
Chris@102:                     typename geometry::tag<Geometry1>::type,
Chris@102:                     typename geometry::tag<Geometry2>::type
Chris@102:                 >,
Chris@102:             implemented_tag
Chris@102:         >::type
Chris@102: {
Chris@102:     static inline bool apply(Geometry1 const& g1, Geometry2 const& g2)
Chris@102:     {
Chris@102:         typedef typename StaticMaskTrait<Geometry1, Geometry2>::type static_mask;
Chris@102:         return detail::relate::relate<static_mask>(g1, g2);
Chris@102:     }
Chris@102: };
Chris@102: 
Chris@102: }} // namespace detail::relate
Chris@102: #endif // DOXYGEN_NO_DETAIL
Chris@102: 
Chris@102: }} // namespace boost::geometry
Chris@102: 
Chris@102: #endif // BOOST_GEOMETRY_ALGORITHMS_DETAIL_RELATE_RELATE_HPP