123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182 |
- // Boost.Geometry (aka GGL, Generic Geometry Library)
- // Unit 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)
- #include <iostream>
- #include <string>
- #define BOOST_GEOMETRY_NO_BOOST_TEST
- #include <test_overlay_p_q.hpp>
- #include <boost/program_options.hpp>
- #include <boost/timer.hpp>
- template <typename Polygon>
- inline void make_star(Polygon& polygon,
- int count, double factor1, double factor2, long double offset = 0)
- {
- typedef typename bg::point_type<Polygon>::type p;
- typedef typename bg::select_most_precise
- <
- typename bg::coordinate_type<Polygon>::type,
- long double
- >::type coordinate_type;
- // Create star
- coordinate_type cx = 25.0;
- coordinate_type cy = 25.0;
- coordinate_type dx = 50.0;
- coordinate_type dy = 50.0;
- coordinate_type half = 0.5;
- coordinate_type two = 2.0;
- coordinate_type a1 = coordinate_type(factor1) * half * dx;
- coordinate_type b1 = coordinate_type(factor1) * half * dy;
- coordinate_type a2 = coordinate_type(factor2) * half * dx;
- coordinate_type b2 = coordinate_type(factor2) * half * dy;
- coordinate_type pi = boost::math::constants::pi<long double>();
- coordinate_type delta = pi * two / coordinate_type(count - 1);
- coordinate_type angle = coordinate_type(offset) * delta;
- for (int i = 0; i < count - 1; i++, angle += delta)
- {
- bool even = i % 2 == 0;
- coordinate_type s = sin(angle);
- coordinate_type c = cos(angle);
- coordinate_type x = cx + (even ? a1 : a2) * s;
- coordinate_type y = cy + (even ? b1 : b2) * c;
- bg::exterior_ring(polygon).push_back(bg::make<p>(x, y));
- }
- bg::exterior_ring(polygon).push_back(bg::exterior_ring(polygon).front());
- }
- template <typename T, typename CalculationType>
- void test_star(int count, int min_points, int max_points, T rotation, p_q_settings const& settings)
- {
- boost::timer t;
- typedef bg::model::d2::point_xy<T> point_type;
- typedef bg::model::polygon<point_type> polygon;
- int n = 0;
- for (int c = 0; c < count; c++)
- {
- for (int i = min_points; i <= max_points; i++)
- {
- std::ostringstream out;
- out << "_" << string_from_type<T>::name() << "_"
- << string_from_type<CalculationType>::name() << "_"
- << i << "_int";
- polygon p;
- make_star(p, i * 2 + 1, 0.5, 1.0);
- polygon q;
- make_star(q, i * 2 + 1, 0.5, 1.0, rotation);
- if (! test_overlay_p_q
- <
- polygon,
- CalculationType
- >(out.str(), p, q, settings))
- {
- return;
- }
- n++;
- }
- }
- std::cout
- << "polygons: " << n
- << " type: " << string_from_type<T>::name()
- << " time: " << t.elapsed() << std::endl;
- }
- template <typename T, typename CalculationType>
- void test_type(int count, int min_points, int max_points, T rotation, p_q_settings const& settings)
- {
- test_star<T, CalculationType>(count, min_points, max_points, rotation, settings);
- }
- template <typename T>
- void test_all(std::string const& type, int count, int min_points, int max_points, T rotation, p_q_settings settings)
- {
- if (type == "float")
- {
- settings.tolerance = 1.0e-3;
- test_type<float, float>(count, min_points, max_points, rotation, settings);
- }
- else if (type == "double")
- {
- test_type<double, double>(count, min_points, max_points, rotation, settings);
- }
- #if defined(HAVE_TTMATH)
- else if (type == "ttmath")
- {
- test_type<ttmath_big, ttmath_big>(count, min_points, max_points, rotation, settings);
- }
- #endif
- }
- 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));
- std::string type = "float";
- int min_points = 9;
- int max_points = 9;
- bool ccw = false;
- bool open = false;
- double rotation = 1.0e-13;
- p_q_settings settings;
- 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")
- ("min_points", po::value<int>(&min_points)->default_value(9), "Minimum number of points")
- ("max_points", po::value<int>(&max_points)->default_value(9), "Maximum number of points")
- ("rotation", po::value<double>(&rotation)->default_value(1.0e-13), "Rotation angle")
- ("ccw", po::value<bool>(&ccw)->default_value(false), "Counter clockwise polygons")
- ("open", po::value<bool>(&open)->default_value(false), "Open polygons")
- ("type", po::value<std::string>(&type)->default_value("float"), "Type (float,double)")
- ("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"))
- {
- std::cout << description << std::endl;
- return 1;
- }
- test_all(type, count, min_points, max_points, rotation, settings);
- }
- catch(std::exception const& e)
- {
- std::cout << "Exception " << e.what() << std::endl;
- }
- catch(...)
- {
- std::cout << "Other exception" << std::endl;
- }
- return 0;
- }
|