Chris@16: /*
Chris@16:   Copyright 2008 Intel Corporation
Chris@16: 
Chris@16:   Use, modification and distribution are subject to the Boost Software License,
Chris@16:   Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
Chris@16:   http://www.boost.org/LICENSE_1_0.txt).
Chris@16: */
Chris@16: #ifndef BOOST_POLYGON_POLYGON_SET_CONCEPT_HPP
Chris@16: #define BOOST_POLYGON_POLYGON_SET_CONCEPT_HPP
Chris@16: #include "polygon_set_data.hpp"
Chris@16: #include "detail/polygon_simplify.hpp"
Chris@16: namespace boost { namespace polygon{
Chris@16: 
Chris@16:   template <typename T, typename T2>
Chris@16:   struct is_either_polygon_set_type {
Chris@16:     typedef typename gtl_or<typename is_polygon_set_type<T>::type, typename is_polygon_set_type<T2>::type >::type type;
Chris@16:   };
Chris@16: 
Chris@16:   template <typename T>
Chris@16:   struct is_any_polygon_set_type {
Chris@16:     typedef typename gtl_or<typename is_polygon_45_or_90_set_type<T>::type, typename is_polygon_set_type<T>::type >::type type;
Chris@16:   };
Chris@16: 
Chris@16:   template <typename polygon_set_type>
Chris@16:   typename enable_if< typename is_any_polygon_set_type<polygon_set_type>::type,
Chris@16:                        typename polygon_set_traits<polygon_set_type>::iterator_type>::type
Chris@16:   begin_polygon_set_data(const polygon_set_type& polygon_set) {
Chris@16:     return polygon_set_traits<polygon_set_type>::begin(polygon_set);
Chris@16:   }
Chris@16: 
Chris@16:   template <typename polygon_set_type>
Chris@16:   typename enable_if< typename is_any_polygon_set_type<polygon_set_type>::type,
Chris@16:                        typename polygon_set_traits<polygon_set_type>::iterator_type>::type
Chris@16:   end_polygon_set_data(const polygon_set_type& polygon_set) {
Chris@16:     return polygon_set_traits<polygon_set_type>::end(polygon_set);
Chris@16:   }
Chris@16: 
Chris@16:   template <typename polygon_set_type>
Chris@16:   typename enable_if< typename is_polygon_set_type<polygon_set_type>::type,
Chris@16:                        bool>::type
Chris@16:   clean(const polygon_set_type& polygon_set) {
Chris@16:     return polygon_set_traits<polygon_set_type>::clean(polygon_set);
Chris@16:   }
Chris@16: 
Chris@16:   //assign
Chris@16:   template <typename polygon_set_type_1, typename polygon_set_type_2>
Chris@16:   typename enable_if< typename gtl_and<
Chris@16:     typename is_mutable_polygon_set_type<polygon_set_type_1>::type,
Chris@16:     typename is_any_polygon_set_type<polygon_set_type_2>::type>::type,
Chris@16:                        polygon_set_type_1>::type &
Chris@16:   assign(polygon_set_type_1& lvalue, const polygon_set_type_2& rvalue) {
Chris@16:     if(clean(rvalue))
Chris@16:       polygon_set_mutable_traits<polygon_set_type_1>::set(lvalue, begin_polygon_set_data(rvalue), end_polygon_set_data(rvalue));
Chris@16:     else {
Chris@16:       polygon_set_data<typename polygon_set_traits<polygon_set_type_2>::coordinate_type> ps;
Chris@16:       ps.insert(begin_polygon_set_data(rvalue), end_polygon_set_data(rvalue));
Chris@16:       ps.clean();
Chris@16:       polygon_set_mutable_traits<polygon_set_type_1>::set(lvalue, ps.begin(), ps.end());
Chris@16:     }
Chris@16:     return lvalue;
Chris@16:   }
Chris@16: 
Chris@16:   //get trapezoids
Chris@16:   template <typename output_container_type, typename polygon_set_type>
Chris@16:   typename enable_if< typename is_mutable_polygon_set_type<polygon_set_type>::type,
Chris@16:                       void>::type
Chris@16:   get_trapezoids(output_container_type& output, const polygon_set_type& polygon_set) {
Chris@16:     polygon_set_data<typename polygon_set_traits<polygon_set_type>::coordinate_type> ps;
Chris@16:     assign(ps, polygon_set);
Chris@16:     ps.get_trapezoids(output);
Chris@16:   }
Chris@16: 
Chris@16:   //get trapezoids
Chris@16:   template <typename output_container_type, typename polygon_set_type>
Chris@16:   typename enable_if< typename is_mutable_polygon_set_type<polygon_set_type>::type,
Chris@16:                       void>::type
Chris@16:   get_trapezoids(output_container_type& output, const polygon_set_type& polygon_set,
Chris@16:                  orientation_2d orient) {
Chris@16:     polygon_set_data<typename polygon_set_traits<polygon_set_type>::coordinate_type> ps;
Chris@16:     assign(ps, polygon_set);
Chris@16:     ps.get_trapezoids(output, orient);
Chris@16:   }
Chris@16: 
Chris@16:   //equivalence
Chris@16:   template <typename polygon_set_type_1, typename polygon_set_type_2>
Chris@16:   typename enable_if< typename gtl_and_3 <
Chris@16:     typename is_any_polygon_set_type<polygon_set_type_1>::type,
Chris@16:     typename is_any_polygon_set_type<polygon_set_type_2>::type,
Chris@16:     typename is_either_polygon_set_type<polygon_set_type_1, polygon_set_type_2>::type>::type,
Chris@16:                        bool>::type
Chris@16:   equivalence(const polygon_set_type_1& lvalue,
Chris@16:               const polygon_set_type_2& rvalue) {
Chris@16:     polygon_set_data<typename polygon_set_traits<polygon_set_type_1>::coordinate_type> ps1;
Chris@16:     assign(ps1, lvalue);
Chris@16:     polygon_set_data<typename polygon_set_traits<polygon_set_type_2>::coordinate_type> ps2;
Chris@16:     assign(ps2, rvalue);
Chris@16:     return ps1 == ps2;
Chris@16:   }
Chris@16: 
Chris@16:   //clear
Chris@16:   template <typename polygon_set_type>
Chris@16:   typename enable_if< typename is_mutable_polygon_set_type<polygon_set_type>::type,
Chris@16:                        void>::type
Chris@16:   clear(polygon_set_type& polygon_set) {
Chris@16:     polygon_set_data<typename polygon_set_traits<polygon_set_type>::coordinate_type> ps;
Chris@16:     assign(polygon_set, ps);
Chris@16:   }
Chris@16: 
Chris@16:   //empty
Chris@16:   template <typename polygon_set_type>
Chris@16:   typename enable_if< typename is_mutable_polygon_set_type<polygon_set_type>::type,
Chris@16:                        bool>::type
Chris@16:   empty(const polygon_set_type& polygon_set) {
Chris@16:     if(clean(polygon_set)) return begin_polygon_set_data(polygon_set) == end_polygon_set_data(polygon_set);
Chris@16:     polygon_set_data<typename polygon_set_traits<polygon_set_type>::coordinate_type> ps;
Chris@16:     assign(ps, polygon_set);
Chris@16:     ps.clean();
Chris@16:     return ps.empty();
Chris@16:   }
Chris@16: 
Chris@16:   //extents
Chris@16:   template <typename polygon_set_type, typename rectangle_type>
Chris@16:   typename enable_if< typename gtl_and<
Chris@16:     typename is_mutable_polygon_set_type<polygon_set_type>::type,
Chris@16:     typename is_mutable_rectangle_concept<typename geometry_concept<rectangle_type>::type>::type>::type,
Chris@16:                        bool>::type
Chris@16:   extents(rectangle_type& extents_rectangle,
Chris@16:           const polygon_set_type& polygon_set) {
Chris@16:     clean(polygon_set);
Chris@16:     polygon_set_data<typename polygon_set_traits<polygon_set_type>::coordinate_type> ps;
Chris@16:     assign(ps, polygon_set);
Chris@16:     return ps.extents(extents_rectangle);
Chris@16:   }
Chris@16: 
Chris@16:   //area
Chris@16:   template <typename polygon_set_type>
Chris@16:   typename enable_if< typename is_mutable_polygon_set_type<polygon_set_type>::type,
Chris@16:                        typename coordinate_traits<typename polygon_set_traits<polygon_set_type>::coordinate_type>::area_type>::type
Chris@16:   area(const polygon_set_type& polygon_set) {
Chris@16:     typedef typename polygon_set_traits<polygon_set_type>::coordinate_type Unit;
Chris@16:     typedef polygon_with_holes_data<Unit> p_type;
Chris@16:     typedef typename coordinate_traits<Unit>::area_type area_type;
Chris@16:     std::vector<p_type> polys;
Chris@16:     assign(polys, polygon_set);
Chris@16:     area_type retval = (area_type)0;
Chris@16:     for(std::size_t i = 0; i < polys.size(); ++i) {
Chris@16:       retval += area(polys[i]);
Chris@16:     }
Chris@16:     return retval;
Chris@16:   }
Chris@16: 
Chris@16:   template <typename polygon_set_type>
Chris@16:   typename enable_if< typename is_mutable_polygon_set_type<polygon_set_type>::type,
Chris@16:                       std::size_t>::type
Chris@16:   simplify(polygon_set_type& polygon_set, typename coordinate_traits<
Chris@16:            typename polygon_set_traits<polygon_set_type>::coordinate_type
Chris@16:            >::coordinate_distance threshold) {
Chris@16:     typedef typename polygon_set_traits<polygon_set_type>::coordinate_type Unit;
Chris@16:     typedef polygon_with_holes_data<Unit> p_type;
Chris@16:     std::vector<p_type> polys;
Chris@16:     assign(polys, polygon_set);
Chris@16:     std::size_t retval = 0;
Chris@16:     for(std::size_t i = 0; i < polys.size(); ++i) {
Chris@16:       retval += detail::simplify_detail::simplify(polys[i].self_.coords_,
Chris@16:                                                   polys[i].self_.coords_, threshold);
Chris@16:       for(typename std::list<polygon_data<Unit> >::iterator itrh =
Chris@16:             polys[i].holes_.begin(); itrh != (polys[i].holes_.end()); ++itrh) {
Chris@16:         retval += detail::simplify_detail::simplify((*itrh).coords_,
Chris@16:                                                     (*itrh).coords_, threshold);
Chris@16:       }
Chris@16:     }
Chris@16:     assign(polygon_set, polys);
Chris@16:     return retval;
Chris@16:   }
Chris@16: 
Chris@16:   template <typename polygon_set_type, typename coord_type>
Chris@16:   typename enable_if< typename is_mutable_polygon_set_type<polygon_set_type>::type,
Chris@16:                        polygon_set_type>::type &
Chris@16:   resize(polygon_set_type& polygon_set, coord_type resizing, bool corner_fill_arcs = false, int num_circle_segments = 0) {
Chris@16:     typedef typename polygon_set_traits<polygon_set_type>::coordinate_type Unit;
Chris@16:     clean(polygon_set);
Chris@16:     polygon_set_data<Unit> ps;
Chris@16:     assign(ps, polygon_set);
Chris@16:     ps.resize(resizing, corner_fill_arcs,num_circle_segments);
Chris@16:     assign(polygon_set, ps);
Chris@16:     return polygon_set;
Chris@16:   }
Chris@16: 
Chris@16:   template <typename polygon_set_type>
Chris@16:   typename enable_if< typename is_mutable_polygon_set_type<polygon_set_type>::type,
Chris@16:                        polygon_set_type>::type &
Chris@16:   bloat(polygon_set_type& polygon_set,
Chris@16:         typename coordinate_traits<typename polygon_set_traits<polygon_set_type>::coordinate_type>::unsigned_area_type bloating) {
Chris@16:     return resize(polygon_set, bloating);
Chris@16:   }
Chris@16: 
Chris@16:   template <typename polygon_set_type>
Chris@16:   typename enable_if< typename is_mutable_polygon_set_type<polygon_set_type>::type,
Chris@16:                        polygon_set_type>::type &
Chris@16:   shrink(polygon_set_type& polygon_set,
Chris@16:         typename coordinate_traits<typename polygon_set_traits<polygon_set_type>::coordinate_type>::unsigned_area_type shrinking) {
Chris@16:     return resize(polygon_set, -(typename polygon_set_traits<polygon_set_type>::coordinate_type)shrinking);
Chris@16:   }
Chris@16: 
Chris@16:   //interact
Chris@16:   template <typename polygon_set_type_1, typename polygon_set_type_2>
Chris@16:   typename enable_if< typename gtl_and_3 <
Chris@16:     typename is_any_polygon_set_type<polygon_set_type_1>::type,
Chris@16:     typename is_any_polygon_set_type<polygon_set_type_2>::type,
Chris@16:     typename is_either_polygon_set_type<polygon_set_type_1, polygon_set_type_2>::type>::type,
Chris@16:     polygon_set_type_1>::type&
Chris@16:   interact(polygon_set_type_1& polygon_set_1, const polygon_set_type_2& polygon_set_2) {
Chris@16:     polygon_set_data<typename polygon_set_traits<polygon_set_type_1>::coordinate_type> ps1;
Chris@16:     assign(ps1, polygon_set_1);
Chris@16:     polygon_set_data<typename polygon_set_traits<polygon_set_type_2>::coordinate_type> ps2;
Chris@16:     assign(ps2, polygon_set_2);
Chris@16:     ps1.interact(ps2);
Chris@16:     assign(polygon_set_1, ps1);
Chris@16:     return polygon_set_1;
Chris@16:   }
Chris@16: 
Chris@16:   template <typename polygon_set_type>
Chris@16:   typename enable_if< typename is_mutable_polygon_set_type<polygon_set_type>::type,
Chris@16:                        polygon_set_type>::type &
Chris@16:   scale_up(polygon_set_type& polygon_set,
Chris@16:            typename coordinate_traits<typename polygon_set_traits<polygon_set_type>::coordinate_type>::unsigned_area_type factor) {
Chris@16:     typedef typename polygon_set_traits<polygon_set_type>::coordinate_type Unit;
Chris@16:     clean(polygon_set);
Chris@16:     polygon_set_data<Unit> ps;
Chris@16:     assign(ps, polygon_set);
Chris@16:     ps.scale_up(factor);
Chris@16:     assign(polygon_set, ps);
Chris@16:     return polygon_set;
Chris@16:   }
Chris@16: 
Chris@16:   template <typename polygon_set_type>
Chris@16:   typename enable_if< typename is_mutable_polygon_set_type<polygon_set_type>::type,
Chris@16:                        polygon_set_type>::type &
Chris@16:   scale_down(polygon_set_type& polygon_set,
Chris@16:              typename coordinate_traits<typename polygon_set_traits<polygon_set_type>::coordinate_type>::unsigned_area_type factor) {
Chris@16:     typedef typename polygon_set_traits<polygon_set_type>::coordinate_type Unit;
Chris@16:     clean(polygon_set);
Chris@16:     polygon_set_data<Unit> ps;
Chris@16:     assign(ps, polygon_set);
Chris@16:     ps.scale_down(factor);
Chris@16:     assign(polygon_set, ps);
Chris@16:     return polygon_set;
Chris@16:   }
Chris@16: 
Chris@16:   //transform
Chris@16:   template <typename polygon_set_type, typename transformation_type>
Chris@16:   typename enable_if< typename is_mutable_polygon_set_type<polygon_set_type>::type,
Chris@16:                        polygon_set_type>::type &
Chris@16:   transform(polygon_set_type& polygon_set,
Chris@16:             const transformation_type& transformation) {
Chris@16:     typedef typename polygon_set_traits<polygon_set_type>::coordinate_type Unit;
Chris@16:     clean(polygon_set);
Chris@16:     polygon_set_data<Unit> ps;
Chris@16:     assign(ps, polygon_set);
Chris@16:     ps.transform(transformation);
Chris@16:     assign(polygon_set, ps);
Chris@16:     return polygon_set;
Chris@16:   }
Chris@16: 
Chris@16:   //keep
Chris@16:   template <typename polygon_set_type>
Chris@16:   typename enable_if< typename is_mutable_polygon_set_type<polygon_set_type>::type,
Chris@16:                        polygon_set_type>::type &
Chris@16:   keep(polygon_set_type& polygon_set,
Chris@16:        typename coordinate_traits<typename polygon_set_traits<polygon_set_type>::coordinate_type>::area_type min_area,
Chris@16:        typename coordinate_traits<typename polygon_set_traits<polygon_set_type>::coordinate_type>::area_type max_area,
Chris@16:        typename coordinate_traits<typename polygon_set_traits<polygon_set_type>::coordinate_type>::unsigned_area_type min_width,
Chris@16:        typename coordinate_traits<typename polygon_set_traits<polygon_set_type>::coordinate_type>::unsigned_area_type max_width,
Chris@16:        typename coordinate_traits<typename polygon_set_traits<polygon_set_type>::coordinate_type>::unsigned_area_type min_height,
Chris@16:        typename coordinate_traits<typename polygon_set_traits<polygon_set_type>::coordinate_type>::unsigned_area_type max_height) {
Chris@16:     typedef typename polygon_set_traits<polygon_set_type>::coordinate_type Unit;
Chris@16:     typedef typename coordinate_traits<Unit>::unsigned_area_type uat;
Chris@16:     std::list<polygon_with_holes_data<Unit> > polys;
Chris@16:     assign(polys, polygon_set);
Chris@16:     typename std::list<polygon_with_holes_data<Unit> >::iterator itr_nxt;
Chris@16:     for(typename std::list<polygon_with_holes_data<Unit> >::iterator itr = polys.begin(); itr != polys.end(); itr = itr_nxt){
Chris@16:       itr_nxt = itr;
Chris@16:       ++itr_nxt;
Chris@16:       rectangle_data<Unit> bbox;
Chris@16:       extents(bbox, *itr);
Chris@16:       uat pwidth = delta(bbox, HORIZONTAL);
Chris@16:       if(pwidth > min_width && pwidth <= max_width){
Chris@16:         uat pheight = delta(bbox, VERTICAL);
Chris@16:         if(pheight > min_height && pheight <= max_height){
Chris@16:           typename coordinate_traits<Unit>::area_type parea = area(*itr);
Chris@16:           if(parea <= max_area && parea >= min_area) {
Chris@16:             continue;
Chris@16:           }
Chris@16:         }
Chris@16:       }
Chris@16:       polys.erase(itr);
Chris@16:     }
Chris@16:     assign(polygon_set, polys);
Chris@16:     return polygon_set;
Chris@16:   }
Chris@16: 
Chris@16:   namespace operators {
Chris@16: 
Chris@16:   struct yes_ps_ob : gtl_yes {};
Chris@16: 
Chris@16:   template <typename geometry_type_1, typename geometry_type_2>
Chris@16:   typename enable_if< typename gtl_and_4 < yes_ps_ob, typename is_any_polygon_set_type<geometry_type_1>::type,
Chris@16:                                             typename is_any_polygon_set_type<geometry_type_2>::type,
Chris@16:                                             typename is_either_polygon_set_type<geometry_type_1, geometry_type_2>::type>::type,
Chris@16:                        polygon_set_view<geometry_type_1, geometry_type_2, 0> >::type
Chris@16:   operator|(const geometry_type_1& lvalue, const geometry_type_2& rvalue) {
Chris@16:     return polygon_set_view<geometry_type_1, geometry_type_2, 0>
Chris@16:       (lvalue, rvalue);
Chris@16:   }
Chris@16: 
Chris@16:   struct yes_ps_op : gtl_yes {};
Chris@16: 
Chris@16:   template <typename geometry_type_1, typename geometry_type_2>
Chris@16:   typename enable_if< typename gtl_and_4 < yes_ps_op,
Chris@16:     typename gtl_if<typename is_any_polygon_set_type<geometry_type_1>::type>::type,
Chris@16:     typename gtl_if<typename is_any_polygon_set_type<geometry_type_2>::type>::type,
Chris@16:     typename gtl_if<typename is_either_polygon_set_type<geometry_type_1, geometry_type_2>::type>::type>
Chris@16:     ::type, polygon_set_view<geometry_type_1, geometry_type_2, 0> >::type
Chris@16:   operator+(const geometry_type_1& lvalue, const geometry_type_2& rvalue) {
Chris@16:     return polygon_set_view<geometry_type_1, geometry_type_2, 0>
Chris@16:       (lvalue, rvalue);
Chris@16:   }
Chris@16: 
Chris@16:   struct yes_ps_os : gtl_yes {};
Chris@16: 
Chris@16:   template <typename geometry_type_1, typename geometry_type_2>
Chris@16:   typename enable_if< typename gtl_and_4 < yes_ps_os,
Chris@16:     typename is_any_polygon_set_type<geometry_type_1>::type,
Chris@16:     typename is_any_polygon_set_type<geometry_type_2>::type,
Chris@16:     typename is_either_polygon_set_type<geometry_type_1, geometry_type_2>::type>::type,
Chris@16:                        polygon_set_view<geometry_type_1, geometry_type_2, 1> >::type
Chris@16:   operator*(const geometry_type_1& lvalue, const geometry_type_2& rvalue) {
Chris@16:     return polygon_set_view<geometry_type_1, geometry_type_2, 1>
Chris@16:       (lvalue, rvalue);
Chris@16:   }
Chris@16: 
Chris@16:   struct yes_ps_oa : gtl_yes {};
Chris@16: 
Chris@16:   template <typename geometry_type_1, typename geometry_type_2>
Chris@16:   typename enable_if< typename gtl_and_4 < yes_ps_oa,
Chris@16:     typename is_any_polygon_set_type<geometry_type_1>::type,
Chris@16:     typename is_any_polygon_set_type<geometry_type_2>::type,
Chris@16:     typename is_either_polygon_set_type<geometry_type_1, geometry_type_2>::type>::type,
Chris@16:                        polygon_set_view<geometry_type_1, geometry_type_2, 1> >::type
Chris@16:   operator&(const geometry_type_1& lvalue, const geometry_type_2& rvalue) {
Chris@16:     return polygon_set_view<geometry_type_1, geometry_type_2, 1>
Chris@16:       (lvalue, rvalue);
Chris@16:   }
Chris@16: 
Chris@16:   struct yes_ps_ox : gtl_yes {};
Chris@16: 
Chris@16:   template <typename geometry_type_1, typename geometry_type_2>
Chris@16:   typename enable_if< typename gtl_and_4 < yes_ps_ox,
Chris@16:     typename is_any_polygon_set_type<geometry_type_1>::type,
Chris@16:     typename is_any_polygon_set_type<geometry_type_2>::type,
Chris@16:     typename is_either_polygon_set_type<geometry_type_1, geometry_type_2>::type>::type,
Chris@16:                        polygon_set_view<geometry_type_1, geometry_type_2, 2> >::type
Chris@16:   operator^(const geometry_type_1& lvalue, const geometry_type_2& rvalue) {
Chris@16:     return polygon_set_view<geometry_type_1, geometry_type_2, 2>
Chris@16:       (lvalue, rvalue);
Chris@16:   }
Chris@16: 
Chris@16:   struct yes_ps_om : gtl_yes {};
Chris@16: 
Chris@16:   template <typename geometry_type_1, typename geometry_type_2>
Chris@16:     typename enable_if< typename gtl_and_4 < yes_ps_om,
Chris@16:     typename gtl_if<typename is_any_polygon_set_type<geometry_type_1>::type>::type,
Chris@16:     typename gtl_if<typename is_any_polygon_set_type<geometry_type_2>::type>::type,
Chris@16:     typename gtl_if<typename is_either_polygon_set_type<geometry_type_1, geometry_type_2>::type>::type>
Chris@16:     ::type, polygon_set_view<geometry_type_1, geometry_type_2, 3> >::type
Chris@16:   operator-(const geometry_type_1& lvalue, const geometry_type_2& rvalue) {
Chris@16:     return polygon_set_view<geometry_type_1, geometry_type_2, 3>
Chris@16:       (lvalue, rvalue);
Chris@16:   }
Chris@16: 
Chris@16:   struct yes_ps_ope : gtl_yes {};
Chris@16: 
Chris@16:   template <typename geometry_type_1, typename geometry_type_2>
Chris@16:   typename enable_if< typename gtl_and_4< yes_ps_ope, gtl_yes, typename is_mutable_polygon_set_type<geometry_type_1>::type,
Chris@16:                                          typename is_any_polygon_set_type<geometry_type_2>::type>::type,
Chris@16:                        geometry_type_1>::type &
Chris@16:   operator+=(geometry_type_1& lvalue, const geometry_type_2& rvalue) {
Chris@16:     return self_assignment_boolean_op<geometry_type_1, geometry_type_2, 0>(lvalue, rvalue);
Chris@16:   }
Chris@16: 
Chris@16:   struct yes_ps_obe : gtl_yes {};
Chris@16: 
Chris@16:   template <typename geometry_type_1, typename geometry_type_2>
Chris@16:   typename enable_if< typename gtl_and_3< yes_ps_obe, typename is_mutable_polygon_set_type<geometry_type_1>::type,
Chris@16:                                          typename is_any_polygon_set_type<geometry_type_2>::type>::type,
Chris@16:                        geometry_type_1>::type &
Chris@16:   operator|=(geometry_type_1& lvalue, const geometry_type_2& rvalue) {
Chris@16:     return self_assignment_boolean_op<geometry_type_1, geometry_type_2, 0>(lvalue, rvalue);
Chris@16:   }
Chris@16: 
Chris@16:   struct yes_ps_ose : gtl_yes {};
Chris@16: 
Chris@16:   template <typename geometry_type_1, typename geometry_type_2>
Chris@16:   typename enable_if< typename gtl_and_3< yes_ps_ose, typename is_mutable_polygon_set_type<geometry_type_1>::type,
Chris@16:                                          typename is_any_polygon_set_type<geometry_type_2>::type>::type,
Chris@16:                        geometry_type_1>::type &
Chris@16:   operator*=(geometry_type_1& lvalue, const geometry_type_2& rvalue) {
Chris@16:     return self_assignment_boolean_op<geometry_type_1, geometry_type_2, 1>(lvalue, rvalue);
Chris@16:   }
Chris@16: 
Chris@16:   struct yes_ps_oae : gtl_yes {};
Chris@16: 
Chris@16:   template <typename geometry_type_1, typename geometry_type_2>
Chris@16:   typename enable_if<
Chris@16:     typename gtl_and_3< yes_ps_oae, typename is_mutable_polygon_set_type<geometry_type_1>::type,
Chris@16:                       typename is_any_polygon_set_type<geometry_type_2>::type>::type,
Chris@16:     geometry_type_1>::type &
Chris@16:   operator&=(geometry_type_1& lvalue, const geometry_type_2& rvalue) {
Chris@16:     return self_assignment_boolean_op<geometry_type_1, geometry_type_2, 1>(lvalue, rvalue);
Chris@16:   }
Chris@16: 
Chris@16:   struct yes_ps_oxe : gtl_yes {};
Chris@16: 
Chris@16:   template <typename geometry_type_1, typename geometry_type_2>
Chris@16:   typename enable_if< typename gtl_and_3< yes_ps_oxe, typename is_mutable_polygon_set_type<geometry_type_1>::type,
Chris@16:                                          typename is_any_polygon_set_type<geometry_type_2>::type>::type,
Chris@16:                        geometry_type_1>::type &
Chris@16:   operator^=(geometry_type_1& lvalue, const geometry_type_2& rvalue) {
Chris@16:     return self_assignment_boolean_op<geometry_type_1, geometry_type_2, 2>(lvalue, rvalue);
Chris@16:   }
Chris@16: 
Chris@16:   struct yes_ps_ome : gtl_yes {};
Chris@16: 
Chris@16:   template <typename geometry_type_1, typename geometry_type_2>
Chris@16:   typename enable_if<
Chris@16:     typename gtl_and_3< yes_ps_ome, typename is_mutable_polygon_set_type<geometry_type_1>::type,
Chris@16:                       typename is_any_polygon_set_type<geometry_type_2>::type>::type,
Chris@16:     geometry_type_1>::type &
Chris@16:   operator-=(geometry_type_1& lvalue, const geometry_type_2& rvalue) {
Chris@16:     return self_assignment_boolean_op<geometry_type_1, geometry_type_2, 3>(lvalue, rvalue);
Chris@16:   }
Chris@16: 
Chris@16:   // TODO: Dafna, test these four resizing operators
Chris@16:   struct y_ps_rpe : gtl_yes {};
Chris@16: 
Chris@16:   template <typename geometry_type_1, typename coordinate_type_1>
Chris@16:   typename enable_if< typename gtl_and_3< y_ps_rpe, typename is_mutable_polygon_set_type<geometry_type_1>::type,
Chris@16:                                          typename gtl_same_type<typename geometry_concept<coordinate_type_1>::type,
Chris@16:                                                                 coordinate_concept>::type>::type,
Chris@16:                        geometry_type_1>::type &
Chris@16:   operator+=(geometry_type_1& lvalue, coordinate_type_1 rvalue) {
Chris@16:     return resize(lvalue, rvalue);
Chris@16:   }
Chris@16: 
Chris@16:   struct y_ps_rme : gtl_yes {};
Chris@16: 
Chris@16:   template <typename geometry_type_1, typename coordinate_type_1>
Chris@16:   typename enable_if< typename gtl_and_3<y_ps_rme, typename gtl_if<typename is_mutable_polygon_set_type<geometry_type_1>::type>::type,
Chris@16:                                          typename gtl_same_type<typename geometry_concept<coordinate_type_1>::type,
Chris@16:                                                                 coordinate_concept>::type>::type,
Chris@16:                        geometry_type_1>::type &
Chris@16:   operator-=(geometry_type_1& lvalue, coordinate_type_1 rvalue) {
Chris@16:     return resize(lvalue, -rvalue);
Chris@16:   }
Chris@16: 
Chris@16:   struct y_ps_rp : gtl_yes {};
Chris@16: 
Chris@16:   template <typename geometry_type_1, typename coordinate_type_1>
Chris@16:   typename enable_if< typename gtl_and_3<y_ps_rp, typename gtl_if<typename is_mutable_polygon_set_type<geometry_type_1>::type>::type,
Chris@16:                                         typename gtl_same_type<typename geometry_concept<coordinate_type_1>::type,
Chris@16:                                                                coordinate_concept>::type>
Chris@16:   ::type, geometry_type_1>::type
Chris@16:   operator+(const geometry_type_1& lvalue, coordinate_type_1 rvalue) {
Chris@16:     geometry_type_1 retval(lvalue);
Chris@16:     retval += rvalue;
Chris@16:     return retval;
Chris@16:   }
Chris@16: 
Chris@16:   struct y_ps_rm : gtl_yes {};
Chris@16: 
Chris@16:   template <typename geometry_type_1, typename coordinate_type_1>
Chris@16:   typename enable_if< typename gtl_and_3<y_ps_rm, typename gtl_if<typename is_mutable_polygon_set_type<geometry_type_1>::type>::type,
Chris@16:                                         typename gtl_same_type<typename geometry_concept<coordinate_type_1>::type,
Chris@16:                                                                coordinate_concept>::type>
Chris@16:   ::type, geometry_type_1>::type
Chris@16:   operator-(const geometry_type_1& lvalue, coordinate_type_1 rvalue) {
Chris@16:     geometry_type_1 retval(lvalue);
Chris@16:     retval -= rvalue;
Chris@16:     return retval;
Chris@16:   }
Chris@16: 
Chris@16: 
Chris@16:   } //end operators namespace
Chris@16: 
Chris@16:   template <typename T>
Chris@16:   struct view_of<polygon_45_set_concept, T> {
Chris@16:     typedef typename get_coordinate_type<T, typename geometry_concept<T>::type >::type coordinate_type;
Chris@16:     T* tp;
Chris@16:     std::vector<polygon_45_with_holes_data<coordinate_type> > polys;
Chris@16:     view_of(const T& obj) : tp(), polys() {
Chris@16:       std::vector<polygon_with_holes_data<coordinate_type> > gpolys;
Chris@16:       assign(gpolys, obj);
Chris@16:       for(typename std::vector<polygon_with_holes_data<coordinate_type> >::iterator itr = gpolys.begin();
Chris@16:           itr != gpolys.end(); ++itr) {
Chris@16:         polys.push_back(polygon_45_with_holes_data<coordinate_type>());
Chris@16:         assign(polys.back(), view_as<polygon_45_with_holes_concept>(*itr));
Chris@16:       }
Chris@16:     }
Chris@16:     view_of(T& obj) : tp(&obj), polys() {
Chris@16:       std::vector<polygon_with_holes_data<coordinate_type> > gpolys;
Chris@16:       assign(gpolys, obj);
Chris@16:       for(typename std::vector<polygon_with_holes_data<coordinate_type> >::iterator itr = gpolys.begin();
Chris@16:           itr != gpolys.end(); ++itr) {
Chris@16:         polys.push_back(polygon_45_with_holes_data<coordinate_type>());
Chris@16:         assign(polys.back(), view_as<polygon_45_with_holes_concept>(*itr));
Chris@16:       }
Chris@16:     }
Chris@16: 
Chris@16:     typedef typename std::vector<polygon_45_with_holes_data<coordinate_type> >::const_iterator iterator_type;
Chris@16:     typedef view_of operator_arg_type;
Chris@16: 
Chris@16:     inline iterator_type begin() const {
Chris@16:       return polys.begin();
Chris@16:     }
Chris@16: 
Chris@16:     inline iterator_type end() const {
Chris@16:       return polys.end();
Chris@16:     }
Chris@16: 
Chris@16:     inline orientation_2d orient() const { return HORIZONTAL; }
Chris@16: 
Chris@16:     inline bool clean() const { return false; }
Chris@16: 
Chris@16:     inline bool sorted() const { return false; }
Chris@16: 
Chris@16:     inline T& get() { return *tp; }
Chris@16:   };
Chris@16: 
Chris@16:   template <typename T>
Chris@16:   struct polygon_45_set_traits<view_of<polygon_45_set_concept, T> > {
Chris@16:     typedef typename view_of<polygon_45_set_concept, T>::coordinate_type coordinate_type;
Chris@16:     typedef typename view_of<polygon_45_set_concept, T>::iterator_type iterator_type;
Chris@16:     typedef view_of<polygon_45_set_concept, T> operator_arg_type;
Chris@16: 
Chris@16:     static inline iterator_type begin(const view_of<polygon_45_set_concept, T>& polygon_set) {
Chris@16:       return polygon_set.begin();
Chris@16:     }
Chris@16: 
Chris@16:     static inline iterator_type end(const view_of<polygon_45_set_concept, T>& polygon_set) {
Chris@16:       return polygon_set.end();
Chris@16:     }
Chris@16: 
Chris@16:     static inline orientation_2d orient(const view_of<polygon_45_set_concept, T>& polygon_set) {
Chris@16:       return polygon_set.orient(); }
Chris@16: 
Chris@16:     static inline bool clean(const view_of<polygon_45_set_concept, T>& polygon_set) {
Chris@16:       return polygon_set.clean(); }
Chris@16: 
Chris@16:     static inline bool sorted(const view_of<polygon_45_set_concept, T>& polygon_set) {
Chris@16:       return polygon_set.sorted(); }
Chris@16: 
Chris@16:   };
Chris@16: 
Chris@16:   template <typename T>
Chris@16:   struct geometry_concept<view_of<polygon_45_set_concept, T> > {
Chris@16:     typedef polygon_45_set_concept type;
Chris@16:   };
Chris@16: 
Chris@16:   template <typename T>
Chris@16:   struct get_coordinate_type<view_of<polygon_45_set_concept, T>, polygon_45_set_concept> {
Chris@16:     typedef typename view_of<polygon_45_set_concept, T>::coordinate_type type;
Chris@16:   };
Chris@16:   template <typename T>
Chris@16:   struct get_iterator_type_2<view_of<polygon_45_set_concept, T>, polygon_45_set_concept> {
Chris@16:     typedef typename view_of<polygon_45_set_concept, T>::iterator_type type;
Chris@16:     static type begin(const view_of<polygon_45_set_concept, T>& t) { return t.begin(); }
Chris@16:     static type end(const view_of<polygon_45_set_concept, T>& t) { return t.end(); }
Chris@16:   };
Chris@16: }
Chris@16: }
Chris@16: #endif