123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523 |
- // Boost.Geometry (aka GGL, Generic Geometry Library)
- // Unit Test
- // Copyright (c) 2007-2014 Barend Gehrels, Amsterdam, the Netherlands.
- // Copyright (c) 2008-2014 Bruno Lalande, Paris, France.
- // Copyright (c) 2009-2014 Mateusz Loskot, London, UK.
- // Copyright (c) 2013-2014 Adam Wulkiewicz, Lodz, Poland.
- // This file was modified by Oracle on 2014, 2017.
- // Modifications copyright (c) 2014-2017, 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
- // Parts of Boost.Geometry are redesigned from Geodan's Geographic Library
- // (geolib/GGL), copyright (c) 1995-2010 Geodan, Amsterdam, the Netherlands.
- // 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 <sstream>
- #include <boost/mpl/if.hpp>
- #include <boost/type_traits/is_integral.hpp>
- #include <boost/type_traits/is_same.hpp>
- #include <geometry_test_common.hpp>
- #include <boost/geometry/algorithms/comparable_distance.hpp>
- #include <boost/geometry/algorithms/make.hpp>
- #include <boost/geometry/geometries/geometries.hpp>
- #include <boost/geometry/geometries/point_xy.hpp>
- #include <boost/geometry/io/wkt/read.hpp>
- #include <boost/geometry/strategies/strategies.hpp>
- template <typename P>
- void test_distance_result()
- {
- typedef typename bg::default_distance_result<P, P>::type distance_type;
- P p1 = bg::make<P>(0, 0);
- P p2 = bg::make<P>(3, 0);
- P p3 = bg::make<P>(0, 4);
- distance_type dr12 = bg::comparable_distance(p1, p2);
- distance_type dr13 = bg::comparable_distance(p1, p3);
- distance_type dr23 = bg::comparable_distance(p2, p3);
- BOOST_CHECK_CLOSE(dr12, 9.000, 0.001);
- BOOST_CHECK_CLOSE(dr13, 16.000, 0.001);
- BOOST_CHECK_CLOSE(dr23, 25.000, 0.001);
- }
- template <typename P>
- void test_distance_point()
- {
- P p1;
- bg::set<0>(p1, 1);
- bg::set<1>(p1, 1);
- P p2;
- bg::set<0>(p2, 2);
- bg::set<1>(p2, 2);
- typename bg::coordinate_type<P>::type d = bg::comparable_distance(p1, p2);
- BOOST_CHECK_CLOSE(d, 2.0, 0.001);
- }
- template <typename P>
- void test_distance_segment()
- {
- typedef typename bg::coordinate_type<P>::type coordinate_type;
- P s1 = bg::make<P>(2, 2);
- P s2 = bg::make<P>(3, 3);
- // Check points left, right, projected-left, projected-right, on segment
- P p1 = bg::make<P>(0, 0);
- P p2 = bg::make<P>(4, 4);
- P p3 = bg::make<P>(2.4, 2.6);
- P p4 = bg::make<P>(2.6, 2.4);
- P p5 = bg::make<P>(2.5, 2.5);
- bg::model::referring_segment<P const> const seg(s1, s2);
- coordinate_type d1 = bg::comparable_distance(p1, seg); BOOST_CHECK_CLOSE(d1, 8.0, 0.001);
- coordinate_type d2 = bg::comparable_distance(p2, seg); BOOST_CHECK_CLOSE(d2, 2.0, 0.001);
- coordinate_type d3 = bg::comparable_distance(p3, seg); BOOST_CHECK_CLOSE(d3, 0.02, 0.001);
- coordinate_type d4 = bg::comparable_distance(p4, seg); BOOST_CHECK_CLOSE(d4, 0.02, 0.001);
- coordinate_type d5 = bg::comparable_distance(p5, seg); BOOST_CHECK_CLOSE(d5, 0.0, 0.001);
- // Reverse case
- coordinate_type dr1 = bg::comparable_distance(seg, p1); BOOST_CHECK_CLOSE(dr1, d1, 0.001);
- coordinate_type dr2 = bg::comparable_distance(seg, p2); BOOST_CHECK_CLOSE(dr2, d2, 0.001);
- }
- template <typename P>
- void test_distance_linestring()
- {
- bg::model::linestring<P> points;
- points.push_back(bg::make<P>(1, 1));
- points.push_back(bg::make<P>(3, 3));
- P p = bg::make<P>(2, 1);
- typename bg::coordinate_type<P>::type d = bg::comparable_distance(p, points);
- BOOST_CHECK_CLOSE(d, 0.5, 0.001);
- p = bg::make<P>(5, 5);
- d = bg::comparable_distance(p, points);
- BOOST_CHECK_CLOSE(d, 8.0, 0.001);
- bg::model::linestring<P> line;
- line.push_back(bg::make<P>(1,1));
- line.push_back(bg::make<P>(2,2));
- line.push_back(bg::make<P>(3,3));
- p = bg::make<P>(5, 5);
- d = bg::comparable_distance(p, line);
- BOOST_CHECK_CLOSE(d, 8.0, 0.001);
- // Reverse case
- d = bg::comparable_distance(line, p);
- BOOST_CHECK_CLOSE(d, 8.0, 0.001);
- }
- template <typename P>
- void test_all()
- {
- test_distance_result<P>();
- test_distance_point<P>();
- test_distance_segment<P>();
- test_distance_linestring<P>();
- }
- template <typename T>
- void test_double_result_from_integer()
- {
- typedef bg::model::point<T, 2, bg::cs::cartesian> point_type;
- point_type point;
- // Check linestring
- bg::model::linestring<point_type> linestring;
- bg::read_wkt("POINT(2 2)", point);
- bg::read_wkt("LINESTRING(4 1,1 4)", linestring);
- double normal_distance = bg::distance(point, linestring);
- double comparable_distance = bg::comparable_distance(point, linestring);
- BOOST_CHECK_CLOSE(normal_distance, std::sqrt(0.5), 0.001);
- BOOST_CHECK_CLOSE(comparable_distance, 0.5, 0.001);
- // Check polygon
- bg::model::polygon<point_type> polygon;
- bg::read_wkt("POLYGON((0 0,1 9,8 1,0 0),(1 1,4 1,1 4,1 1))", polygon);
- normal_distance = bg::distance(point, polygon);
- comparable_distance = bg::comparable_distance(point, polygon);
- BOOST_CHECK_CLOSE(normal_distance, std::sqrt(0.5), 0.001);
- BOOST_CHECK_CLOSE(comparable_distance, 0.5, 0.001);
- }
- template <typename T>
- struct test_variant_different_default_strategy
- {
- static inline void apply()
- {
- typedef bg::model::point<T, 2, bg::cs::cartesian> point_type;
- typedef bg::model::segment<point_type> segment_type;
- typedef bg::model::box<point_type> box_type;
- typedef boost::variant<point_type, segment_type, box_type> variant_type;
- point_type point;
- bg::read_wkt("POINT(1 3)", point);
- segment_type seg;
- bg::read_wkt("LINESTRING(1 1,4 4)", seg);
- box_type box;
- bg::read_wkt("BOX(-1 -1,0 0)", box);
- variant_type v1, v2;
-
- BOOST_MPL_ASSERT((
- boost::is_same
- <
- typename bg::comparable_distance_result
- <
- variant_type, variant_type, bg::default_strategy
- >::type,
- typename bg::comparable_distance_result
- <
- point_type, point_type, bg::default_strategy
- >::type
- >
- ));
- // Default strategy
- v1 = point;
- v2 = point;
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2),
- bg::comparable_distance(point, point),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, point),
- bg::comparable_distance(point, point),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(point, v2),
- bg::comparable_distance(point, point),
- 0.0001);
- v1 = point;
- v2 = seg;
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2),
- bg::comparable_distance(point, seg),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, seg),
- bg::comparable_distance(point, seg),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(point, v2),
- bg::comparable_distance(point, seg), 0.0001);
- v1 = point;
- v2 = box;
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2),
- bg::comparable_distance(point, box),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, box),
- bg::comparable_distance(point, box),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(point, v2),
- bg::comparable_distance(point, box), 0.0001);
- }
- };
- template <typename T, typename ExpectedResultType = double>
- struct test_variant_same_default_strategy
- {
- static inline void apply()
- {
- typedef bg::model::point<T, 2, bg::cs::cartesian> point_type;
- typedef bg::model::segment<point_type> segment_type;
- typedef bg::model::linestring<point_type> linestring_type;
- typedef boost::variant
- <
- point_type, segment_type, linestring_type
- > variant_type;
- point_type point;
- bg::read_wkt("POINT(1 3)", point);
- segment_type seg;
- bg::read_wkt("LINESTRING(1 1,4 4)", seg);
- linestring_type linestring;
- bg::read_wkt("LINESTRING(-1 -1,-1 0,0 0,0 -1,-1 -1)", linestring);
- variant_type v1, v2;
-
- BOOST_MPL_ASSERT((
- boost::is_same
- <
- typename bg::comparable_distance_result
- <
- variant_type, variant_type, bg::default_strategy
- >::type,
- ExpectedResultType
- >
- ));
- BOOST_MPL_ASSERT((
- boost::is_same
- <
- typename bg::comparable_distance_result
- <
- point_type, point_type, bg::default_strategy
- >::type,
- ExpectedResultType
- >
- ));
- // Default strategy
- v1 = point;
- v2 = point;
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2),
- bg::comparable_distance(point, point),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, point),
- bg::comparable_distance(point, point),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(point, v2),
- bg::comparable_distance(point, point),
- 0.0001);
- v1 = point;
- v2 = seg;
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2),
- bg::comparable_distance(point, seg),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, seg),
- bg::comparable_distance(point, seg),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(point, v2),
- bg::comparable_distance(point, seg),
- 0.0001);
- v1 = point;
- v2 = linestring;
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2),
- bg::comparable_distance(point, linestring),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, linestring),
- bg::comparable_distance(point, linestring),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(point, v2),
- bg::comparable_distance(point, linestring),
- 0.0001);
- }
- };
- template <typename T, typename ExpectedResultType = T>
- struct test_variant_with_strategy
- {
- static inline void apply()
- {
- typedef bg::strategy::distance::projected_point<T> strategy_type;
- typedef bg::model::point<T, 2, bg::cs::cartesian> point_type;
- typedef bg::model::segment<point_type> segment_type;
- typedef bg::model::linestring<point_type> linestring_type;
- typedef bg::model::multi_linestring
- <
- linestring_type
- > multi_linestring_type;
- typedef boost::variant
- <
- segment_type, linestring_type, multi_linestring_type
- > variant_type;
- segment_type seg;
- bg::read_wkt("LINESTRING(1 1,4 4)", seg);
- linestring_type ls;
- bg::read_wkt("LINESTRING(-1 -1,-1 0,0 0,0 -1,-1 -1)", ls);
- multi_linestring_type mls;
- bg::read_wkt("MULTILINESTRING((10 0,20 0),(30 0,40 0))", mls);
- variant_type v1, v2;
- strategy_type strategy;
- BOOST_MPL_ASSERT((
- boost::is_same
- <
- typename bg::comparable_distance_result
- <
- variant_type, variant_type, strategy_type
- >::type,
- ExpectedResultType
- >
- ));
- BOOST_MPL_ASSERT((
- boost::is_same
- <
- typename bg::comparable_distance_result
- <
- segment_type, linestring_type, strategy_type
- >::type,
- ExpectedResultType
- >
- ));
- // Passed strategy
- v1 = seg;
- v2 = seg;
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2, strategy),
- bg::comparable_distance(seg, seg, strategy),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, seg, strategy),
- bg::comparable_distance(seg, seg, strategy),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(seg, v2, strategy),
- bg::comparable_distance(seg, seg, strategy),
- 0.0001);
- v1 = seg;
- v2 = ls;
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2, strategy),
- bg::comparable_distance(seg, ls, strategy),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, ls, strategy),
- bg::comparable_distance(seg, ls, strategy),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(seg, v2, strategy),
- bg::comparable_distance(seg, ls, strategy),
- 0.0001);
- v1 = seg;
- v2 = mls;
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2, strategy),
- bg::comparable_distance(seg, mls, strategy),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, mls, strategy),
- bg::comparable_distance(seg, mls, strategy),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(seg, v2, strategy),
- bg::comparable_distance(seg, mls, strategy),
- 0.0001);
- v1 = ls;
- v2 = mls;
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2, strategy),
- bg::comparable_distance(ls, mls, strategy),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(v1, mls, strategy),
- bg::comparable_distance(ls, mls, strategy),
- 0.0001);
- BOOST_CHECK_CLOSE(bg::comparable_distance(ls, v2, strategy),
- bg::comparable_distance(ls, mls, strategy),
- 0.0001);
- }
- };
- template <typename T, bool IsIntergral = boost::is_integral<T>::value>
- struct check_result
- {
- template <typename ExpectedResult>
- static inline void apply(T const& value,
- ExpectedResult const& expected_value)
- {
- BOOST_CHECK_EQUAL(value, expected_value);
- }
- };
- template <typename T>
- struct check_result<T, false>
- {
- template <typename ExpectedResult>
- static inline void apply(T const& value,
- ExpectedResult const& expected_value)
- {
- BOOST_CHECK_CLOSE(value, expected_value, 0.0001);
- }
- };
- template <typename T>
- struct test_variant_boxes
- {
- static inline void apply()
- {
- typedef bg::model::point<T, 2, bg::cs::cartesian> point_type;
- typedef bg::model::box<point_type> box_type;
- typedef boost::variant<box_type> variant_type;
- box_type box1, box2;
- bg::read_wkt("BOX(-1 -1,0 0)", box1);
- bg::read_wkt("BOX(1 1,2 2)", box2);
- variant_type v1 = box1, v2 = box2;
- typedef typename boost::mpl::if_c
- <
- boost::is_float<T>::value,
- double,
- typename bg::util::detail::default_integral::type
- >::type expected_result_type;
- BOOST_MPL_ASSERT((
- boost::is_same
- <
- typename bg::comparable_distance_result
- <
- variant_type, variant_type, bg::default_strategy
- >::type,
- expected_result_type
- >
- ));
- // Default strategy
- check_result<T>::apply(bg::comparable_distance(v1, v2),
- bg::comparable_distance(box1, box2));
- check_result<T>::apply(bg::comparable_distance(v1, box2),
- bg::comparable_distance(box1, box2));
- check_result<T>::apply(bg::comparable_distance(box1, v2),
- bg::comparable_distance(box1, box2));
- }
- };
- int test_main(int, char* [])
- {
- test_double_result_from_integer<int>();
- test_double_result_from_integer<boost::long_long_type>();
- test_all<bg::model::d2::point_xy<float> >();
- test_all<bg::model::d2::point_xy<double> >();
- #ifdef HAVE_TTMATH
- test_all<bg::model::d2::point_xy<ttmath_big> >();
- #endif
- // test variant support
- test_variant_different_default_strategy<double>::apply();
- test_variant_same_default_strategy<double>::apply();
- test_variant_same_default_strategy<int>::apply();
- test_variant_same_default_strategy<long>::apply();
- test_variant_with_strategy<double>::apply();
- test_variant_with_strategy<float>::apply();
- test_variant_with_strategy<long double>::apply();
- test_variant_with_strategy<int, double>::apply();
- #ifdef HAVE_TTMATH
- test_variant_with_strategy<ttmath_big>::apply();
- #endif
- test_variant_boxes<double>::apply();
- test_variant_boxes<int>::apply();
- test_variant_boxes<long>::apply();
- return 0;
- }
|