123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863864865866867868869870871872873874875876877878879880881882883884885886887888889890891892893894895896897898899900901 |
- // Boost.Geometry (aka GGL, Generic Geometry Library)
- // Unit Test
- // Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands.
- // This file was modified by Oracle on 2017.
- // Modifications copyright (c) 2017, Oracle and/or its affiliates.
- // 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)
- #include <iostream>
- #include <geometry_test_common.hpp>
- #include <boost/array.hpp>
- #include <boost/foreach.hpp>
- #include <boost/geometry/algorithms/intersection.hpp>
- #include <boost/geometry/algorithms/make.hpp>
- #include <boost/geometry/algorithms/detail/overlay/get_turn_info.hpp>
- #include <boost/geometry/algorithms/detail/overlay/debug_turn_info.hpp>
- #include <boost/geometry/geometries/point_xy.hpp>
- #if defined(TEST_WITH_SVG)
- # include <boost/geometry/io/svg/svg_mapper.hpp>
- #endif
- // For test purposes, returns the point specified in the constructor
- template <typename Point>
- struct sub_range_from_points
- {
- typedef Point point_type;
- sub_range_from_points(Point const& i, Point const& j, Point const& k)
- {
- m_points[0] = i;
- m_points[1] = j;
- m_points[2] = k;
- }
- static inline bool is_first_segment() { return false; }
- static inline bool is_last_segment() { return false; }
- static inline std::size_t size() { return 3; }
- inline Point const& at(std::size_t index) const
- {
- return m_points[index % 3];
- }
- private :
- boost::array<Point, 3> m_points;
- };
- template <typename P, typename T>
- void test_with_point(std::string const& caseid,
- T pi_x, T pi_y, T pj_x, T pj_y, T pk_x, T pk_y,
- T qi_x, T qi_y, T qj_x, T qj_y, T qk_x, T qk_y,
- bg::detail::overlay::method_type expected_method,
- bool expected_touch_only,
- T ip_x, T ip_y,
- std::string const& expected,
- T ip_x2, T ip_y2)
- {
- P pi = bg::make<P>(pi_x, pi_y);
- P pj = bg::make<P>(pj_x, pj_y);
- P pk = bg::make<P>(pk_x, pk_y);
- P qi = bg::make<P>(qi_x, qi_y);
- P qj = bg::make<P>(qj_x, qj_y);
- P qk = bg::make<P>(qk_x, qk_y);
- typedef typename bg::strategy::intersection::services::default_strategy
- <
- typename bg::cs_tag<P>::type
- >::type strategy_type;
- typedef typename bg::detail::no_rescale_policy rescale_policy_type;
- typedef bg::detail::overlay::turn_info
- <
- P,
- typename bg::detail::segment_ratio_type<P, rescale_policy_type>::type
- > turn_info;
- typedef std::vector<turn_info> tp_vector;
- turn_info model;
- tp_vector info;
- strategy_type strategy;
- rescale_policy_type rescale_policy;
- sub_range_from_points<P> sub_range_p(pi, pj, pk);
- sub_range_from_points<P> sub_range_q(qi, qj, qk);
- bg::detail::overlay::get_turn_info
- <
- bg::detail::overlay::assign_null_policy
- >::apply(sub_range_p, sub_range_q, model, strategy, rescale_policy, std::back_inserter(info));
- if (info.size() == 0)
- {
- BOOST_CHECK_EQUAL(expected_method,
- bg::detail::overlay::method_none);
- }
- std::string detected;
- std::string method;
- for (typename tp_vector::const_iterator it = info.begin(); it != info.end(); ++it)
- {
- for (int t = 0; t < 2; t++)
- {
- detected += bg::operation_char(it->operations[t].operation);
- method += bg::method_char(it->method);
- }
- }
- BOOST_CHECK_MESSAGE(detected == expected,
- caseid
- << (caseid.find("_") == std::string::npos ? " " : "")
- << " method: " << method
- << " detected: " << detected
- << " expected: " << expected);
- if (! info.empty())
- {
- BOOST_CHECK_EQUAL(info[0].method, expected_method);
- BOOST_CHECK_MESSAGE(info[0].touch_only == expected_touch_only,
- caseid
- << " detected: " << info[0].touch_only
- << " expected: " << expected_touch_only);
- BOOST_CHECK_CLOSE(bg::get<0>(info[0].point), ip_x, 0.001);
- BOOST_CHECK_CLOSE(bg::get<1>(info[0].point), ip_y, 0.001);
- if (info.size() > 1)
- {
- BOOST_CHECK_EQUAL(info.size(), 2u);
- BOOST_CHECK_EQUAL(info[1].method, expected_method);
- BOOST_CHECK_CLOSE(bg::get<0>(info[1].point), ip_x2, 0.001);
- BOOST_CHECK_CLOSE(bg::get<1>(info[1].point), ip_y2, 0.001);
- }
- }
- #if defined(TEST_WITH_SVG)
- {
- std::ostringstream filename;
- filename << "get_turn_info_" << caseid
- << "_" << string_from_type<typename bg::coordinate_type<P>::type>::name()
- << ".svg";
- std::ofstream svg(filename.str().c_str());
- bg::svg_mapper<P> mapper(svg, 500, 500);
- mapper.add(bg::make<P>(0, 0));
- mapper.add(bg::make<P>(10, 10));
- bg::model::linestring<P> p; p.push_back(pi); p.push_back(pj); p.push_back(pk);
- bg::model::linestring<P> q; q.push_back(qi); q.push_back(qj); q.push_back(qk);
- mapper.map(p, "opacity:0.8;stroke:rgb(0,192,0);stroke-width:3");
- mapper.map(q, "opacity:0.8;stroke:rgb(0,0,255);stroke-width:3");
- std::string style = ";font-family='Verdana';font-weight:bold";
- std::string align = ";text-anchor:end;text-align:end";
- int offset = 8;
- mapper.text(pi, "pi", "fill:rgb(0,192,0)" + style, offset, offset);
- mapper.text(pj, "pj", "fill:rgb(0,192,0)" + style, offset, offset);
- mapper.text(pk, "pk", "fill:rgb(0,192,0)" + style, offset, offset);
- mapper.text(qi, "qi", "fill:rgb(0,0,255)" + style + align, -offset, offset);
- mapper.text(qj, "qj", "fill:rgb(0,0,255)" + style + align, -offset, offset);
- mapper.text(qk, "qk", "fill:rgb(0,0,255)" + style + align, -offset, offset);
- int factor = 1; // second info, if any, will go left by factor -1
- int ch = '1';
- for (typename tp_vector::const_iterator it = info.begin();
- it != info.end();
- ++it, factor *= -1, ch++)
- {
- bool at_j = it->method == bg::detail::overlay::method_crosses;
- std::string op;
- op += bg::operation_char(it->operations[0].operation);
- align = ";text-anchor:middle;text-align:center";
- mapper.text(at_j ? pj : pk, op, "fill:rgb(255,128,0)" + style + align, offset * factor, -offset);
- op.clear();
- op += bg::operation_char(it->operations[1].operation);
- mapper.text(at_j ? qj : qk, op, "fill:rgb(255,128,0)" + style + align, offset * factor, -offset);
- // Map intersection point + method
- mapper.map(it->point, "opacity:0.8;fill:rgb(255,0,0);stroke:rgb(0,0,100);stroke-width:1");
- op.clear();
- op += bg::method_char(it->method);
- op += ' ';
- op += (it->touch_only ? 'o' : '*');
- if (info.size() != 1)
- {
- op += ch;
- op += " p:"; op += bg::operation_char(it->operations[0].operation);
- op += " q:"; op += bg::operation_char(it->operations[1].operation);
- }
- mapper.text(it->point, op, "fill:rgb(255,0,0)" + style, offset, -offset);
- }
- }
- #endif
- }
- template <typename P, typename T>
- void test_both(std::string const& caseid,
- T pi_x, T pi_y, T pj_x, T pj_y, T pk_x, T pk_y,
- T qi_x, T qi_y, T qj_x, T qj_y, T qk_x, T qk_y,
- bg::detail::overlay::method_type method
- = bg::detail::overlay::method_none,
- bool expected_touch_only = false,
- T ip_x = -1, T ip_y = -1,
- std::string const& expected = "",
- T ip_x2 = -1, T ip_y2 = -1)
- {
- test_with_point<P, double>(caseid,
- pi_x, pi_y, pj_x, pj_y, pk_x, pk_y,
- qi_x, qi_y, qj_x, qj_y, qk_x, qk_y,
- method, expected_touch_only, ip_x, ip_y, expected, ip_x2, ip_y2);
- std::string reversed(expected.rbegin(), expected.rend());
-
- if (ip_x2 >= 0 && ip_y2 >= 0)
- {
- std::swap(ip_x, ip_x2);
- std::swap(ip_y, ip_y2);
- }
- test_with_point<P, double>(caseid + "_r",
- qi_x, qi_y, qj_x, qj_y, qk_x, qk_y, // q
- pi_x, pi_y, pj_x, pj_y, pk_x, pk_y, // p
- method, expected_touch_only, ip_x, ip_y, reversed, ip_x2, ip_y2);
- }
- template <typename P>
- void test_all()
- {
- using namespace bg::detail::overlay;
- // See powerpoint "doc/testcases/get_turn_info.ppt"
- // ------------------------------------------------------------------------
- // "Real" intersections ("i"), or, crossing
- // ------------------------------------------------------------------------
- test_both<P, double>("il1",
- 5, 1, 5, 6, 7, 8, // p
- 3, 3, 7, 5, 8, 3, // q
- method_crosses, false, 5, 4, "ui");
- test_both<P, double>("il2",
- 5, 1, 5, 6, 7, 8, // p
- 3, 5, 7, 5, 3, 3, // q
- method_crosses, false, 5, 5, "ui");
- test_both<P, double>("il3",
- 5, 1, 5, 6, 7, 8, // p
- 3, 3, 7, 5, 3, 5, // q
- method_crosses, false, 5, 4, "ui");
- test_both<P, double>("il4",
- 5, 1, 5, 6, 7, 8, // p
- 3, 3, 7, 5, 4, 8, // q
- method_crosses, false, 5, 4, "ui");
- test_both<P, double>("ir1",
- 5, 1, 5, 6, 7, 8, // p
- 7, 5, 3, 3, 2, 5, // q
- method_crosses, false, 5, 4, "iu");
- // ------------------------------------------------------------------------
- // TOUCH INTERIOR or touch in the middle ("m")
- // ------------------------------------------------------------------------
- test_both<P, double>("ml1",
- 5, 1, 5, 6, 7, 8, // p
- 3, 3, 5, 4, 7, 3, // q
- method_touch_interior, false, 5, 4, "ui");
- test_both<P, double>("ml2",
- 5, 1, 5, 6, 7, 8, // p
- 3, 3, 5, 4, 3, 6, // q
- method_touch_interior, true, 5, 4, "iu");
- test_both<P, double>("ml3",
- 5, 1, 5, 6, 7, 8, // p
- 3, 6, 5, 4, 3, 3, // q
- method_touch_interior, true, 5, 4, "uu");
- test_both<P, double>("mr1",
- 5, 1, 5, 6, 7, 8, // p
- 7, 3, 5, 4, 3, 3, // q
- method_touch_interior, false, 5, 4, "iu");
- test_both<P, double>("mr2",
- 5, 1, 5, 6, 7, 8, // p
- 7, 3, 5, 4, 7, 6, // q
- method_touch_interior, true, 5, 4, "ui");
- test_both<P, double>("mr3",
- 5, 1, 5, 6, 7, 8, // p
- 7, 6, 5, 4, 7, 3, // q
- method_touch_interior, true, 5, 4, "ii");
- test_both<P, double>("mcl",
- 5, 1, 5, 6, 7, 8, // p
- 3, 2, 5, 3, 5, 5, // q
- method_touch_interior, false, 5, 3, "cc");
- test_both<P, double>("mcr",
- 5, 1, 5, 6, 7, 8, // p
- 7, 2, 5, 3, 5, 5, // q
- method_touch_interior, false, 5, 3, "cc");
- test_both<P, double>("mclo",
- 5, 1, 5, 6, 7, 8, // p
- 3, 4, 5, 5, 5, 3, // q
- method_touch_interior, false, 5, 5, "ux");
- test_both<P, double>("mcro",
- 5, 1, 5, 6, 7, 8, // p
- 7, 4, 5, 5, 5, 3, // q
- method_touch_interior, false, 5, 5, "ix");
- // ------------------------------------------------------------------------
- // COLLINEAR
- // ------------------------------------------------------------------------
- test_both<P, double>("cll1",
- 5, 1, 5, 6, 3, 8, // p
- 5, 5, 5, 7, 3, 8, // q
- method_collinear, false, 5, 6, "ui");
- test_both<P, double>("cll2",
- 5, 1, 5, 6, 3, 8, // p
- 5, 3, 5, 5, 3, 6, // q
- method_collinear, false, 5, 5, "iu");
- test_both<P, double>("clr1",
- 5, 1, 5, 6, 3, 8, // p
- 5, 5, 5, 7, 6, 8, // q
- method_collinear, false, 5, 6, "ui");
- test_both<P, double>("clr2",
- 5, 1, 5, 6, 3, 8, // p
- 5, 3, 5, 5, 6, 6, // q
- method_collinear, false, 5, 5, "ui");
- test_both<P, double>("crl1",
- 5, 1, 5, 6, 7, 8, // p
- 5, 5, 5, 7, 3, 8, // q
- method_collinear, false, 5, 6, "iu");
- test_both<P, double>("crl2",
- 5, 1, 5, 6, 7, 8, // p
- 5, 3, 5, 5, 3, 6, // q
- method_collinear, false, 5, 5, "iu");
- test_both<P, double>("crr1",
- 5, 1, 5, 6, 7, 8, // p
- 5, 5, 5, 7, 6, 8, // q
- method_collinear, false, 5, 6, "iu");
- test_both<P, double>("crr2",
- 5, 1, 5, 6, 7, 8, // p
- 5, 3, 5, 5, 6, 6, // q
- method_collinear, false, 5, 5, "ui");
- // The next two cases are changed (BSG 2013-09-24), they contain turn info (#buffer_rt_g)
- // In new approach they are changed back (BSG 2013-10-20)
- test_both<P, double>("ccx1",
- 5, 1, 5, 6, 5, 8, // p
- 5, 5, 5, 7, 3, 8, // q
- method_collinear, false, 5, 6, "cc"); // "iu");
- test_both<P, double>("cxc1",
- 5, 1, 5, 6, 7, 8, // p
- 5, 3, 5, 5, 5, 7, // q
- method_collinear, false, 5, 5, "cc"); // "iu");
- // Bug in case #54 of "overlay_cases.hpp"
- test_both<P, double>("c_bug1",
- 5, 0, 2, 0, 2, 2, // p
- 4, 0, 1, 0, 1, 2, // q
- method_collinear, false, 2, 0, "iu");
- // ------------------------------------------------------------------------
- // COLLINEAR OPPOSITE
- // ------------------------------------------------------------------------
- test_both<P, double>("clo1",
- 5, 2, 5, 6, 3, 8, // p
- 5, 7, 5, 5, 3, 3, // q
- method_collinear, false, 5, 6, "ixxu", 5, 5);
- test_both<P, double>("clo2",
- 5, 2, 5, 6, 3, 8, // p
- 5, 7, 5, 5, 5, 2, // q
- method_collinear, false, 5, 6, "ix");
- // actually "xxix", xx is skipped everywhere
- test_both<P, double>("clo3",
- 5, 2, 5, 6, 3, 8, // p
- 5, 7, 5, 5, 7, 3, // q
- method_collinear, false, 5, 6, "ixxi", 5, 5);
- test_both<P, double>("cco1",
- 5, 2, 5, 6, 5, 8, // p
- 5, 7, 5, 5, 3, 3, // q
- method_collinear, false, 5, 5, "xu"); // "xuxx"
- test_both<P, double>("cco2",
- 5, 2, 5, 6, 5, 8, // p
- 5, 7, 5, 5, 5, 2); // q "xxxx"
- test_both<P, double>("cco3",
- 5, 2, 5, 6, 5, 8, // p
- 5, 7, 5, 5, 7, 3, // q
- method_collinear, false, 5, 5, "xi"); // "xixx"
- test_both<P, double>("cro1",
- 5, 2, 5, 6, 7, 8, // p
- 5, 7, 5, 5, 3, 3, // q
- method_collinear, false, 5, 6, "uxxu", 5, 5);
- test_both<P, double>("cro2",
- 5, 2, 5, 6, 7, 8, // p
- 5, 7, 5, 5, 5, 2, // q
- method_collinear, false, 5, 6, "ux"); // "xxux"
- test_both<P, double>("cro3",
- 5, 2, 5, 6, 7, 8, // p
- 5, 7, 5, 5, 7, 3, // q
- method_collinear, false, 5, 6, "uxxi", 5, 5);
- test_both<P, double>("cxo1",
- 5, 2, 5, 6, 3, 8, // p
- 5, 5, 5, 3, 3, 1, // q
- method_collinear, false, 5, 3, "xu");
- test_both<P, double>("cxo2",
- 5, 2, 5, 6, 3, 8, // p
- 5, 5, 5, 3, 5, 0); // q "xx"
- test_both<P, double>("cxo3",
- 5, 2, 5, 6, 3, 8, // p
- 5, 5, 5, 3, 7, 1, // q
- method_collinear, false, 5, 3, "xi");
- test_both<P, double>("cxo4",
- 5, 2, 5, 6, 3, 8, // p
- 5, 7, 5, 1, 3, 0, // q
- method_collinear, false, 5, 6, "ix");
- test_both<P, double>("cxo5",
- 5, 2, 5, 6, 5, 8, // p
- 5, 7, 5, 1, 3, 0); // q "xx"
- test_both<P, double>("cxo6",
- 5, 2, 5, 6, 7, 8, // p
- 5, 7, 5, 1, 3, 0, // q
- method_collinear, false, 5, 6, "ux");
- // Verify
- test_both<P, double>("cvo1",
- 5, 3, 5, 7, 7, 9, // p
- 5, 5, 5, 3, 3, 1 // q
- );
- test_both<P, double>("cvo2",
- 5, 3, 5, 7, 7, 9, // p
- 5, 4, 5, 2, 3, 0 // q
- );
- // ------------------------------------------------------------------------
- // TOUCH - both same
- // ------------------------------------------------------------------------
- // Both left, Q turns right
- test_both<P, double>("blr1",
- 5, 1, 5, 6, 4, 4, // p
- 3, 7, 5, 6, 3, 5, // q
- method_touch, true, 5, 6, "ui");
- test_both<P, double>("blr2",
- 5, 1, 5, 6, 1, 4, // p
- 3, 7, 5, 6, 3, 5, // q
- method_touch, false, 5, 6, "cc");
- test_both<P, double>("blr3",
- 5, 1, 5, 6, 3, 6, // p
- 3, 7, 5, 6, 3, 5, // q
- method_touch, false, 5, 6, "iu");
- test_both<P, double>("blr4",
- 5, 1, 5, 6, 1, 8, // p
- 3, 7, 5, 6, 3, 5, // q
- method_touch, false, 5, 6, "xu");
- test_both<P, double>("blr5",
- 5, 1, 5, 6, 4, 8, // p
- 3, 7, 5, 6, 3, 5, // q
- method_touch, true, 5, 6, "uu");
- test_both<P, double>("blr6",
- 5, 1, 5, 6, 6, 4, // p
- 3, 7, 5, 6, 3, 5, // q
- method_touch, true, 5, 6, "uu");
- test_both<P, double>("blr7",
- 5, 1, 5, 6, 3, 6, // p
- 3, 7, 5, 6, 5, 3, // q
- method_touch, false, 5, 6, "ix");
- test_both<P, double>("blr8",
- 5, 1, 5, 6, 3, 6, // p
- 3, 6, 5, 6, 5, 3, // q
- method_touch, false, 5, 6, "xx");
- test_both<P, double>("blr9",
- 5, 1, 5, 6, 3, 6, // p
- 3, 5, 5, 6, 5, 3, // q
- method_touch, false, 5, 6, "ux");
- // Variants
- test_both<P, double>("blr7-a",
- 5, 1, 5, 6, 3, 6, // p
- 5, 8, 5, 6, 5, 3, // q
- method_touch, false, 5, 6, "ix");
- test_both<P, double>("blr7-b", // in fact NOT "both-left"
- 5, 1, 5, 6, 3, 6, // p
- 6, 8, 5, 6, 5, 3, // q
- method_touch, false, 5, 6, "ix");
- // To check if "collinear-check" on other side
- // does not apply to this side
- test_both<P, double>("blr6-c1",
- 5, 1, 5, 6, 7, 5, // p
- 3, 7, 5, 6, 3, 5, // q
- method_touch, true, 5, 6, "uu");
- test_both<P, double>("blr6-c2",
- 5, 1, 5, 6, 7, 7, // p
- 3, 7, 5, 6, 3, 5, // q
- method_touch, true, 5, 6, "uu");
- // Both right, Q turns right
- test_both<P, double>("brr1",
- 5, 1, 5, 6, 6, 4, // p
- 7, 5, 5, 6, 7, 7, // q
- method_touch, true, 5, 6, "uu");
- test_both<P, double>("brr2",
- 5, 1, 5, 6, 9, 4, // p
- 7, 5, 5, 6, 7, 7, // q
- method_touch, false, 5, 6, "xu");
- test_both<P, double>("brr3",
- 5, 1, 5, 6, 7, 6, // p
- 7, 5, 5, 6, 7, 7, // q
- method_touch, false, 5, 6, "iu");
- test_both<P, double>("brr4",
- 5, 1, 5, 6, 9, 8, // p
- 7, 5, 5, 6, 7, 7, // q
- method_touch, false, 5, 6, "cc");
- test_both<P, double>("brr5",
- 5, 1, 5, 6, 6, 8, // p
- 7, 5, 5, 6, 7, 7, // q
- method_touch, true, 5, 6, "ui");
- test_both<P, double>("brr6",
- 5, 1, 5, 6, 4, 4, // p
- 7, 5, 5, 6, 7, 7, // q
- method_touch, true, 5, 6, "ui");
- // Both right, Q turns left
- test_both<P, double>("brl1",
- 5, 1, 5, 6, 6, 4, // p
- 7, 7, 5, 6, 7, 5, // q
- method_touch, true, 5, 6, "iu");
- test_both<P, double>("brl2",
- 5, 1, 5, 6, 9, 4, // p
- 7, 7, 5, 6, 7, 5, // q
- method_touch, false, 5, 6, "cc");
- test_both<P, double>("brl3",
- 5, 1, 5, 6, 7, 6, // p
- 7, 7, 5, 6, 7, 5, // q
- method_touch, false, 5, 6, "ui");
- test_both<P, double>("brl4",
- 5, 1, 5, 6, 9, 8, // p
- 7, 7, 5, 6, 7, 5, // q
- method_touch, false, 5, 6, "xi");
- test_both<P, double>("brl5",
- 5, 1, 5, 6, 6, 8, // p
- 7, 7, 5, 6, 7, 5, // q
- method_touch, true, 5, 6, "ii");
- test_both<P, double>("brl6",
- 5, 1, 5, 6, 4, 4, // p
- 7, 7, 5, 6, 7, 5, // q
- method_touch, true, 5, 6, "ii");
- test_both<P, double>("brl7",
- 5, 1, 5, 6, 7, 6, // p
- 7, 7, 5, 6, 5, 3, // q
- method_touch, false, 5, 6, "ux");
- test_both<P, double>("brl8",
- 5, 1, 5, 6, 7, 6, // p
- 7, 6, 5, 6, 5, 3, // q
- method_touch, false, 5, 6, "xx");
- test_both<P, double>("brl9",
- 5, 1, 5, 6, 7, 6, // p
- 7, 5, 5, 6, 5, 3, // q
- method_touch, false, 5, 6, "ix");
- // Variants
- test_both<P, double>("brl7-a",
- 5, 1, 5, 6, 7, 6, // p
- 5, 8, 5, 6, 5, 3, // q
- method_touch, false, 5, 6, "ux");
- test_both<P, double>("brl7-b", // in fact NOT "both right"
- 5, 1, 5, 6, 7, 6, // p
- 4, 8, 5, 6, 5, 3, // q
- method_touch, false, 5, 6, "ux");
- // Both left, Q turns left
- test_both<P, double>("bll1",
- 5, 1, 5, 6, 4, 4, // p
- 3, 5, 5, 6, 3, 7, // q
- method_touch, true, 5, 6, "ii");
- test_both<P, double>("bll2",
- 5, 1, 5, 6, 1, 4, // p
- 3, 5, 5, 6, 3, 7, // q
- method_touch, false, 5, 6, "xi");
- test_both<P, double>("bll3",
- 5, 1, 5, 6, 3, 6, // p
- 3, 5, 5, 6, 3, 7, // q
- method_touch, false, 5, 6, "ui");
- test_both<P, double>("bll4",
- 5, 1, 5, 6, 1, 8, // p
- 3, 5, 5, 6, 3, 7, // q
- method_touch, false, 5, 6, "cc");
- test_both<P, double>("bll5",
- 5, 1, 5, 6, 4, 8, // p
- 3, 5, 5, 6, 3, 7, // q
- method_touch, true, 5, 6, "iu");
- test_both<P, double>("bll6",
- 5, 1, 5, 6, 6, 4, // p
- 3, 5, 5, 6, 3, 7, // q
- method_touch, true, 5, 6, "iu");
- // TOUCH - COLLINEAR + one side
- // Collinear/left, Q turns right
- test_both<P, double>("t-clr1",
- 5, 1, 5, 6, 4, 4, // p
- 5, 8, 5, 6, 3, 5, // q
- method_touch, true, 5, 6, "ui");
- test_both<P, double>("t-clr2",
- 5, 1, 5, 6, 1, 4, // p
- 5, 8, 5, 6, 3, 5, // q
- method_touch, false, 5, 6, "cc");
- test_both<P, double>("t-clr3",
- 5, 1, 5, 6, 3, 6, // p
- 5, 8, 5, 6, 3, 5, // q
- method_touch, false, 5, 6, "iu");
- test_both<P, double>("t-clr4",
- 5, 1, 5, 6, 5, 8, // p
- 5, 8, 5, 6, 3, 5, // q
- method_touch, false, 5, 6, "xu");
- // 5 n.a.
- test_both<P, double>("t-clr6",
- 5, 1, 5, 6, 6, 4, // p
- 5, 8, 5, 6, 3, 5, // q
- method_touch, true, 5, 6, "uu");
- // Collinear/right, Q turns right
- test_both<P, double>("t-crr1",
- 5, 1, 5, 6, 6, 4, // p
- 7, 5, 5, 6, 5, 8, // q
- method_touch, true, 5, 6, "uu");
- test_both<P, double>("t-crr2",
- 5, 1, 5, 6, 9, 4, // p
- 7, 5, 5, 6, 5, 8, // q
- method_touch, false, 5, 6, "xu");
- test_both<P, double>("t-crr3",
- 5, 1, 5, 6, 7, 6, // p
- 7, 5, 5, 6, 5, 8, // q
- method_touch, false, 5, 6, "iu");
- test_both<P, double>("t-crr4",
- 5, 1, 5, 6, 5, 9, // p
- 7, 5, 5, 6, 5, 8, // q
- method_touch, false, 5, 6, "cc");
- // 5 n.a.
- test_both<P, double>("t-crr6",
- 5, 1, 5, 6, 4, 4, // p
- 7, 5, 5, 6, 5, 8, // q
- method_touch, true, 5, 6, "ui");
- // Collinear/right, Q turns left
- test_both<P, double>("t-crl1",
- 5, 1, 5, 6, 6, 4, // p
- 5, 7, 5, 6, 7, 5, // q
- method_touch, true, 5, 6, "iu");
- test_both<P, double>("t-crl2",
- 5, 1, 5, 6, 9, 4, // p
- 5, 7, 5, 6, 7, 5, // q
- method_touch, false, 5, 6, "cc");
- test_both<P, double>("t-crl3",
- 5, 1, 5, 6, 7, 6, // p
- 5, 7, 5, 6, 7, 5, // q
- method_touch, false, 5, 6, "ui");
- test_both<P, double>("t-crl4",
- 5, 1, 5, 6, 5, 8, // p
- 5, 7, 5, 6, 7, 5, // q
- method_touch, false, 5, 6, "xi");
- // 5 n.a.
- test_both<P, double>("t-crl6",
- 5, 1, 5, 6, 4, 4, // p
- 5, 7, 5, 6, 7, 5, // q
- method_touch, true, 5, 6, "ii");
- // Collinear/left, Q turns left
- test_both<P, double>("t-cll1",
- 5, 1, 5, 6, 4, 4, // p
- 3, 5, 5, 6, 5, 8, // q
- method_touch, true, 5, 6, "ii");
- test_both<P, double>("t-cll2",
- 5, 1, 5, 6, 1, 4, // p
- 3, 5, 5, 6, 5, 8, // q
- method_touch, false, 5, 6, "xi");
- test_both<P, double>("t-cll3",
- 5, 1, 5, 6, 3, 6, // p
- 3, 5, 5, 6, 5, 8, // q
- method_touch, false, 5, 6, "ui");
- test_both<P, double>("t-cll4",
- 5, 1, 5, 6, 5, 9, // p
- 3, 5, 5, 6, 5, 8, // q
- method_touch, false, 5, 6, "cc");
- // 5 n.a.
- test_both<P, double>("t-cll6",
- 5, 1, 5, 6, 6, 4, // p
- 3, 5, 5, 6, 5, 8, // q
- method_touch, true, 5, 6, "iu");
- // Left to right
- test_both<P, double>("lr1",
- 5, 1, 5, 6, 3, 3, // p
- 1, 5, 5, 6, 9, 5, // q
- method_touch, true, 5, 6, "ii");
- test_both<P, double>("lr2",
- 5, 1, 5, 6, 1, 5, // p
- 1, 5, 5, 6, 9, 5, // q
- method_touch, false, 5, 6, "xi");
- test_both<P, double>("lr3",
- 5, 1, 5, 6, 4, 8, // p
- 1, 5, 5, 6, 9, 5, // q
- method_touch, false, 5, 6, "ui");
- test_both<P, double>("lr4",
- 5, 1, 5, 6, 9, 5, // p
- 1, 5, 5, 6, 9, 5, // q
- method_touch, false, 5, 6, "cc");
- test_both<P, double>("lr5",
- 5, 1, 5, 6, 7, 3, // p
- 1, 5, 5, 6, 9, 5, // q
- method_touch, true, 5, 6, "iu");
- // otherwise case more thoroughly
- test_both<P, double>("lr3a",
- 5, 1, 5, 6, 1, 6, // p
- 1, 5, 5, 6, 9, 5, // q
- method_touch, false, 5, 6, "ui");
- test_both<P, double>("lr3b",
- 5, 1, 5, 6, 5, 10, // p
- 1, 5, 5, 6, 9, 5, // q
- method_touch, false, 5, 6, "ui");
- test_both<P, double>("lr3c",
- 5, 1, 5, 6, 8, 9, // p
- 1, 5, 5, 6, 9, 5, // q
- method_touch, false, 5, 6, "ui");
- test_both<P, double>("lr3d",
- 5, 1, 5, 6, 9, 7, // p
- 1, 5, 5, 6, 9, 5, // q
- method_touch, false, 5, 6, "ui");
- test_both<P, double>("lr3e",
- 5, 1, 5, 6, 9, 6, // p
- 1, 5, 5, 6, 9, 5, // q
- method_touch, false, 5, 6, "ui");
- // Right to left
- test_both<P, double>("rl1",
- 5, 1, 5, 6, 3, 3, // p
- 9, 5, 5, 6, 1, 5, // q
- method_touch, true, 5, 6, "ui");
- test_both<P, double>("rl2",
- 5, 1, 5, 6, 1, 5, // p
- 9, 5, 5, 6, 1, 5, // q
- method_touch, false, 5, 6, "cc");
- test_both<P, double>("rl3",
- 5, 1, 5, 6, 4, 8, // p
- 9, 5, 5, 6, 1, 5, // q
- method_touch, false, 5, 6, "iu");
- test_both<P, double>("rl4",
- 5, 1, 5, 6, 9, 5, // p
- 9, 5, 5, 6, 1, 5, // q
- method_touch, false, 5, 6, "xu");
- test_both<P, double>("rl5",
- 5, 1, 5, 6, 7, 3, // p
- 9, 5, 5, 6, 1, 5, // q
- method_touch, true, 5, 6, "uu");
- // Equal (p1/q1 are equal)
- test_both<P, double>("ebl1",
- 5, 1, 5, 6, 3, 4, // p
- 5, 1, 5, 6, 3, 8, // q
- method_equal, false, 5, 6, "ui");
- test_both<P, double>("ebl2",
- 5, 1, 5, 6, 3, 8, // p
- 5, 1, 5, 6, 3, 4, // q
- method_equal, false, 5, 6, "iu");
- test_both<P, double>("ebl3",
- 5, 1, 5, 6, 3, 8, // p
- 5, 1, 5, 6, 3, 8, // q
- method_equal, false, 5, 6, "cc");
- test_both<P, double>("ebl3-c1",
- 5, 1, 5, 6, 10, 1, // p
- 5, 1, 5, 6, 3, 8, // q
- method_equal, false, 5, 6, "iu");
- test_both<P, double>("ebr1",
- 5, 1, 5, 6, 7, 4, // p
- 5, 1, 5, 6, 7, 8, // q
- method_equal, false, 5, 6, "iu");
- test_both<P, double>("ebr2",
- 5, 1, 5, 6, 7, 8, // p
- 5, 1, 5, 6, 7, 4, // q
- method_equal, false, 5, 6, "ui");
- test_both<P, double>("ebr3",
- 5, 1, 5, 6, 7, 8, // p
- 5, 1, 5, 6, 7, 8, // q
- method_equal, false, 5, 6, "cc");
- test_both<P, double>("ebr3-c1",
- 5, 1, 5, 6, 0, 1, // p
- 5, 1, 5, 6, 7, 8, // q
- method_equal, false, 5, 6, "ui");
- test_both<P, double>("elr1",
- 5, 1, 5, 6, 7, 8, // p
- 5, 1, 5, 6, 3, 8, // q
- method_equal, false, 5, 6, "iu");
- test_both<P, double>("elr2",
- 5, 1, 5, 6, 3, 8, // p
- 5, 1, 5, 6, 7, 8, // q
- method_equal, false, 5, 6, "ui");
- test_both<P, double>("ec1",
- 5, 1, 5, 6, 5, 8, // p
- 5, 1, 5, 6, 5, 8, // q
- method_equal, false, 5, 6, "cc");
- test_both<P, double>("ec2",
- 5, 1, 5, 6, 5, 8, // p
- 5, 1, 5, 6, 5, 7, // q
- method_equal, false, 5, 6, "cc");
- test_both<P, double>("snl-1",
- 0, 3, 2, 3, 4, 3, // p
- 4, 3, 2, 3, 0, 3, // q
- method_touch, false, 2, 3, "xx");
- // BSG 2012-05-26 to be decided what's the problem here and what it tests...
- // Anyway, test results are not filled out.
- //test_both<P, double>("issue_buffer_mill",
- // 5.1983614873206241 , 6.7259025813913107 , 5.0499999999999998 , 6.4291796067500622 , 5.1983614873206241 , 6.7259025813913107, // p
- // 5.0499999999999998 , 6.4291796067500622 , 5.0499999999999998 , 6.4291796067500622 , 5.1983614873206241 , 6.7259025813913107, // q
- // method_collinear, false, 2, 0, "tt");
- }
- /***
- #include <boost/geometry/geometries/adapted/c_array.hpp>
- BOOST_GEOMETRY_REGISTER_C_ARRAY_CS(cs::cartesian)
- template <typename G>
- void test2(G const& geometry)
- {
- typedef typename bg::point_type<G>::type P;
- typedef typename bg::tag<G>::type T;
- typedef typename bg::tag<P>::type PT;
- std::cout << typeid(G).name() << std::endl;
- std::cout << typeid(T).name() << std::endl;
- std::cout << typeid(P).name() << std::endl;
- std::cout << typeid(PT).name() << std::endl;
- std::cout << bg::length(geometry) << std::endl;
- typedef bg::model::point<float, 3, bg::cs::cartesian> P2;
- bg::model::linestring<P2> out;
- bg::strategy::transform::scale_transformer<float[3], P2> scaler(5);
- bg::transform(geometry, out, scaler);
- std::cout << bg::dsv(out) << std::endl;
- }
- void test_f3()
- {
- float vertices[][3] = {
- {-1, -1, 1}, {1, -1, 1}, {1, 1, 1}, {-1, 1, 1},
- {-1, -1, -1}, {1, -1, -1}, {1, 1, -1}, {-1, 1, -1}
- };
- test2(vertices);
- }
- ***/
- int test_main(int, char* [])
- {
- test_all<bg::model::d2::point_xy<double> >();
- return 0;
- }
|