123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172 |
- // Boost.Geometry (aka GGL, Generic Geometry Library)
- // Copyright (c) 2007-2015 Barend Gehrels, Amsterdam, the Netherlands.
- // Copyright (c) 2013-2015 Adam Wulkiewicz, Lodz, Poland.
- // This file was modified by Oracle on 2013, 2015, 2016.
- // Modifications copyright (c) 2013-2016, Oracle and/or its affiliates.
- // 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)
- #include "test_intersects.hpp"
- #include <boost/geometry/geometries/geometries.hpp>
- #include <boost/geometry/geometries/point_xy.hpp>
- #include <boost/geometry/util/rational.hpp>
- template <typename P1, typename P2>
- void test_all()
- {
- typedef bg::model::polygon<P1> polygon;
- typedef bg::model::ring<P1> ring;
- // intersect <=> ! disjoint (in most cases)
- // so most tests are done in disjoint test.
- // We only test compilation of a few cases.
- test_geometry<P1, bg::model::box<P2> >("POINT(1 1)", "BOX(0 0,2 2)", true);
- test_geometry<polygon, bg::model::box<P2> >(
- "POLYGON((1992 3240,1992 1440,3792 1800,3792 3240,1992 3240))",
- "BOX(1941 2066, 2055 2166)", true);
- test_geometry<ring, bg::model::box<P2> >(
- "POLYGON((1992 3240,1992 1440,3792 1800,3792 3240,1992 3240))",
- "BOX(1941 2066, 2055 2166)", true);
- test_geometry<polygon, bg::model::box<P2> >(
- "POLYGON((1941 2066,2055 2066,2055 2166,1941 2166))",
- "BOX(1941 2066, 2055 2166)", true);
- test_geometry<P1, bg::model::box<P2> >(
- "POINT(0 0)",
- "BOX(0 0,4 4)",
- true);
- }
- template <typename P>
- void test_all()
- {
- test_all<P, P>();
- }
- // Those tests won't pass for rational<> because numeric_limits<> isn't specialized for this type
- template <typename P>
- void test_additional()
- {
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(0 0,3 3)",
- "BOX(1 2,3 5)",
- true);
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(1 1,2 3)",
- "BOX(0 0,4 4)",
- true);
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(1 1,1 1)",
- "BOX(1 0,3 5)",
- true);
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(0 1,0 1)",
- "BOX(1 0,3 5)",
- false);
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(2 1,2 1)",
- "BOX(1 0,3 5)",
- true);
- test_geometry<bg::model::linestring<P>, bg::model::box<P> >(
- "LINESTRING(0 0,1 0,10 10)",
- "BOX(1 2,3 5)",
- true);
- test_geometry<bg::model::linestring<P>, bg::model::box<P> >(
- "LINESTRING(1 2)",
- "BOX(0 0,3 5)",
- true);
- // http://stackoverflow.com/questions/32457920/boost-rtree-of-box-gives-wrong-intersection-with-segment
- // http://lists.boost.org/geometry/2015/09/3476.php
- typedef bg::model::point<typename bg::coordinate_type<P>::type, 3, bg::cs::cartesian> point3d_t;
- test_geometry<bg::model::segment<point3d_t>, bg::model::box<point3d_t> >(
- "SEGMENT(2 1 0,2 1 10)",
- "BOX(0 0 0,10 0 10)",
- false);
- test_geometry<bg::model::segment<point3d_t>, bg::model::box<point3d_t> >(
- "SEGMENT(2 1 0,2 1 10)",
- "BOX(0 -5 0,10 0 10)",
- false);
- // and derived from the cases above
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(12 0,20 10)",
- "BOX(0 0,10 10)",
- false);
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(2 0,8 6)",
- "BOX(0 0,10 10)",
- true);
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(2 0,18 8)",
- "BOX(0 0,10 10)",
- true);
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(1 0,1 10)",
- "BOX(0 0,0 10)",
- false);
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(-1 0,-1 10)",
- "BOX(0 0,0 10)",
- false);
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(2 1,2 1)",
- "BOX(0 0,10 0)",
- false);
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(2 3,2 3)",
- "BOX(0 0,10 0)",
- false);
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(0 0,10 0)",
- "BOX(0 0,10 0)",
- true);
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(10 0,10 0)",
- "BOX(0 0,10 0)",
- true);
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(1 0,1 0)",
- "BOX(0 0,10 0)",
- true);
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(0 0,0 10)",
- "BOX(0 0,0 10)",
- true);
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(0 10,0 10)",
- "BOX(0 0,0 10)",
- true);
- test_geometry<bg::model::segment<P>, bg::model::box<P> >(
- "SEGMENT(0 1,0 1)",
- "BOX(0 0,0 10)",
- true);
- }
- int test_main( int , char* [] )
- {
- test_all<bg::model::d2::point_xy<float>, bg::model::point<double, 2, bg::cs::cartesian> >();
- test_all<bg::model::d2::point_xy<double> >();
- test_additional<bg::model::d2::point_xy<double> >();
- #if ! defined(BOOST_GEOMETRY_RESCALE_TO_ROBUST)
- test_all<bg::model::d2::point_xy<boost::rational<int> > >();
- #endif
- #if defined(HAVE_TTMATH)
- test_all<bg::model::d2::point_xy<ttmath_big> >();
- #endif
- return 0;
- }
|