boost/geometry
2
0
Fork 0
mirror of https://github.com/boostorg/geometry.git synced 2026-02-13 00:22:10 +00:00
Code Issues Projects Releases Wiki Activity

[algorithms][disjoint] add alternate implementation for

Browse Source 
disjoint(multipoint, multipoint) that uses BG's partition;
This commit is contained in:
Menelaos Karavelas
2014-07-31 15:40:58 +03:00
parent 79d8a2c59b
commit 9c96fc5e94
1 changed files with 76 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

76
include/boost/geometry/algorithms/detail/disjoint/multipoint_geometry.hpp
View File

@@ -10,6 +10,8 @@
#ifndef BOOST_GEOMETRY_ALGORITHMS_DETAIL_DISJOINT_MULTIPOINT_GEOMETRY_HPP
#define BOOST_GEOMETRY_ALGORITHMS_DETAIL_DISJOINT_MULTIPOINT_GEOMETRY_HPP
#include <algorithm>
#include <boost/assert.hpp>
#include <boost/range.hpp>
@@ -18,6 +20,14 @@
#include <boost/geometry/algorithms/detail/disjoint/point_geometry.hpp>
#include <boost/geometry/algorithms/detail/relate/less.hpp>
#if defined(BOOST_GEOMETRY_DISJOINT_USE_PARTITION)
#include <cstddef>
#include <boost/geometry/geometries/box.hpp>
#include <boost/geometry/algorithms/equals.hpp>
#include <boost/geometry/algorithms/expand.hpp>
#include <boost/geometry/algorithms/detail/partition.hpp>
#endif
#include <boost/geometry/algorithms/dispatch/disjoint.hpp>
@@ -75,6 +85,71 @@ struct multipoint_geometry
};
#if defined(BOOST_GEOMETRY_DISJOINT_USE_PARTITION)
template<typename MultiPoint1, typename MultiPoint2>
class multipoint_multipoint
{
private:
struct point_expand
{
template <typename Box, typename Point>
static inline void apply(Box& total, Point const& point)
{
geometry::expand(total, point);
}
};
struct overlaps_point
{
template <typename Box, typename Point>
static inline bool apply(Box const& box, Point const& point)
{
return !dispatch::disjoint<Point, Box>::apply(point, box);
}
};
struct two_point_visitor
{
std::size_t count;
two_point_visitor()
: count(0)
{}
template <typename Point1, typename Point2>
inline void apply(Point1 const& item1, Point2 const& item2)
{
if ( count == 0 && geometry::equals(item1, item2) )
{
++count;
}
}
};
public:
static inline bool apply(MultiPoint1 const& multipoint1,
MultiPoint2 const& multipoint2)
{
BOOST_ASSERT( boost::size(multipoint1) <= boost::size(multipoint2) );
typedef typename boost::range_value<MultiPoint1>::type point1_type;
two_point_visitor visitor;
geometry::partition
<
geometry::model::box
<
typename boost::range_value<MultiPoint1>::type
>,
point_expand,
overlaps_point
>::apply(multipoint1, multipoint2, visitor);
return visitor.count == 0;
}
};
#else
template<typename MultiPoint1, typename MultiPoint2>
struct multipoint_multipoint
{
@@ -108,6 +183,7 @@ struct multipoint_multipoint
return true;
}
};
#endif
}} // namespace detail::disjoint