diff options
Diffstat (limited to 'libs/geometry/test/algorithms/relational_operations/disjoint/disjoint_coverage.cpp')
| -rw-r--r-- | libs/geometry/test/algorithms/relational_operations/disjoint/disjoint_coverage.cpp | 1723 |
1 files changed, 1723 insertions, 0 deletions
diff --git a/libs/geometry/test/algorithms/relational_operations/disjoint/disjoint_coverage.cpp b/libs/geometry/test/algorithms/relational_operations/disjoint/disjoint_coverage.cpp new file mode 100644 index 000000000..4708fd6db --- /dev/null +++ b/libs/geometry/test/algorithms/relational_operations/disjoint/disjoint_coverage.cpp @@ -0,0 +1,1723 @@ +// Boost.Geometry (aka GGL, Generic Geometry Library) + +// Copyright (c) 2014, Oracle and/or its affiliates. + +// Licensed under the Boost Software License version 1.0. +// http://www.boost.org/users/license.html + +// Contributed and/or modified by Menelaos Karavelas, on behalf of Oracle + +#ifndef BOOST_TEST_MODULE +#define BOOST_TEST_MODULE test_disjoint_coverage +#endif + +// unit test to test disjoint for all geometry combinations + +#include <iostream> + +#include <boost/test/included/unit_test.hpp> + +#include <boost/geometry/core/tag.hpp> +#include <boost/geometry/core/tags.hpp> +#include <boost/geometry/multi/core/tags.hpp> + +#include <boost/geometry/strategies/strategies.hpp> + +#include <boost/geometry/io/wkt/read.hpp> +#include <boost/geometry/io/wkt/write.hpp> +#include <boost/geometry/multi/io/wkt/read.hpp> +#include <boost/geometry/multi/io/wkt/write.hpp> +#include <boost/geometry/io/dsv/write.hpp> + +#include <boost/geometry/geometries/point.hpp> +#include <boost/geometry/geometries/segment.hpp> +#include <boost/geometry/geometries/linestring.hpp> +#include <boost/geometry/geometries/ring.hpp> +#include <boost/geometry/geometries/polygon.hpp> + +#include <boost/geometry/multi/geometries/multi_point.hpp> +#include <boost/geometry/multi/geometries/multi_linestring.hpp> +#include <boost/geometry/multi/geometries/multi_polygon.hpp> + +#include <boost/geometry/algorithms/disjoint.hpp> + +#include <from_wkt.hpp> + + +#ifdef HAVE_TTMATH +#include <boost/geometry/extensions/contrib/ttmath_stub.hpp> +#endif + +namespace bg = ::boost::geometry; + +//============================================================================ + +struct test_disjoint +{ + template <typename Geometry1, typename Geometry2> + static inline void apply(std::string const& case_id, + Geometry1 const& geometry1, + Geometry2 const& geometry2, + bool expected_result) + { + bool result = bg::disjoint(geometry1, geometry2); + BOOST_CHECK_MESSAGE(result == expected_result, + "case ID: " << case_id << ", G1: " << bg::wkt(geometry1) + << ", G2: " << bg::wkt(geometry2) << " -> Expected: " + << expected_result << ", detected: " << result); + + result = bg::disjoint(geometry2, geometry1); + BOOST_CHECK_MESSAGE(result == expected_result, + "case ID: " << case_id << ", G1: " << bg::wkt(geometry2) + << ", G2: " << bg::wkt(geometry1) << " -> Expected: " + << expected_result << ", detected: " << result); + +#ifdef BOOST_GEOMETRY_TEST_DEBUG + std::cout << "case ID: " << case_id << "; G1 - G2: "; + std::cout << bg::wkt(geometry1) << " - "; + std::cout << bg::wkt(geometry2) << std::endl; + std::cout << std::boolalpha; + std::cout << "expected/computed result: " + << expected_result << " / " << result << std::endl; + std::cout << std::endl; + std::cout << std::noboolalpha; +#endif + } +}; + +//============================================================================ + +// pointlike-pointlike geometries +template <typename P> +inline void test_point_point() +{ + typedef test_disjoint tester; + + tester::apply("p-p-01", + from_wkt<P>("POINT(0 0)"), + from_wkt<P>("POINT(0 0)"), + false); + + tester::apply("p-p-02", + from_wkt<P>("POINT(0 0)"), + from_wkt<P>("POINT(1 1)"), + true); +} + +template <typename P> +inline void test_point_multipoint() +{ + typedef bg::model::multi_point<P> MP; + + typedef test_disjoint tester; + + tester::apply("p-mp-01", + from_wkt<P>("POINT(0 0)"), + from_wkt<MP>("MULTIPOINT(0 0,1 1)"), + false); + + tester::apply("p-mp-02", + from_wkt<P>("POINT(0 0)"), + from_wkt<MP>("MULTIPOINT(1 1,2 2)"), + true); + + tester::apply("p-mp-03", + from_wkt<P>("POINT(0 0)"), + from_wkt<MP>("MULTIPOINT()"), + true); +} + +template <typename P> +inline void test_multipoint_multipoint() +{ + typedef bg::model::multi_point<P> MP; + + typedef test_disjoint tester; + + tester::apply("mp-mp-01", + from_wkt<MP>("MULTIPOINT(0 0,1 0)"), + from_wkt<MP>("MULTIPOINT(0 0,1 1)"), + false); + + tester::apply("mp-mp-02", + from_wkt<MP>("MULTIPOINT(0 0,1 0)"), + from_wkt<MP>("MULTIPOINT(1 1,2 2)"), + true); + + tester::apply("mp-mp-03", + from_wkt<MP>("MULTIPOINT()"), + from_wkt<MP>("MULTIPOINT(1 1,2 2)"), + true); + + tester::apply("mp-mp-04", + from_wkt<MP>("MULTIPOINT(0 0,1 0)"), + from_wkt<MP>("MULTIPOINT()"), + true); +} + +//============================================================================ + +// pointlike-linear geometries +template <typename P> +inline void test_point_segment() +{ + typedef test_disjoint tester; + typedef bg::model::segment<P> S; + + tester::apply("p-s-01", + from_wkt<P>("POINT(0 0)"), + from_wkt<S>("SEGMENT(0 0,2 0)"), + false); + + tester::apply("p-s-02", + from_wkt<P>("POINT(1 0)"), + from_wkt<S>("SEGMENT(0 0,2 0)"), + false); + + tester::apply("p-s-03", + from_wkt<P>("POINT(1 1)"), + from_wkt<S>("SEGMENT(0 0,2 0)"), + true); + + tester::apply("p-s-04", + from_wkt<P>("POINT(3 0)"), + from_wkt<S>("SEGMENT(0 0,2 0)"), + true); + + tester::apply("p-s-05", + from_wkt<P>("POINT(-1 0)"), + from_wkt<S>("SEGMENT(0 0,2 0)"), + true); +} + +template <typename P> +inline void test_point_linestring() +{ + typedef bg::model::linestring<P> L; + + typedef test_disjoint tester; + + tester::apply("p-l-01", + from_wkt<P>("POINT(0 0)"), + from_wkt<L>("LINESTRING(0 0,2 2,4 4)"), + false); + + tester::apply("p-l-02", + from_wkt<P>("POINT(1 1)"), + from_wkt<L>("LINESTRING(0 0,2 2,4 4)"), + false); + + tester::apply("p-l-03", + from_wkt<P>("POINT(3 3)"), + from_wkt<L>("LINESTRING(0 0,2 2,4 4)"), + false); + + tester::apply("p-l-04", + from_wkt<P>("POINT(1 0)"), + from_wkt<L>("LINESTRING(0 0,2 2,4 4)"), + true); + + tester::apply("p-l-05", + from_wkt<P>("POINT(5 5)"), + from_wkt<L>("LINESTRING(0 0,2 2,4 4)"), + true); + + tester::apply("p-l-06", + from_wkt<P>("POINT(5 5)"), + from_wkt<L>("LINESTRING(0 0,2 2)"), + true); +} + +template <typename P> +inline void test_point_multilinestring() +{ + typedef bg::model::linestring<P> L; + typedef bg::model::multi_linestring<L> ML; + + typedef test_disjoint tester; + + tester::apply("p-ml-01", + from_wkt<P>("POINT(0 1)"), + from_wkt<ML>("MULTILINESTRING((0 0,2 2,4 4),(0 0,2 0,4 0))"), + true); + + tester::apply("p-ml-02", + from_wkt<P>("POINT(0 0)"), + from_wkt<ML>("MULTILINESTRING((0 0,2 2,4 4),(0 0,2 0,4 0))"), + false); + + tester::apply("p-ml-03", + from_wkt<P>("POINT(1 1)"), + from_wkt<ML>("MULTILINESTRING((0 0,2 2,4 4),(0 0,2 0,4 0))"), + false); + + tester::apply("p-ml-04", + from_wkt<P>("POINT(1 0)"), + from_wkt<ML>("MULTILINESTRING((0 0,2 2,4 4),(0 0,2 0,4 0))"), + false); + + tester::apply("p-ml-05", + from_wkt<P>("POINT(0 0)"), + from_wkt<ML>("MULTILINESTRING((1 1,2 2,4 4),(3 0,4 0))"), + true); + + tester::apply("p-ml-06", + from_wkt<P>("POINT(0 0)"), + from_wkt<ML>("MULTILINESTRING((1 1,2 2,4 4),(0 0,4 0))"), + false); + + tester::apply("p-ml-07", + from_wkt<P>("POINT(0 0)"), + from_wkt<ML>("MULTILINESTRING((1 1,2 2,4 4),(-1 0,4 0))"), + false); +} + +template <typename P> +inline void test_multipoint_segment() +{ + typedef test_disjoint tester; + typedef bg::model::multi_point<P> MP; + typedef bg::model::segment<P> S; + + tester::apply("mp-s-01", + from_wkt<MP>("MULTIPOINT(0 0,1 1)"), + from_wkt<S>("SEGMENT(0 0,2 0)"), + false); + + tester::apply("mp-s-02", + from_wkt<MP>("MULTIPOINT(1 0,1 1)"), + from_wkt<S>("SEGMENT(0 0,2 0)"), + false); + + tester::apply("mp-s-03", + from_wkt<MP>("MULTIPOINT(1 1,2 2)"), + from_wkt<S>("SEGMENT(0 0,2 0)"), + true); + + tester::apply("mp-s-04", + from_wkt<MP>("MULTIPOINT()"), + from_wkt<S>("SEGMENT(0 0,2 0)"), + true); + + tester::apply("mp-s-05", + from_wkt<MP>("MULTIPOINT(3 0,4 0)"), + from_wkt<S>("SEGMENT(0 0,2 0)"), + true); + + tester::apply("mp-s-06", + from_wkt<MP>("MULTIPOINT(1 0,4 0)"), + from_wkt<S>("SEGMENT(0 0,2 0)"), + false); + + // segments that degenerate to a point + tester::apply("mp-s-07", + from_wkt<MP>("MULTIPOINT(1 1,2 2)"), + from_wkt<S>("SEGMENT(0 0,0 0)"), + true); + + tester::apply("mp-s-08", + from_wkt<MP>("MULTIPOINT(1 1,2 2)"), + from_wkt<S>("SEGMENT(1 1,1 1)"), + false); +} + +template <typename P> +inline void test_multipoint_linestring() +{ + typedef bg::model::multi_point<P> MP; + typedef bg::model::linestring<P> L; + + typedef test_disjoint tester; + + tester::apply("mp-l-01", + from_wkt<MP>("MULTIPOINT(0 0,1 0)"), + from_wkt<L>("LINESTRING(0 0,2 2,4 4)"), + false); + + tester::apply("mp-l-02", + from_wkt<MP>("MULTIPOINT(1 0,1 1)"), + from_wkt<L>("LINESTRING(0 0,2 2,4 4)"), + false); + + tester::apply("mp-l-03", + from_wkt<MP>("MULTIPOINT(1 0,3 3)"), + from_wkt<L>("LINESTRING(0 0,2 2,4 4)"), + false); + + tester::apply("mp-l-04", + from_wkt<MP>("MULTIPOINT(1 0,2 0)"), + from_wkt<L>("LINESTRING(0 0,2 2,4 4)"), + true); + + tester::apply("mp-l-05", + from_wkt<MP>("MULTIPOINT(-1 -1,2 0)"), + from_wkt<L>("LINESTRING(0 0,2 2,4 4)"), + true); + + tester::apply("mp-l-06", + from_wkt<MP>("MULTIPOINT(-1 -1,2 0)"), + from_wkt<L>("LINESTRING(1 0,3 0)"), + false); + + tester::apply("mp-l-07", + from_wkt<MP>("MULTIPOINT(-1 -1,2 0)"), + from_wkt<L>("LINESTRING(1 0,3 0)"), + false); +} + +template <typename P> +inline void test_multipoint_multilinestring() +{ + typedef bg::model::multi_point<P> MP; + typedef bg::model::linestring<P> L; + typedef bg::model::multi_linestring<L> ML; + + typedef test_disjoint tester; + + tester::apply("mp-ml-01", + from_wkt<MP>("MULTIPOINT(0 1,0 2)"), + from_wkt<ML>("MULTILINESTRING((0 0,2 2,4 4),(0 0,2 0,4 0))"), + true); + + tester::apply("mp-ml-02", + from_wkt<MP>("MULTIPOINT(0 0,1 0)"), + from_wkt<ML>("MULTILINESTRING((0 0,2 2,4 4),(0 0,2 0,4 0))"), + false); + + tester::apply("mp-ml-03", + from_wkt<MP>("MULTIPOINT(0 1,1 1)"), + from_wkt<ML>("MULTILINESTRING((0 0,2 2,4 4),(0 0,2 0,4 0))"), + false); + + tester::apply("mp-ml-04", + from_wkt<MP>("MULTIPOINT(0 1,1 0)"), + from_wkt<ML>("MULTILINESTRING((0 0,2 2,4 4),(0 0,2 0,4 0))"), + false); +} + +//============================================================================ + +// pointlike-areal geometries +template <typename P> +inline void test_point_box() +{ + typedef test_disjoint tester; + typedef bg::model::box<P> B; + + tester::apply("p-b-01", + from_wkt<P>("POINT(0 0)"), + from_wkt<B>("BOX(0 0,1 1)"), + false); + + tester::apply("p-b-02", + from_wkt<P>("POINT(2 2)"), + from_wkt<B>("BOX(0 0,1 0)"), + true); +} + +template <typename P> +inline void test_point_ring() +{ + typedef bg::model::ring<P, false, false> R; // ccw, open + + typedef test_disjoint tester; + + tester::apply("p-r-01", + from_wkt<P>("POINT(0 0)"), + from_wkt<R>("POLYGON((0 0,1 0,0 1))"), + false); + + tester::apply("p-r-02", + from_wkt<P>("POINT(1 1)"), + from_wkt<R>("POLYGON((0 0,1 0,0 1))"), + true); +} + +template <typename P> +inline void test_point_polygon() +{ + typedef bg::model::polygon<P, false, false> PL; // ccw, open + + typedef test_disjoint tester; + + tester::apply("p-pg-01", + from_wkt<P>("POINT(0 0)"), + from_wkt<PL>("POLYGON((0 0,1 0,0 1))"), + false); + + tester::apply("p-pg-02", + from_wkt<P>("POINT(1 1)"), + from_wkt<PL>("POLYGON((0 0,1 0,0 1))"), + true); +} + +template <typename P> +inline void test_point_multipolygon() +{ + typedef bg::model::polygon<P, false, false> PL; // ccw, open + typedef bg::model::multi_polygon<PL> MPL; + + typedef test_disjoint tester; + + tester::apply("p-mpg-01", + from_wkt<P>("POINT(0 0)"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,1 0,0 1)),((2 0,3 0,2 1)))"), + false); + + tester::apply("p-mpg-02", + from_wkt<P>("POINT(1 1)"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,1 0,0 1)),((2 0,3 0,2 1)))"), + true); +} + +template <typename P> +inline void test_multipoint_box() +{ + typedef test_disjoint tester; + typedef bg::model::multi_point<P> MP; + typedef bg::model::box<P> B; + + tester::apply("mp-b-01", + from_wkt<MP>("MULTIPOINT(0 0,1 1)"), + from_wkt<B>("BOX(0 0,2 2)"), + false); + + tester::apply("mp-b-02", + from_wkt<MP>("MULTIPOINT(1 1,3 3)"), + from_wkt<B>("BOX(0 0,2 2)"), + false); + + tester::apply("mp-b-03", + from_wkt<MP>("MULTIPOINT(3 3,4 4)"), + from_wkt<B>("BOX(0 0,2 2)"), + true); + + tester::apply("mp-b-04", + from_wkt<MP>("MULTIPOINT()"), + from_wkt<B>("BOX(0 0,2 2)"), + true); +} + +template <typename P> +inline void test_multipoint_ring() +{ + typedef bg::model::multi_point<P> MP; + typedef bg::model::ring<P, false, false> R; // ccw, open + + typedef test_disjoint tester; + + tester::apply("mp-r-01", + from_wkt<MP>("MULTIPOINT(0 0,1 0)"), + from_wkt<R>("POLYGON((0 0,1 0,0 1))"), + false); + + tester::apply("mp-r-02", + from_wkt<MP>("MULTIPOINT(1 0,1 1)"), + from_wkt<R>("POLYGON((0 0,1 0,0 1))"), + false); + + tester::apply("mp-r-03", + from_wkt<MP>("MULTIPOINT(1 1,2 2)"), + from_wkt<R>("POLYGON((0 0,1 0,0 1))"), + true); +} + +template <typename P> +inline void test_multipoint_polygon() +{ + typedef bg::model::multi_point<P> MP; + typedef bg::model::polygon<P, false, false> PL; // ccw, open + + typedef test_disjoint tester; + + tester::apply("mp-pg-01", + from_wkt<MP>("MULTIPOINT(0 0,1 0)"), + from_wkt<PL>("POLYGON(((0 0,1 0,0 1)))"), + false); + + tester::apply("mp-pg-02", + from_wkt<MP>("MULTIPOINT(0 0,2 0)"), + from_wkt<PL>("POLYGON(((0 0,1 0,0 1)))"), + false); + + tester::apply("mp-pg-03", + from_wkt<MP>("MULTIPOINT(1 1,2 0)"), + from_wkt<PL>("POLYGON(((0 0,1 0,0 1)))"), + true); + + tester::apply("mp-pg-04", + from_wkt<MP>("MULTIPOINT(1 1,2 3)"), + from_wkt<PL>("POLYGON(((0 0,1 0,0 1)))"), + true); +} + +template <typename P> +inline void test_multipoint_multipolygon() +{ + typedef bg::model::multi_point<P> MP; + typedef bg::model::polygon<P, false, false> PL; // ccw, open + typedef bg::model::multi_polygon<PL> MPL; + + typedef test_disjoint tester; + + tester::apply("mp-mp-01", + from_wkt<MP>("MULTIPOINT(0 0,2 0)"), + from_wkt<MPL>("MULTIPOLYGON((0 0,1 0,0 1)),(2 0,3 0,2 1))"), + false); + + tester::apply("mp-mp-02", + from_wkt<MP>("MULTIPOINT(0 0,1 0)"), + from_wkt<MPL>("MULTIPOLYGON((0 0,1 0,0 1)),(2 0,3 0,2 1))"), + false); + + tester::apply("mp-mp-03", + from_wkt<MP>("MULTIPOINT(1 1,2 0)"), + from_wkt<MPL>("MULTIPOLYGON((0 0,1 0,0 1)),(2 0,3 0,2 1))"), + false); + + tester::apply("mp-mp-04", + from_wkt<MP>("MULTIPOINT(1 1,2 3)"), + from_wkt<MPL>("MULTIPOLYGON((0 0,1 0,0 1)),(2 0,3 0,2 1))"), + true); +} + +//============================================================================ + +// linear-linear geometries +template <typename P> +inline void test_segment_segment() +{ + typedef bg::model::segment<P> S; + + typedef test_disjoint tester; + + tester::apply("s-s-01", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<S>("SEGMENT(0 0,0 2)"), + false); + + tester::apply("s-s-02", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<S>("SEGMENT(2 0,3 0)"), + false); + + tester::apply("s-s-03", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<S>("SEGMENT(1 0,3 0)"), + false); + + tester::apply("s-s-04", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<S>("SEGMENT(1 0,1 1)"), + false); + + tester::apply("s-s-05", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<S>("SEGMENT(1 1,2 2)"), + true); + + tester::apply("s-s-06", + from_wkt<S>("SEGMENT(0 0,1 1)"), + from_wkt<S>("SEGMENT(1 1,1 1)"), + false); + + tester::apply("s-s-07", + from_wkt<S>("SEGMENT(0 0,1 1)"), + from_wkt<S>("SEGMENT(2 2,2 2)"), + true); + + tester::apply("s-s-08", + from_wkt<S>("SEGMENT(0 0,1 1)"), + from_wkt<S>("SEGMENT(2 2,3 3)"), + true); +} + +template <typename P> +inline void test_linestring_segment() +{ + typedef bg::model::segment<P> S; + typedef bg::model::linestring<P> L; + + typedef test_disjoint tester; + + tester::apply("l-s-01", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<L>("LINESTRING(0 0,0 2)"), + false); + + tester::apply("l-s-02", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<L>("LINESTRING(2 0,3 0)"), + false); + + tester::apply("l-s-03", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<L>("LINESTRING(1 0,3 0)"), + false); + + tester::apply("l-s-04", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<L>("LINESTRING(1 0,1 1)"), + false); + + tester::apply("l-s-05", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<L>("LINESTRING(1 1,2 2)"), + true); + + tester::apply("l-s-06", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<L>("LINESTRING(1 1,1 1,2 2)"), + true); + + tester::apply("l-s-07", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<L>("LINESTRING(1 0,1 0,1 1,2 2)"), + false); + + tester::apply("l-s-08", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<L>("LINESTRING(1 0,1 0,3 0)"), + false); + + tester::apply("l-s-09", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<L>("LINESTRING(3 0,3 0,4 0)"), + true); + + tester::apply("l-s-10", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<L>("LINESTRING(3 0,3 0)"), + true); + + tester::apply("l-s-11", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<L>("LINESTRING(-1 0,-1 0)"), + true); + + tester::apply("l-s-12", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<L>("LINESTRING(1 0,1 0)"), + false); + + tester::apply("l-s-13", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<L>("LINESTRING(1 1,1 1)"), + true); +} + +template <typename P> +inline void test_multilinestring_segment() +{ + typedef bg::model::segment<P> S; + typedef bg::model::linestring<P> L; + typedef bg::model::multi_linestring<L> ML; + + typedef test_disjoint tester; + + tester::apply("s-ml-01", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<ML>("MULTILINESTRING((0 0,0 2))"), + false); + + tester::apply("s-ml-02", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<ML>("MULTILINESTRING((2 0,3 0))"), + false); + + tester::apply("s-ml-03", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<ML>("MULTILINESTRING((1 0,3 0))"), + false); + + tester::apply("s-ml-04", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<ML>("MULTILINESTRING((1 0,1 1))"), + false); + + tester::apply("s-ml-05", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<ML>("MULTILINESTRING((1 1,2 2))"), + true); + + tester::apply("s-ml-06", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<ML>("MULTILINESTRING((1 1,2 2),(3 3,3 3))"), + true); + + tester::apply("s-ml-07", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<ML>("MULTILINESTRING((1 1,2 2),(1 0,1 0))"), + false); + + tester::apply("s-ml-08", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<ML>("MULTILINESTRING((1 1,2 2),(3 0,3 0))"), + true); +} + +template <typename P> +inline void test_linestring_linestring() +{ + typedef bg::model::linestring<P> L; + + typedef test_disjoint tester; + + tester::apply("l-l-01", + from_wkt<L>("LINESTRING(0 0,2 0)"), + from_wkt<L>("LINESTRING(0 0,0 2)"), + false); + + tester::apply("l-l-02", + from_wkt<L>("LINESTRING(0 0,2 0)"), + from_wkt<L>("LINESTRING(2 0,3 0)"), + false); + + tester::apply("l-l-03", + from_wkt<L>("LINESTRING(0 0,2 0)"), + from_wkt<L>("LINESTRING(1 0,3 0)"), + false); + + tester::apply("l-l-04", + from_wkt<L>("LINESTRING(0 0,2 0)"), + from_wkt<L>("LINESTRING(1 0,1 1)"), + false); + + tester::apply("l-l-05", + from_wkt<L>("LINESTRING(0 0,2 0)"), + from_wkt<L>("LINESTRING(1 1,2 2)"), + true); +} + +template <typename P> +inline void test_linestring_multilinestring() +{ + typedef bg::model::linestring<P> L; + typedef bg::model::multi_linestring<L> ML; + + typedef test_disjoint tester; + + tester::apply("l-ml-01", + from_wkt<L>("LINESTRING(0 0,2 0)"), + from_wkt<ML>("MULTILINESTRING((0 0,0 2))"), + false); + + tester::apply("l-ml-02", + from_wkt<L>("LINESTRING(0 0,2 0)"), + from_wkt<ML>("MULTILINESTRING((2 0,3 0))"), + false); + + tester::apply("l-ml-03", + from_wkt<L>("LINESTRING(0 0,2 0)"), + from_wkt<ML>("MULTILINESTRING((1 0,3 0))"), + false); + + tester::apply("l-ml-04", + from_wkt<L>("LINESTRING(0 0,2 0)"), + from_wkt<ML>("MULTILINESTRING((1 0,1 1))"), + false); + + tester::apply("l-ml-05", + from_wkt<L>("LINESTRING(0 0,2 0)"), + from_wkt<ML>("MULTILINESTRING((1 1,2 2))"), + true); +} + +template <typename P> +inline void test_multilinestring_multilinestring() +{ + typedef bg::model::linestring<P> L; + typedef bg::model::multi_linestring<L> ML; + + typedef test_disjoint tester; + + tester::apply("ml-ml-01", + from_wkt<ML>("MULTILINESTRING((0 0,2 0))"), + from_wkt<ML>("MULTILINESTRING((0 0,0 2))"), + false); + + tester::apply("ml-ml-02", + from_wkt<ML>("MULTILINESTRING((0 0,2 0))"), + from_wkt<ML>("MULTILINESTRING((2 0,3 0))"), + false); + + tester::apply("ml-ml-03", + from_wkt<ML>("MULTILINESTRING((0 0,2 0))"), + from_wkt<ML>("MULTILINESTRING((1 0,3 0))"), + false); + + tester::apply("ml-ml-04", + from_wkt<ML>("MULTILINESTRING((0 0,2 0))"), + from_wkt<ML>("MULTILINESTRING((1 0,1 1))"), + false); + + tester::apply("ml-ml-05", + from_wkt<ML>("MULTILINESTRING((0 0,2 0))"), + from_wkt<ML>("MULTILINESTRING((1 1,2 2))"), + true); +} + +//============================================================================ + +// linear-areal geometries +template <typename P> +inline void test_segment_box() +{ + typedef bg::model::segment<P> S; + typedef bg::model::box<P> B; + + typedef test_disjoint tester; + + tester::apply("s-b-01", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<B>("BOX(0 0,2 2)"), + false); + + tester::apply("s-b-02", + from_wkt<S>("SEGMENT(1 1,3 3)"), + from_wkt<B>("BOX(0 0,2 2)"), + false); + + tester::apply("s-b-03", + from_wkt<S>("SEGMENT(2 2,3 3)"), + from_wkt<B>("BOX(0 0,2 2)"), + false); + + tester::apply("s-b-04", + from_wkt<S>("SEGMENT(4 4,3 3)"), + from_wkt<B>("BOX(0 0,2 2)"), + true); + + tester::apply("s-b-05", + from_wkt<S>("SEGMENT(0 4,4 4)"), + from_wkt<B>("BOX(0 0,2 2)"), + true); + + tester::apply("s-b-06", + from_wkt<S>("SEGMENT(4 0,4 4)"), + from_wkt<B>("BOX(0 0,2 2)"), + true); + + tester::apply("s-b-07", + from_wkt<S>("SEGMENT(0 -2,0 -1)"), + from_wkt<B>("BOX(0 0,1 1)"), + true); + + tester::apply("s-b-08", + from_wkt<S>("SEGMENT(-2 -2,-2 -1)"), + from_wkt<B>("BOX(0 0,1 1)"), + true); + + tester::apply("s-b-09", + from_wkt<S>("SEGMENT(-2 -2,-2 -2)"), + from_wkt<B>("BOX(0 0,1 1)"), + true); + + tester::apply("s-b-10", + from_wkt<S>("SEGMENT(-2 0,-2 0)"), + from_wkt<B>("BOX(0 0,1 1)"), + true); + + tester::apply("s-b-11", + from_wkt<S>("SEGMENT(0 -2,0 -2)"), + from_wkt<B>("BOX(0 0,1 1)"), + true); + + tester::apply("s-b-12", + from_wkt<S>("SEGMENT(-2 0,-1 0)"), + from_wkt<B>("BOX(0 0,1 1)"), + true); + + // segment degenerates to a point + tester::apply("s-b-13", + from_wkt<S>("SEGMENT(0 0,0 0)"), + from_wkt<B>("BOX(0 0,1 1)"), + false); + + tester::apply("s-b-14", + from_wkt<S>("SEGMENT(1 1,1 1)"), + from_wkt<B>("BOX(0 0,2 2)"), + false); + + tester::apply("s-b-15", + from_wkt<S>("SEGMENT(2 2,2 2)"), + from_wkt<B>("BOX(0 0,2 2)"), + false); + + tester::apply("s-b-16", + from_wkt<S>("SEGMENT(2 0,2 0)"), + from_wkt<B>("BOX(0 0,2 2)"), + false); + + tester::apply("s-b-17", + from_wkt<S>("SEGMENT(0 2,0 2)"), + from_wkt<B>("BOX(0 0,2 2)"), + false); + + tester::apply("s-b-18", + from_wkt<S>("SEGMENT(2 2,2 2)"), + from_wkt<B>("BOX(0 0,1 1)"), + true); +} + +template <typename P> +inline void test_segment_ring() +{ + typedef bg::model::segment<P> S; + typedef bg::model::ring<P, false, false> R; // ccw, open + + typedef test_disjoint tester; + + tester::apply("s-r-01", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<R>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("s-r-02", + from_wkt<S>("SEGMENT(1 0,3 3)"), + from_wkt<R>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("s-r-03", + from_wkt<S>("SEGMENT(1 1,3 3)"), + from_wkt<R>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("s-r-04", + from_wkt<S>("SEGMENT(2 2,3 3)"), + from_wkt<R>("POLYGON((0 0,2 0,0 2))"), + true); +} + +template <typename P> +inline void test_segment_polygon() +{ + typedef bg::model::segment<P> S; + typedef bg::model::polygon<P, false, false> PL; // ccw, open + + typedef test_disjoint tester; + + tester::apply("s-pg-01", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<PL>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("s-pg-02", + from_wkt<S>("SEGMENT(1 0,3 3)"), + from_wkt<PL>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("s-pg-03", + from_wkt<S>("SEGMENT(1 1,3 3)"), + from_wkt<PL>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("s-pg-04", + from_wkt<S>("SEGMENT(2 2,3 3)"), + from_wkt<PL>("POLYGON((0 0,2 0,0 2))"), + true); +} + +template <typename P> +inline void test_segment_multipolygon() +{ + typedef bg::model::segment<P> S; + typedef bg::model::polygon<P, false, false> PL; // ccw, open + typedef bg::model::multi_polygon<PL> MPL; + + typedef test_disjoint tester; + + tester::apply("s-mpg-01", + from_wkt<S>("SEGMENT(0 0,2 0)"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,0 2)))"), + false); + + tester::apply("s-mpg-02", + from_wkt<S>("SEGMENT(1 0,3 3)"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,0 2)))"), + false); + + tester::apply("s-mpg-03", + from_wkt<S>("SEGMENT(1 1,3 3)"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,0 2)))"), + false); + + tester::apply("s-mpg-04", + from_wkt<S>("SEGMENT(2 2,3 3)"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,0 2)))"), + true); +} + +template <typename P> +inline void test_linestring_box() +{ + typedef bg::model::linestring<P> L; + typedef bg::model::box<P> B; + + typedef test_disjoint tester; + + tester::apply("l-b-01", + from_wkt<L>("LINESTRING(0 0,2 0)"), + from_wkt<B>("BOX(0 0,2 2)"), + false); + + tester::apply("l-b-02", + from_wkt<L>("LINESTRING(1 1,3 3)"), + from_wkt<B>("BOX(0 0,2 2)"), + false); + + tester::apply("l-b-03", + from_wkt<L>("LINESTRING(2 2,3 3)"), + from_wkt<B>("BOX(0 0,2 2)"), + false); + + tester::apply("l-b-04", + from_wkt<L>("LINESTRING(4 4,3 3)"), + from_wkt<B>("BOX(0 0,2 2)"), + true); +} + +template <typename P> +inline void test_linestring_ring() +{ + typedef bg::model::linestring<P> L; + typedef bg::model::ring<P, false, false> R; // ccw, open + + typedef test_disjoint tester; + + tester::apply("l-r-01", + from_wkt<L>("LINESTRING(0 0,2 0)"), + from_wkt<R>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("l-r-02", + from_wkt<L>("LINESTRING(1 0,3 3)"), + from_wkt<R>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("l-r-03", + from_wkt<L>("LINESTRING(1 1,3 3)"), + from_wkt<R>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("l-r-04", + from_wkt<L>("LINESTRING(2 2,3 3)"), + from_wkt<R>("POLYGON((0 0,2 0,0 2))"), + true); +} + +template <typename P> +inline void test_linestring_polygon() +{ + typedef bg::model::linestring<P> L; + typedef bg::model::polygon<P, false, false> PL; // ccw, open + + typedef test_disjoint tester; + + tester::apply("l-pg-01", + from_wkt<L>("LINESTRING(0 0,2 0)"), + from_wkt<PL>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("l-pg-02", + from_wkt<L>("LINESTRING(1 0,3 3)"), + from_wkt<PL>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("l-pg-03", + from_wkt<L>("LINESTRING(1 1,3 3)"), + from_wkt<PL>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("l-pg-04", + from_wkt<L>("LINESTRING(2 2,3 3)"), + from_wkt<PL>("POLYGON((0 0,2 0,0 2))"), + true); +} + +template <typename P> +inline void test_linestring_multipolygon() +{ + typedef bg::model::linestring<P> L; + typedef bg::model::polygon<P, false, false> PL; // ccw, open + typedef bg::model::multi_polygon<PL> MPL; + + typedef test_disjoint tester; + + tester::apply("l-mpg-01", + from_wkt<L>("LINESTRING(0 0,2 0)"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,0 2)))"), + false); + + tester::apply("l-mpg-02", + from_wkt<L>("LINESTRING(1 0,3 3)"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,0 2)))"), + false); + + tester::apply("l-mpg-03", + from_wkt<L>("LINESTRING(1 1,3 3)"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,0 2)))"), + false); + + tester::apply("l-mpg-04", + from_wkt<L>("LINESTRING(2 2,3 3)"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,0 2)))"), + true); +} + +template <typename P> +inline void test_multilinestring_box() +{ + typedef bg::model::linestring<P> L; + typedef bg::model::multi_linestring<L> ML; + typedef bg::model::box<P> B; + + typedef test_disjoint tester; + + tester::apply("ml-b-01", + from_wkt<ML>("MULTILINESTRING((0 0,2 0))"), + from_wkt<B>("BOX(0 0,2 2)"), + false); + + tester::apply("ml-b-02", + from_wkt<ML>("MULTILINESTRING((1 1,3 3))"), + from_wkt<B>("BOX(0 0,2 2)"), + false); + + tester::apply("ml-b-03", + from_wkt<ML>("MULTILINESTRING((2 2,3 3))"), + from_wkt<B>("BOX(0 0,2 2)"), + false); + + tester::apply("ml-b-04", + from_wkt<ML>("MULTILINESTRING((4 4,3 3))"), + from_wkt<B>("BOX(0 0,2 2)"), + true); +} + +template <typename P> +inline void test_multilinestring_ring() +{ + typedef bg::model::linestring<P> L; + typedef bg::model::multi_linestring<L> ML; + typedef bg::model::ring<P, false, false> R; // ccw, open + + typedef test_disjoint tester; + + tester::apply("ml-r-01", + from_wkt<ML>("MULTILINESTRING((0 0,2 0))"), + from_wkt<R>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("ml-r-02", + from_wkt<ML>("MULTILINESTRING((1 0,3 3))"), + from_wkt<R>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("ml-r-03", + from_wkt<ML>("MULTILINESTRING((1 1,3 3))"), + from_wkt<R>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("ml-r-04", + from_wkt<ML>("MULTILINESTRING((2 2,3 3))"), + from_wkt<R>("POLYGON((0 0,2 0,0 2))"), + true); +} + +template <typename P> +inline void test_multilinestring_polygon() +{ + typedef bg::model::linestring<P> L; + typedef bg::model::multi_linestring<L> ML; + typedef bg::model::polygon<P, false, false> PL; // ccw, open + + typedef test_disjoint tester; + + tester::apply("ml-pg-01", + from_wkt<ML>("MULTILINESTRING((0 0,2 0))"), + from_wkt<PL>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("ml-pg-02", + from_wkt<ML>("MULTILINESTRING((1 0,3 3))"), + from_wkt<PL>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("ml-pg-03", + from_wkt<ML>("MULTILINESTRING((1 1,3 3))"), + from_wkt<PL>("POLYGON((0 0,2 0,0 2))"), + false); + + tester::apply("ml-pg-04", + from_wkt<ML>("MULTILINESTRING((2 2,3 3))"), + from_wkt<PL>("POLYGON((0 0,2 0,0 2))"), + true); +} + +template <typename P> +inline void test_multilinestring_multipolygon() +{ + typedef bg::model::linestring<P> L; + typedef bg::model::multi_linestring<L> ML; + typedef bg::model::polygon<P, false, false> PL; // ccw, open + typedef bg::model::multi_polygon<PL> MPL; + + typedef test_disjoint tester; + + tester::apply("ml-mpg-01", + from_wkt<ML>("MULTILINESTRING((0 0,2 0))"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,0 2)))"), + false); + + tester::apply("ml-mpg-02", + from_wkt<ML>("MULTILINESTRING((1 0,3 3))"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,0 2)))"), + false); + + tester::apply("ml-mpg-03", + from_wkt<ML>("MULTILINESTRING((1 1,3 3))"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,0 2)))"), + false); + + tester::apply("ml-mpg-04", + from_wkt<ML>("MULTILINESTRING((2 2,3 3))"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,0 2)))"), + true); +} + +//============================================================================ + +// areal-areal geometries +template <typename P> +inline void test_box_box() +{ + typedef bg::model::box<P> B; + + typedef test_disjoint tester; + + tester::apply("b-b-01", + from_wkt<B>("BOX(2 2,3 3)"), + from_wkt<B>("BOX(0 0,2 2)"), + false); + + tester::apply("b-b-02", + from_wkt<B>("BOX(1 1,3 3)"), + from_wkt<B>("BOX(0 0,2 2)"), + false); + + tester::apply("b-b-03", + from_wkt<B>("BOX(3 3,4 4)"), + from_wkt<B>("BOX(0 0,2 2)"), + true); +} + +template <typename P> +inline void test_ring_box() +{ + typedef bg::model::box<P> B; + typedef bg::model::ring<P, false, false> R; // ccw, open + + typedef test_disjoint tester; + + tester::apply("r-b-01", + from_wkt<B>("BOX(2 2,3 3)"), + from_wkt<R>("POLYGON((0 0,2 0,2 2,0 2))"), + false); + + tester::apply("r-b-02", + from_wkt<B>("BOX(1 1,3 3)"), + from_wkt<R>("POLYGON((0 0,2 0,2 2,0 2))"), + false); + + tester::apply("r-b-03", + from_wkt<B>("BOX(3 3,4 4)"), + from_wkt<R>("POLYGON((0 0,2 0,2 2,0 2))"), + true); +} + +template <typename P> +inline void test_polygon_box() +{ + typedef bg::model::box<P> B; + typedef bg::model::polygon<P, false, false> PL; // ccw, open + + typedef test_disjoint tester; + + tester::apply("pg-b-01", + from_wkt<B>("BOX(2 2,3 3)"), + from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"), + false); + + tester::apply("pg-b-02", + from_wkt<B>("BOX(1 1,3 3)"), + from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"), + false); + + tester::apply("pg-b-03", + from_wkt<B>("BOX(3 3,4 4)"), + from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"), + true); +} + +template <typename P> +inline void test_multipolygon_box() +{ + typedef bg::model::box<P> B; + typedef bg::model::polygon<P, false, false> PL; // ccw, open + typedef bg::model::multi_polygon<PL> MPL; + + typedef test_disjoint tester; + + tester::apply("mpg-b-01", + from_wkt<B>("BOX(2 2,3 3)"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"), + false); + + tester::apply("mpg-b-02", + from_wkt<B>("BOX(1 1,3 3)"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"), + false); + + tester::apply("mpg-b-03", + from_wkt<B>("BOX(3 3,4 4)"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"), + true); +} + +template <typename P> +inline void test_ring_ring() +{ + typedef bg::model::ring<P, false, false> R; // ccw, open + + typedef test_disjoint tester; + + tester::apply("r-r-01", + from_wkt<R>("POLYGON((2 2,2 3,3 3,3 2))"), + from_wkt<R>("POLYGON((0 0,2 0,2 2,0 2))"), + false); + + tester::apply("r-r-02", + from_wkt<R>("POLYGON((1 1,1 3,3 3,3 1))"), + from_wkt<R>("POLYGON((0 0,2 0,2 2,0 2))"), + false); + + tester::apply("r-r-03", + from_wkt<R>("POLYGON((3 3,3 4,4 4,4 3))"), + from_wkt<R>("POLYGON((0 0,2 0,2 2,0 2))"), + true); +} + +template <typename P> +inline void test_polygon_ring() +{ + typedef bg::model::ring<P, false, false> R; // ccw, open + typedef bg::model::polygon<P, false, false> PL; // ccw, open + + typedef test_disjoint tester; + + tester::apply("pg-r-01", + from_wkt<R>("POLYGON((2 2,2 3,3 3,3 2))"), + from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"), + false); + + tester::apply("pg-r-02", + from_wkt<R>("POLYGON((1 1,1 3,3 3,3 1))"), + from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"), + false); + + tester::apply("pg-r-03", + from_wkt<R>("POLYGON((3 3,3 4,4 4,4 3))"), + from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"), + true); +} + +template <typename P> +inline void test_multipolygon_ring() +{ + typedef bg::model::ring<P, false, false> R; // ccw, open + typedef bg::model::polygon<P, false, false> PL; // ccw, open + typedef bg::model::multi_polygon<PL> MPL; + + typedef test_disjoint tester; + + tester::apply("mpg-r-01", + from_wkt<R>("POLYGON((2 2,2 3,3 3,3 2))"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"), + false); + + tester::apply("mpg-r-02", + from_wkt<R>("POLYGON((1 1,1 3,3 3,3 1))"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"), + false); + + tester::apply("mpg-r-03", + from_wkt<R>("POLYGON((3 3,3 4,4 4,4 3))"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"), + true); +} + +template <typename P> +inline void test_polygon_polygon() +{ + typedef bg::model::polygon<P, false, false> PL; // ccw, open + + typedef test_disjoint tester; + + tester::apply("pg-pg-01", + from_wkt<PL>("POLYGON((2 2,2 3,3 3,3 2))"), + from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"), + false); + + tester::apply("pg-pg-02", + from_wkt<PL>("POLYGON((1 1,1 3,3 3,3 1))"), + from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"), + false); + + tester::apply("pg-pg-03", + from_wkt<PL>("POLYGON((3 3,3 4,4 4,4 3))"), + from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"), + true); + + tester::apply("pg-pg-04", + from_wkt<PL>("POLYGON((0 0,9 0,9 9,0 9))"), + from_wkt<PL>("POLYGON((3 3,6 3,6 6,3 6))"), + false); + // polygon with a hole which entirely contains the other polygon + tester::apply("pg-pg-05", + from_wkt<PL>("POLYGON((0 0,9 0,9 9,0 9),(2 2,2 7,7 7,7 2))"), + from_wkt<PL>("POLYGON((3 3,6 3,6 6,3 6))"), + true); + // polygon with a hole, but the inner ring intersects the other polygon + tester::apply("pg-pg-06", + from_wkt<PL>("POLYGON((0 0,9 0,9 9,0 9),(3 2,3 7,7 7,7 2))"), + from_wkt<PL>("POLYGON((2 3,6 3,6 6,2 6))"), + false); + // polygon with a hole, but the other polygon is entirely contained + // between the inner and outer rings. + tester::apply("pg-pg-07", + from_wkt<PL>("POLYGON((0 0,9 0,9 9,0 9),(6 2,6 7,7 7,7 2))"), + from_wkt<PL>("POLYGON((3 3,5 3,5 6,3 6))"), + false); + // polygon with a hole and the outer ring of the other polygon lies + // between the inner and outer, but without touching either. + tester::apply("pg-pg-08", + from_wkt<PL>("POLYGON((0 0,9 0,9 9,0 9),(3 3,3 6,6 6,6 3))"), + from_wkt<PL>("POLYGON((2 2,7 2,7 7,2 7))"), + false); + + { + typedef bg::model::polygon<P> PL; // cw, closed + + // https://svn.boost.org/trac/boost/ticket/10647 + tester::apply("ticket-10647", + from_wkt<PL>("POLYGON((0 0, 0 5, 5 5, 5 0, 0 0)(1 1, 4 1, 4 4, 1 4, 1 1))"), + from_wkt<PL>("POLYGON((2 2, 2 3, 3 3, 3 2, 2 2))"), + true); + } +} + +template <typename P> +inline void test_polygon_multipolygon() +{ + typedef bg::model::polygon<P, false, false> PL; // ccw, open + typedef bg::model::multi_polygon<PL> MPL; + + typedef test_disjoint tester; + + tester::apply("pg-mpg-01", + from_wkt<PL>("POLYGON((2 2,2 3,3 3,3 2))"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"), + false); + + tester::apply("pg-mpg-02", + from_wkt<PL>("POLYGON((1 1,1 3,3 3,3 1))"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"), + false); + + tester::apply("pg-mpg-03", + from_wkt<PL>("POLYGON((3 3,3 4,4 4,4 3))"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"), + true); +} + +template <typename P> +inline void test_multipolygon_multipolygon() +{ + typedef bg::model::polygon<P, false, false> PL; // ccw, open + typedef bg::model::multi_polygon<PL> MPL; + + typedef test_disjoint tester; + + tester::apply("mpg-mpg-01", + from_wkt<MPL>("MULTIPOLYGON(((2 2,2 3,3 3,3 2)))"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"), + false); + + tester::apply("mpg-mpg-02", + from_wkt<MPL>("MULTIPOLYGON(((1 1,1 3,3 3,3 1)))"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"), + false); + + tester::apply("mpg-mpg-03", + from_wkt<MPL>("MULTIPOLYGON(((3 3,3 4,4 4,4 3)))"), + from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"), + true); +} + +//============================================================================ + +template <typename CoordinateType> +inline void test_pointlike_pointlike() +{ + typedef bg::model::point<CoordinateType, 2, bg::cs::cartesian> point_type; + + test_point_point<point_type>(); + test_point_multipoint<point_type>(); + + test_multipoint_multipoint<point_type>(); +} + +template <typename CoordinateType> +inline void test_pointlike_linear() +{ + typedef bg::model::point<CoordinateType, 2, bg::cs::cartesian> point_type; + + test_point_linestring<point_type>(); + test_point_multilinestring<point_type>(); + test_point_segment<point_type>(); + + // not implemented yet + // test_multipoint_linestring<point_type>(); + // test_multipoint_multilinestring<point_type>(); + test_multipoint_segment<point_type>(); +} + +template <typename CoordinateType> +inline void test_pointlike_areal() +{ + typedef bg::model::point<CoordinateType, 2, bg::cs::cartesian> point_type; + + test_point_polygon<point_type>(); + test_point_multipolygon<point_type>(); + test_point_ring<point_type>(); + test_point_box<point_type>(); + + // not implemented yet + // test_multipoint_polygon<point_type>(); + // test_multipoint_multipolygon<point_type>(); + // test_multipoint_ring<point_type>(); + test_multipoint_box<point_type>(); +} + +template <typename CoordinateType> +inline void test_linear_linear() +{ + typedef bg::model::point<CoordinateType, 2, bg::cs::cartesian> point_type; + + test_linestring_linestring<point_type>(); + test_linestring_multilinestring<point_type>(); + test_linestring_segment<point_type>(); + + test_multilinestring_multilinestring<point_type>(); + test_multilinestring_segment<point_type>(); + + test_segment_segment<point_type>(); +} + +template <typename CoordinateType> +inline void test_linear_areal() +{ + typedef bg::model::point<CoordinateType, 2, bg::cs::cartesian> point_type; + + test_segment_polygon<point_type>(); + test_segment_multipolygon<point_type>(); + test_segment_ring<point_type>(); + test_segment_box<point_type>(); + + test_linestring_polygon<point_type>(); + test_linestring_multipolygon<point_type>(); + test_linestring_ring<point_type>(); + test_linestring_box<point_type>(); + + test_multilinestring_polygon<point_type>(); + test_multilinestring_multipolygon<point_type>(); + test_multilinestring_ring<point_type>(); + test_multilinestring_box<point_type>(); +} + +template <typename CoordinateType> +inline void test_areal_areal() +{ + typedef bg::model::point<CoordinateType, 2, bg::cs::cartesian> point_type; + + test_polygon_polygon<point_type>(); + test_polygon_multipolygon<point_type>(); + test_polygon_ring<point_type>(); + test_polygon_box<point_type>(); + + test_multipolygon_multipolygon<point_type>(); + test_multipolygon_ring<point_type>(); + test_multipolygon_box<point_type>(); + + test_ring_ring<point_type>(); + test_ring_box<point_type>(); + + test_box_box<point_type>(); +} + +//============================================================================ + +BOOST_AUTO_TEST_CASE( test_pointlike_pointlike_all ) +{ + test_pointlike_pointlike<double>(); + test_pointlike_pointlike<int>(); +#ifdef HAVE_TTMATH + test_pointlike_pointlike<ttmath_big>(); +#endif +} + +BOOST_AUTO_TEST_CASE( test_pointlike_linear_all ) +{ + test_pointlike_linear<double>(); + test_pointlike_linear<int>(); +#ifdef HAVE_TTMATH + test_pointlike_linear<ttmath_big>(); +#endif +} + +BOOST_AUTO_TEST_CASE( test_pointlike_areal_all ) +{ + test_pointlike_areal<double>(); + test_pointlike_areal<int>(); +#ifdef HAVE_TTMATH + test_pointlike_areal<ttmath_big>(); +#endif +} + +BOOST_AUTO_TEST_CASE( test_linear_linear_all ) +{ + test_linear_linear<double>(); + test_linear_linear<int>(); +#ifdef HAVE_TTMATH + test_linear_linear<ttmath_big>(); +#endif +} + +BOOST_AUTO_TEST_CASE( test_linear_areal_all ) +{ + test_linear_areal<double>(); + test_linear_areal<int>(); +#ifdef HAVE_TTMATH + test_linear_areal<ttmath_big>(); +#endif +} + +BOOST_AUTO_TEST_CASE( test_areal_areal_all ) +{ + test_areal_areal<double>(); + test_areal_areal<int>(); +#ifdef HAVE_TTMATH + test_areal_areal<ttmath_big>(); +#endif +} |