123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854 |
- // Boost.Geometry (aka GGL, Generic Geometry Library)
- // Copyright (c) 2011-2015 Barend Gehrels, Amsterdam, the Netherlands.
- // Copyright (c) 2017 Adam Wulkiewicz, Lodz, Poland.
- // This file was modified by Oracle on 2015-2020.
- // Modifications copyright (c) 2015-2020 Oracle and/or its affiliates.
- // Contributed and/or modified by Menelaos Karavelas, on behalf of Oracle
- // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
- // Use, modification and distribution is subject to the Boost Software License,
- // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
- // http://www.boost.org/LICENSE_1_0.txt)
- #ifndef BOOST_GEOMETRY_ALGORITHMS_DETAIL_PARTITION_HPP
- #define BOOST_GEOMETRY_ALGORITHMS_DETAIL_PARTITION_HPP
- #include <cstddef>
- #include <type_traits>
- #include <vector>
- #include <boost/range/begin.hpp>
- #include <boost/range/empty.hpp>
- #include <boost/range/end.hpp>
- #include <boost/range/size.hpp>
- #include <boost/geometry/algorithms/assign.hpp>
- #include <boost/geometry/core/access.hpp>
- #include <boost/geometry/core/coordinate_type.hpp>
- namespace boost { namespace geometry
- {
- namespace detail { namespace partition
- {
- template <typename T, bool IsIntegral = std::is_integral<T>::value>
- struct divide_interval
- {
- static inline T apply(T const& mi, T const& ma)
- {
- static T const two = 2;
- return (mi + ma) / two;
- }
- };
- template <typename T>
- struct divide_interval<T, true>
- {
- static inline T apply(T const& mi, T const& ma)
- {
- // Avoid overflow
- return mi / 2 + ma / 2 + (mi % 2 + ma % 2) / 2;
- }
- };
- struct visit_no_policy
- {
- template <typename Box>
- static inline void apply(Box const&, std::size_t )
- {}
- };
- struct include_all_policy
- {
- template <typename Item>
- static inline bool apply(Item const&)
- {
- return true;
- }
- };
- template <std::size_t Dimension, typename Box>
- inline void divide_box(Box const& box, Box& lower_box, Box& upper_box)
- {
- using coor_t = typename coordinate_type<Box>::type;
- // Divide input box into two halves
- // either left/right (Dimension 0)
- // or top/bottom (Dimension 1)
- coor_t const mid
- = divide_interval<coor_t>::apply(geometry::get<min_corner, Dimension>(box),
- geometry::get<max_corner, Dimension>(box));
- lower_box = box;
- upper_box = box;
- geometry::set<max_corner, Dimension>(lower_box, mid);
- geometry::set<min_corner, Dimension>(upper_box, mid);
- }
- // Divide forward_range into three subsets: lower, upper and oversized
- // (not-fitting)
- // (lower == left or bottom, upper == right or top)
- template <typename Box, typename IteratorVector, typename OverlapsPolicy>
- inline void divide_into_subsets(Box const& lower_box,
- Box const& upper_box,
- IteratorVector const& input,
- IteratorVector& lower,
- IteratorVector& upper,
- IteratorVector& exceeding,
- OverlapsPolicy const& overlaps_policy)
- {
- for (auto it = boost::begin(input); it != boost::end(input); ++it)
- {
- bool const lower_overlapping = overlaps_policy.apply(lower_box, **it);
- bool const upper_overlapping = overlaps_policy.apply(upper_box, **it);
- if (lower_overlapping && upper_overlapping)
- {
- exceeding.push_back(*it);
- }
- else if (lower_overlapping)
- {
- lower.push_back(*it);
- }
- else if (upper_overlapping)
- {
- upper.push_back(*it);
- }
- else
- {
- // Is nowhere. That is (since 1.58) possible, it might be
- // skipped by the OverlapsPolicy to enhance performance
- }
- }
- }
- template
- <
- typename Box,
- typename IteratorVector,
- typename ExpandPolicy
- >
- inline void expand_with_elements(Box& total, IteratorVector const& input,
- ExpandPolicy const& expand_policy)
- {
- for (auto const& it : input)
- {
- expand_policy.apply(total, *it);
- }
- }
- // Match forward_range with itself
- template <typename IteratorVector, typename VisitPolicy>
- inline bool handle_one(IteratorVector const& input, VisitPolicy& visitor)
- {
- if (boost::empty(input))
- {
- return true;
- }
- // Quadratic behaviour at lowest level (lowest quad, or all exceeding)
- for (auto it1 = boost::begin(input); it1 != boost::end(input); ++it1)
- {
- auto it2 = it1;
- for (++it2; it2 != boost::end(input); ++it2)
- {
- if (! visitor.apply(**it1, **it2))
- {
- return false; // Bail out if visitor returns false
- }
- }
- }
- return true;
- }
- // Match forward range 1 with forward range 2
- template
- <
- typename IteratorVector1,
- typename IteratorVector2,
- typename VisitPolicy
- >
- inline bool handle_two(IteratorVector1 const& input1,
- IteratorVector2 const& input2,
- VisitPolicy& visitor)
- {
- if (boost::empty(input1) || boost::empty(input2))
- {
- return true;
- }
- for (auto const& it1 : input1)
- {
- for (auto const& it2 : input2)
- {
- if (! visitor.apply(*it1, *it2))
- {
- return false; // Bail out if visitor returns false
- }
- }
- }
- return true;
- }
- template <typename IteratorVector>
- inline bool recurse_ok(IteratorVector const& input,
- std::size_t min_elements, std::size_t level)
- {
- return boost::size(input) >= min_elements
- && level < 100;
- }
- template <typename IteratorVector1, typename IteratorVector2>
- inline bool recurse_ok(IteratorVector1 const& input1,
- IteratorVector2 const& input2,
- std::size_t min_elements, std::size_t level)
- {
- return boost::size(input1) >= min_elements
- && recurse_ok(input2, min_elements, level);
- }
- template
- <
- typename IteratorVector1,
- typename IteratorVector2,
- typename IteratorVector3
- >
- inline bool recurse_ok(IteratorVector1 const& input1,
- IteratorVector2 const& input2,
- IteratorVector3 const& input3,
- std::size_t min_elements, std::size_t level)
- {
- return boost::size(input1) >= min_elements
- && recurse_ok(input2, input3, min_elements, level);
- }
- template <std::size_t Dimension, typename Box>
- class partition_two_ranges;
- template <std::size_t Dimension, typename Box>
- class partition_one_range
- {
- template <typename IteratorVector, typename ExpandPolicy>
- static inline Box get_new_box(IteratorVector const& input,
- ExpandPolicy const& expand_policy)
- {
- Box box;
- geometry::assign_inverse(box);
- expand_with_elements(box, input, expand_policy);
- return box;
- }
- template
- <
- typename IteratorVector,
- typename VisitPolicy,
- typename ExpandPolicy,
- typename OverlapsPolicy,
- typename VisitBoxPolicy
- >
- static inline bool next_level(Box const& box,
- IteratorVector const& input,
- std::size_t level, std::size_t min_elements,
- VisitPolicy& visitor,
- ExpandPolicy const& expand_policy,
- OverlapsPolicy const& overlaps_policy,
- VisitBoxPolicy& box_policy)
- {
- if (recurse_ok(input, min_elements, level))
- {
- return partition_one_range
- <
- 1 - Dimension,
- Box
- >::apply(box, input, level + 1, min_elements,
- visitor, expand_policy, overlaps_policy, box_policy);
- }
- else
- {
- return handle_one(input, visitor);
- }
- }
- // Function to switch to two forward ranges if there are
- // geometries exceeding the separation line
- template
- <
- typename IteratorVector,
- typename VisitPolicy,
- typename ExpandPolicy,
- typename OverlapsPolicy,
- typename VisitBoxPolicy
- >
- static inline bool next_level2(Box const& box,
- IteratorVector const& input1,
- IteratorVector const& input2,
- std::size_t level, std::size_t min_elements,
- VisitPolicy& visitor,
- ExpandPolicy const& expand_policy,
- OverlapsPolicy const& overlaps_policy,
- VisitBoxPolicy& box_policy)
- {
- if (recurse_ok(input1, input2, min_elements, level))
- {
- return partition_two_ranges
- <
- 1 - Dimension, Box
- >::apply(box, input1, input2, level + 1, min_elements,
- visitor, expand_policy, overlaps_policy,
- expand_policy, overlaps_policy, box_policy);
- }
- else
- {
- return handle_two(input1, input2, visitor);
- }
- }
- public :
- template
- <
- typename IteratorVector,
- typename VisitPolicy,
- typename ExpandPolicy,
- typename OverlapsPolicy,
- typename VisitBoxPolicy
- >
- static inline bool apply(Box const& box,
- IteratorVector const& input,
- std::size_t level,
- std::size_t min_elements,
- VisitPolicy& visitor,
- ExpandPolicy const& expand_policy,
- OverlapsPolicy const& overlaps_policy,
- VisitBoxPolicy& box_policy)
- {
- box_policy.apply(box, level);
- Box lower_box, upper_box;
- divide_box<Dimension>(box, lower_box, upper_box);
- IteratorVector lower, upper, exceeding;
- divide_into_subsets(lower_box, upper_box,
- input, lower, upper, exceeding,
- overlaps_policy);
- if (! boost::empty(exceeding))
- {
- // Get the box of exceeding-only
- Box const exceeding_box = get_new_box(exceeding, expand_policy);
- // Recursively do exceeding elements only, in next dimension they
- // will probably be less exceeding within the new box
- if (! (next_level(exceeding_box, exceeding, level, min_elements,
- visitor, expand_policy, overlaps_policy, box_policy)
- // Switch to two forward ranges, combine exceeding with
- // lower resp upper, but not lower/lower, upper/upper
- && next_level2(exceeding_box, exceeding, lower, level, min_elements,
- visitor, expand_policy, overlaps_policy, box_policy)
- && next_level2(exceeding_box, exceeding, upper, level, min_elements,
- visitor, expand_policy, overlaps_policy, box_policy)) )
- {
- return false; // Bail out if visitor returns false
- }
- }
- // Recursively call operation both parts
- return next_level(lower_box, lower, level, min_elements,
- visitor, expand_policy, overlaps_policy, box_policy)
- && next_level(upper_box, upper, level, min_elements,
- visitor, expand_policy, overlaps_policy, box_policy);
- }
- };
- template
- <
- std::size_t Dimension,
- typename Box
- >
- class partition_two_ranges
- {
- template
- <
- typename IteratorVector1,
- typename IteratorVector2,
- typename VisitPolicy,
- typename ExpandPolicy1,
- typename OverlapsPolicy1,
- typename ExpandPolicy2,
- typename OverlapsPolicy2,
- typename VisitBoxPolicy
- >
- static inline bool next_level(Box const& box,
- IteratorVector1 const& input1,
- IteratorVector2 const& input2,
- std::size_t level, std::size_t min_elements,
- VisitPolicy& visitor,
- ExpandPolicy1 const& expand_policy1,
- OverlapsPolicy1 const& overlaps_policy1,
- ExpandPolicy2 const& expand_policy2,
- OverlapsPolicy2 const& overlaps_policy2,
- VisitBoxPolicy& box_policy)
- {
- return partition_two_ranges
- <
- 1 - Dimension, Box
- >::apply(box, input1, input2, level + 1, min_elements,
- visitor, expand_policy1, overlaps_policy1,
- expand_policy2, overlaps_policy2, box_policy);
- }
- template <typename IteratorVector, typename ExpandPolicy>
- static inline Box get_new_box(IteratorVector const& input,
- ExpandPolicy const& expand_policy)
- {
- Box box;
- geometry::assign_inverse(box);
- expand_with_elements(box, input, expand_policy);
- return box;
- }
- template
- <
- typename IteratorVector1, typename IteratorVector2,
- typename ExpandPolicy1, typename ExpandPolicy2
- >
- static inline Box get_new_box(IteratorVector1 const& input1,
- IteratorVector2 const& input2,
- ExpandPolicy1 const& expand_policy1,
- ExpandPolicy2 const& expand_policy2)
- {
- Box box = get_new_box(input1, expand_policy1);
- expand_with_elements(box, input2, expand_policy2);
- return box;
- }
- public :
- template
- <
- typename IteratorVector1,
- typename IteratorVector2,
- typename VisitPolicy,
- typename ExpandPolicy1,
- typename OverlapsPolicy1,
- typename ExpandPolicy2,
- typename OverlapsPolicy2,
- typename VisitBoxPolicy
- >
- static inline bool apply(Box const& box,
- IteratorVector1 const& input1,
- IteratorVector2 const& input2,
- std::size_t level,
- std::size_t min_elements,
- VisitPolicy& visitor,
- ExpandPolicy1 const& expand_policy1,
- OverlapsPolicy1 const& overlaps_policy1,
- ExpandPolicy2 const& expand_policy2,
- OverlapsPolicy2 const& overlaps_policy2,
- VisitBoxPolicy& box_policy)
- {
- box_policy.apply(box, level);
- Box lower_box, upper_box;
- divide_box<Dimension>(box, lower_box, upper_box);
- IteratorVector1 lower1, upper1, exceeding1;
- IteratorVector2 lower2, upper2, exceeding2;
- divide_into_subsets(lower_box, upper_box,
- input1, lower1, upper1, exceeding1,
- overlaps_policy1);
- divide_into_subsets(lower_box, upper_box,
- input2, lower2, upper2, exceeding2,
- overlaps_policy2);
- if (! boost::empty(exceeding1))
- {
- // All exceeding from 1 with 2:
- if (recurse_ok(exceeding1, exceeding2, min_elements, level))
- {
- Box const exceeding_box = get_new_box(exceeding1, exceeding2,
- expand_policy1, expand_policy2);
- if (! next_level(exceeding_box, exceeding1, exceeding2, level,
- min_elements, visitor, expand_policy1, overlaps_policy1,
- expand_policy2, overlaps_policy2, box_policy))
- {
- return false; // Bail out if visitor returns false
- }
- }
- else
- {
- if (! handle_two(exceeding1, exceeding2, visitor))
- {
- return false; // Bail out if visitor returns false
- }
- }
- // All exceeding from 1 with lower and upper of 2:
- // (Check sizes of all three forward ranges to avoid recurse into
- // the same combinations again and again)
- if (recurse_ok(lower2, upper2, exceeding1, min_elements, level))
- {
- Box const exceeding_box = get_new_box(exceeding1, expand_policy1);
- if (! (next_level(exceeding_box, exceeding1, lower2, level,
- min_elements, visitor, expand_policy1, overlaps_policy1,
- expand_policy2, overlaps_policy2, box_policy)
- && next_level(exceeding_box, exceeding1, upper2, level,
- min_elements, visitor, expand_policy1, overlaps_policy1,
- expand_policy2, overlaps_policy2, box_policy)) )
- {
- return false; // Bail out if visitor returns false
- }
- }
- else
- {
- if (! (handle_two(exceeding1, lower2, visitor)
- && handle_two(exceeding1, upper2, visitor)) )
- {
- return false; // Bail out if visitor returns false
- }
- }
- }
- if (! boost::empty(exceeding2))
- {
- // All exceeding from 2 with lower and upper of 1:
- if (recurse_ok(lower1, upper1, exceeding2, min_elements, level))
- {
- Box const exceeding_box = get_new_box(exceeding2, expand_policy2);
- if (! (next_level(exceeding_box, lower1, exceeding2, level,
- min_elements, visitor, expand_policy1, overlaps_policy1,
- expand_policy2, overlaps_policy2, box_policy)
- && next_level(exceeding_box, upper1, exceeding2, level,
- min_elements, visitor, expand_policy1, overlaps_policy1,
- expand_policy2, overlaps_policy2, box_policy)) )
- {
- return false; // Bail out if visitor returns false
- }
- }
- else
- {
- if (! (handle_two(lower1, exceeding2, visitor)
- && handle_two(upper1, exceeding2, visitor)) )
- {
- return false; // Bail out if visitor returns false
- }
- }
- }
- if (recurse_ok(lower1, lower2, min_elements, level))
- {
- if (! next_level(lower_box, lower1, lower2, level,
- min_elements, visitor, expand_policy1, overlaps_policy1,
- expand_policy2, overlaps_policy2, box_policy) )
- {
- return false; // Bail out if visitor returns false
- }
- }
- else
- {
- if (! handle_two(lower1, lower2, visitor))
- {
- return false; // Bail out if visitor returns false
- }
- }
- if (recurse_ok(upper1, upper2, min_elements, level))
- {
- if (! next_level(upper_box, upper1, upper2, level,
- min_elements, visitor, expand_policy1, overlaps_policy1,
- expand_policy2, overlaps_policy2, box_policy) )
- {
- return false; // Bail out if visitor returns false
- }
- }
- else
- {
- if (! handle_two(upper1, upper2, visitor))
- {
- return false; // Bail out if visitor returns false
- }
- }
- return true;
- }
- };
- }} // namespace detail::partition
- template
- <
- typename Box,
- typename IncludePolicy1 = detail::partition::include_all_policy,
- typename IncludePolicy2 = detail::partition::include_all_policy
- >
- class partition
- {
- static const std::size_t default_min_elements = 16;
- template
- <
- typename IncludePolicy,
- typename ForwardRange,
- typename IteratorVector,
- typename ExpandPolicy
- >
- static inline void expand_to_range(ForwardRange const& forward_range,
- Box& total,
- IteratorVector& iterator_vector,
- ExpandPolicy const& expand_policy)
- {
- for (auto it = boost::begin(forward_range); it != boost::end(forward_range); ++it)
- {
- if (IncludePolicy::apply(*it))
- {
- expand_policy.apply(total, *it);
- iterator_vector.push_back(it);
- }
- }
- }
- public:
- template
- <
- typename ForwardRange,
- typename VisitPolicy,
- typename ExpandPolicy,
- typename OverlapsPolicy
- >
- static inline bool apply(ForwardRange const& forward_range,
- VisitPolicy& visitor,
- ExpandPolicy const& expand_policy,
- OverlapsPolicy const& overlaps_policy)
- {
- return apply(forward_range, visitor, expand_policy, overlaps_policy,
- default_min_elements, detail::partition::visit_no_policy());
- }
- template
- <
- typename ForwardRange,
- typename VisitPolicy,
- typename ExpandPolicy,
- typename OverlapsPolicy
- >
- static inline bool apply(ForwardRange const& forward_range,
- VisitPolicy& visitor,
- ExpandPolicy const& expand_policy,
- OverlapsPolicy const& overlaps_policy,
- std::size_t min_elements)
- {
- return apply(forward_range, visitor, expand_policy, overlaps_policy,
- min_elements, detail::partition::visit_no_policy());
- }
- template
- <
- typename ForwardRange,
- typename VisitPolicy,
- typename ExpandPolicy,
- typename OverlapsPolicy,
- typename VisitBoxPolicy
- >
- static inline bool apply(ForwardRange const& forward_range,
- VisitPolicy& visitor,
- ExpandPolicy const& expand_policy,
- OverlapsPolicy const& overlaps_policy,
- std::size_t min_elements,
- VisitBoxPolicy box_visitor)
- {
- using iterator_t = typename boost::range_iterator
- <
- ForwardRange const
- >::type;
- if (std::size_t(boost::size(forward_range)) > min_elements)
- {
- std::vector<iterator_t> iterator_vector;
- Box total;
- assign_inverse(total);
- expand_to_range<IncludePolicy1>(forward_range, total,
- iterator_vector, expand_policy);
- return detail::partition::partition_one_range
- <
- 0, Box
- >::apply(total, iterator_vector, 0, min_elements,
- visitor, expand_policy, overlaps_policy, box_visitor);
- }
- else
- {
- for (auto it1 = boost::begin(forward_range);
- it1 != boost::end(forward_range);
- ++it1)
- {
- auto it2 = it1;
- for (++it2; it2 != boost::end(forward_range); ++it2)
- {
- if (! visitor.apply(*it1, *it2))
- {
- return false; // Bail out if visitor returns false
- }
- }
- }
- }
- return true;
- }
- template
- <
- typename ForwardRange1,
- typename ForwardRange2,
- typename VisitPolicy,
- typename ExpandPolicy1,
- typename OverlapsPolicy1
- >
- static inline bool apply(ForwardRange1 const& forward_range1,
- ForwardRange2 const& forward_range2,
- VisitPolicy& visitor,
- ExpandPolicy1 const& expand_policy1,
- OverlapsPolicy1 const& overlaps_policy1)
- {
- return apply(forward_range1, forward_range2, visitor,
- expand_policy1, overlaps_policy1, expand_policy1, overlaps_policy1,
- default_min_elements, detail::partition::visit_no_policy());
- }
- template
- <
- typename ForwardRange1,
- typename ForwardRange2,
- typename VisitPolicy,
- typename ExpandPolicy1,
- typename OverlapsPolicy1,
- typename ExpandPolicy2,
- typename OverlapsPolicy2
- >
- static inline bool apply(ForwardRange1 const& forward_range1,
- ForwardRange2 const& forward_range2,
- VisitPolicy& visitor,
- ExpandPolicy1 const& expand_policy1,
- OverlapsPolicy1 const& overlaps_policy1,
- ExpandPolicy2 const& expand_policy2,
- OverlapsPolicy2 const& overlaps_policy2)
- {
- return apply(forward_range1, forward_range2, visitor,
- expand_policy1, overlaps_policy1, expand_policy2, overlaps_policy2,
- default_min_elements, detail::partition::visit_no_policy());
- }
- template
- <
- typename ForwardRange1,
- typename ForwardRange2,
- typename VisitPolicy,
- typename ExpandPolicy1,
- typename OverlapsPolicy1,
- typename ExpandPolicy2,
- typename OverlapsPolicy2
- >
- static inline bool apply(ForwardRange1 const& forward_range1,
- ForwardRange2 const& forward_range2,
- VisitPolicy& visitor,
- ExpandPolicy1 const& expand_policy1,
- OverlapsPolicy1 const& overlaps_policy1,
- ExpandPolicy2 const& expand_policy2,
- OverlapsPolicy2 const& overlaps_policy2,
- std::size_t min_elements)
- {
- return apply(forward_range1, forward_range2, visitor,
- expand_policy1, overlaps_policy1, expand_policy2, overlaps_policy2,
- min_elements, detail::partition::visit_no_policy());
- }
- template
- <
- typename ForwardRange1,
- typename ForwardRange2,
- typename VisitPolicy,
- typename ExpandPolicy1,
- typename OverlapsPolicy1,
- typename ExpandPolicy2,
- typename OverlapsPolicy2,
- typename VisitBoxPolicy
- >
- static inline bool apply(ForwardRange1 const& forward_range1,
- ForwardRange2 const& forward_range2,
- VisitPolicy& visitor,
- ExpandPolicy1 const& expand_policy1,
- OverlapsPolicy1 const& overlaps_policy1,
- ExpandPolicy2 const& expand_policy2,
- OverlapsPolicy2 const& overlaps_policy2,
- std::size_t min_elements,
- VisitBoxPolicy box_visitor)
- {
- using iterator1_t = typename boost::range_iterator
- <
- ForwardRange1 const
- >::type;
- using iterator2_t = typename boost::range_iterator
- <
- ForwardRange2 const
- >::type;
- if (std::size_t(boost::size(forward_range1)) > min_elements
- && std::size_t(boost::size(forward_range2)) > min_elements)
- {
- std::vector<iterator1_t> iterator_vector1;
- std::vector<iterator2_t> iterator_vector2;
- Box total;
- assign_inverse(total);
- expand_to_range<IncludePolicy1>(forward_range1, total,
- iterator_vector1, expand_policy1);
- expand_to_range<IncludePolicy2>(forward_range2, total,
- iterator_vector2, expand_policy2);
- return detail::partition::partition_two_ranges
- <
- 0, Box
- >::apply(total, iterator_vector1, iterator_vector2,
- 0, min_elements, visitor, expand_policy1,
- overlaps_policy1, expand_policy2, overlaps_policy2,
- box_visitor);
- }
- else
- {
- for (auto it1 = boost::begin(forward_range1);
- it1 != boost::end(forward_range1);
- ++it1)
- {
- for (auto it2 = boost::begin(forward_range2);
- it2 != boost::end(forward_range2);
- ++it2)
- {
- if (! visitor.apply(*it1, *it2))
- {
- return false; // Bail out if visitor returns false
- }
- }
- }
- }
- return true;
- }
- };
- }} // namespace boost::geometry
- #endif // BOOST_GEOMETRY_ALGORITHMS_DETAIL_PARTITION_HPP
|