123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218 |
- // Boost.Geometry (aka GGL, Generic Geometry Library)
- // Robustness Test
- // Copyright (c) 2009-2012 Barend Gehrels, 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)
- #define BOOST_GEOMETRY_REPORT_OVERLAY_ERROR
- #define BOOST_GEOMETRY_NO_BOOST_TEST
- #include <test_overlay_p_q.hpp>
- #include <boost/program_options.hpp>
- #include <boost/random/linear_congruential.hpp>
- #include <boost/random/uniform_int.hpp>
- #include <boost/random/uniform_real.hpp>
- #include <boost/random/variate_generator.hpp>
- #include <boost/timer.hpp>
- template <typename Polygon, typename Generator>
- inline void make_polygon(Polygon& polygon, Generator& generator, bool triangular)
- {
- typedef typename bg::point_type<Polygon>::type point_type;
- typedef typename bg::coordinate_type<Polygon>::type coordinate_type;
- coordinate_type x, y;
- x = generator();
- y = generator();
- typename bg::ring_type<Polygon>::type& ring = bg::exterior_ring(polygon);
- point_type p;
- bg::set<0>(p, x); bg::set<1>(p, y); ring.push_back(p);
- bg::set<0>(p, x); bg::set<1>(p, y + 1); ring.push_back(p);
- bg::set<0>(p, x + 1); bg::set<1>(p, y + 1); ring.push_back(p);
- bg::set<0>(p, x + 1); bg::set<1>(p, y); ring.push_back(p);
- bg::set<0>(p, x); bg::set<1>(p, y); ring.push_back(p);
- if (triangular)
- {
- // Remove a point depending on generator
- int c = generator() % 4;
- if (c >= 1 && c <= 3)
- {
- ring.erase(ring.begin() + c);
- }
- }
- }
- template <typename MultiPolygon, typename Generator>
- bool test_recursive_boxes(MultiPolygon& result, int& index,
- Generator& generator,
- int level, bool triangular, p_q_settings const& settings)
- {
- MultiPolygon p, q;
- // Generate two boxes
- if (level == 0)
- {
- p.resize(1);
- q.resize(1);
- make_polygon(p.front(), generator, triangular);
- make_polygon(q.front(), generator, triangular);
- bg::correct(p);
- bg::correct(q);
- }
- else
- {
- bg::correct(p);
- bg::correct(q);
- if (! test_recursive_boxes(p, index, generator, level - 1, triangular, settings)
- || ! test_recursive_boxes(q, index, generator, level - 1, triangular, settings))
- {
- return false;
- }
- }
- typedef typename boost::range_value<MultiPolygon>::type polygon;
- std::ostringstream out;
- out << "recursive_box_" << index++ << "_" << level;
- if (! test_overlay_p_q
- <
- polygon,
- typename bg::coordinate_type<MultiPolygon>::type
- >(out.str(), p, q, settings))
- {
- return false;
- }
- MultiPolygon mp;
- bg::detail::union_::union_insert
- <
- polygon
- >(p, q, std::back_inserter(mp));
- bg::unique(mp);
- bg::simplify(mp, result, 0.01);
- bg::correct(mp);
- return true;
- }
- template <typename T, bool Clockwise, bool Closed>
- void test_all(int seed, int count, int field_size, int level, bool triangular, p_q_settings const& settings)
- {
- boost::timer t;
- typedef boost::minstd_rand base_generator_type;
- base_generator_type generator(seed);
- boost::uniform_int<> random_coordinate(0, field_size - 1);
- boost::variate_generator<base_generator_type&, boost::uniform_int<> >
- coordinate_generator(generator, random_coordinate);
- typedef bg::model::polygon
- <
- bg::model::d2::point_xy<T>, Clockwise, Closed
- > polygon;
- typedef bg::model::multi_polygon<polygon> mp;
- int index = 0;
- for(int i = 0; i < count; i++)
- {
- mp p;
- test_recursive_boxes<mp>(p, index, coordinate_generator, level, triangular, settings);
- }
- std::cout
- << "polygons: " << index
- << " type: " << string_from_type<T>::name()
- << " time: " << t.elapsed() << std::endl;
- }
- int main(int argc, char** argv)
- {
- try
- {
- namespace po = boost::program_options;
- po::options_description description("=== recursive_polygons ===\nAllowed options");
- int count = 1;
- int seed = static_cast<unsigned int>(std::time(0));
- int level = 3;
- int field_size = 10;
- bool ccw = false;
- bool open = false;
- p_q_settings settings;
- std::string form = "box";
- description.add_options()
- ("help", "Help message")
- ("seed", po::value<int>(&seed), "Initialization seed for random generator")
- ("count", po::value<int>(&count)->default_value(1), "Number of tests")
- ("diff", po::value<bool>(&settings.also_difference)->default_value(false), "Include testing on difference")
- ("validity", po::value<bool>(&settings.validity)->default_value(true), "Include testing on validity")
- ("level", po::value<int>(&level)->default_value(3), "Level to reach (higher->slower)")
- ("size", po::value<int>(&field_size)->default_value(10), "Size of the field")
- ("form", po::value<std::string>(&form)->default_value("box"), "Form of the polygons (box, triangle)")
- ("ccw", po::value<bool>(&ccw)->default_value(false), "Counter clockwise polygons")
- ("open", po::value<bool>(&open)->default_value(false), "Open polygons")
- ("wkt", po::value<bool>(&settings.wkt)->default_value(false), "Create a WKT of the inputs, for all tests")
- ("svg", po::value<bool>(&settings.svg)->default_value(false), "Create a SVG for all tests")
- ;
- po::variables_map varmap;
- po::store(po::parse_command_line(argc, argv, description), varmap);
- po::notify(varmap);
- if (varmap.count("help")
- || (form != "box" && form != "triangle"))
- {
- std::cout << description << std::endl;
- return 1;
- }
- bool triangular = form != "box";
- if (ccw && open)
- {
- test_all<double, false, false>(seed, count, field_size, level, triangular, settings);
- }
- else if (ccw)
- {
- test_all<double, false, true>(seed, count, field_size, level, triangular, settings);
- }
- else if (open)
- {
- test_all<double, true, false>(seed, count, field_size, level, triangular, settings);
- }
- else
- {
- test_all<double, true, true>(seed, count, field_size, level, triangular, settings);
- }
- #if defined(HAVE_TTMATH)
- // test_all<ttmath_big, true, true>(seed, count, max, svg, level);
- #endif
- }
- catch(std::exception const& e)
- {
- std::cout << "Exception " << e.what() << std::endl;
- }
- catch(...)
- {
- std::cout << "Other exception" << std::endl;
- }
- return 0;
- }
|