#### /Src/Dependencies/Boost/libs/geometry/test/algorithms/overlay/robustness/star_comb.hpp

C++ Header | 130 lines | 93 code | 26 blank | 11 comment | 9 complexity | c5e551696445797c0142d96b1f12bcd5 MD5 | raw file
```  1// Boost.Geometry (aka GGL, Generic Geometry Library)
2
3// Copyright (c) 2009-2011 Barend Gehrels, Amsterdam, the Netherlands.
4
5// Use, modification and distribution is subject to the Boost Software License,
6// Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
8
9#ifndef BOOST_GEOMETRY_STAR_COMB_HPP
10#define BOOST_GEOMETRY_STAR_COMB_HPP
11
12#include <iostream>
13#include <string>
14#include <vector>
15
16#include <boost/timer.hpp>
17#include <boost/tuple/tuple.hpp>
18#include <boost/math/constants/constants.hpp>
19
20#include <test_overlay_p_q.hpp>
21
22
24inline void make_star(Polygon& polygon, AddFunctor functor,
25            int count, double factor1, double factor2,
26            double offset = 0.0,
27            bool close = true,
28            double orientation = 1.0)
29{
30    // Create star
31    double cx = 25.0;
32    double cy = 25.0;
33
34    double dx = 50.0;
35    double dy = 50.0;
36
37    double a1 = factor1 * 0.5 * dx;
38    double b1 = factor1 * 0.5 * dy;
39    double a2 = factor2 * 0.5 * dx;
40    double b2 = factor2 * 0.5 * dy;
41
42    double delta = orientation * boost::math::constants::pi<double>() * 2.0 / (count - 1);
43    double angle = offset * delta;
44    double x0, y0;
45    bool first = true;
46    for (int i = 0; i < count - 1; i++, angle += delta)
47    {
48        bool even = i % 2 == 0;
49        double x = cx + (even ? a1 : a2) * sin(angle);
50        double y = cy + (even ? b1 : b2) * cos(angle);
51        functor(polygon, x, y, i);
52        if (first)
53        {
54            x0 = x;
55            y0 = y;
56            first = false;
57        }
58
59    }
60    if (close)
61    {
62        functor(polygon, x0, y0, count);
63    }
64}
65
66template <typename Vector>
67void ccw_pushback(Vector& vector, double x, double y, int)
68{
69    vector.push_back(boost::make_tuple(x, y));
70}
71
73inline void make_comb(Polygon& polygon, AddFunctor functor,
74            int count,
75            bool close = true,
76            bool clockwise = true)
77{
78    int n = 0;
79
80    if (! clockwise)
81    {
82        typedef boost::tuple<double, double>  tup;
83        typedef std::vector<tup> vector_type;
84        vector_type vec;
85
86        // Create in clockwise order
87        make_comb(vec, ccw_pushback<vector_type>, count, close, true);
88
90        // (For GCC 3.4 have it const)
91        vector_type const& v = vec;
92        for (vector_type::const_reverse_iterator
93            it = v.rbegin(); it != v.rend(); ++it)
94        {
95            functor(polygon, it->get<0>(), it->get<1>(), n++);
96        }
97        return;
98    }
99
100    // Create comb
101    functor(polygon, 25.0, 0.0, n++);
102    functor(polygon, 0.0, 25.0, n++);
103    functor(polygon, 25.0, 50.0, n++);
104
105    // Function parameters
106    double diff = (25.0 / (count - 0.5)) / 2.0;
107
108    double b1 = -25.0;
109    double b2 = 25.0 - diff * 2.0;
110
111    double x1 = 50.0, x2 = 25.0;
112
113    for (int i = 0; i < count - 1; i++)
114    {
115        functor(polygon, x1, (x1 + b1), n++); x1 -= diff;
116        functor(polygon, x1, (x1 + b1), n++); x1 -= diff;
117        functor(polygon, x2, (x2 + b2), n++); x2 -= diff;
118        functor(polygon, x2, (x2 + b2), n++); x2 -= diff;
119    }
120    functor(polygon, x1, (x1 + b1), n++);
121
122    if (close)
123    {
124        functor(polygon, 25.0, 0.0, 4);
125    }
126}
127
128
129#endif // BOOST_GEOMETRY_STAR_COMB_HPP
130
```