#include <boost/ggl/ggl.hpp>
#include <boost/ggl/algorithms/union.hpp>
#include <boost/ggl/geometries/adapted/c_array_cartesian.hpp>
#include <boost/ggl/geometries/cartesian2d.hpp>
#include <boost/ggl/multi/geometries/multi_polygon.hpp>
#include <iostream>

using ggl::assign;
using ggl::correct;
using ggl::dsv;
using ggl::polygon_2d;
using ggl::multi_polygon;
using ggl::union_inserter;
using std::cout;

typedef multi_polygon<polygon_2d> polygon_set;

int main()
{
  polygon_2d poly1;
  {
    const double coor[][2] =
      { {0.0, 0.0}, {0.0, 2.0}, {2.0, 2.0}, {2.0, 0.0}, {0.0, 0.0} };
    assign(poly1, coor);
  }
  polygon_2d poly2;
  {
    const double coor[][2] =
      { {0.5, 1.5}, {0.5, 2.5}, {1.5, 2.5}, {1.5, 1.5}, {0.5, 1.5} };
    assign(poly2, coor);
  }

  cout << "poly1: " << dsv(poly1) << "\n";
  cout << "poly2: " << dsv(poly2) << "\n";
  polygon_set ps1;
  union_inserter<polygon_2d>(poly1, poly2, std::back_inserter(ps1));
  cout << "poly1|poly2:\n";
  for (polygon_set::const_iterator it = ps1.begin(); it != ps1.end(); ++it)
    cout << dsv(*it) << "\n";

  return EXIT_SUCCESS;
}