#include <iostream>
#include <boost/geometry.hpp>
#include <boost/geometry/io/wkt/wkt.hpp>



namespace bg = boost::geometry;

void intersectsTest(size_t maxtimes) {
  typedef bg::model::point<double, 2, bg::cs::spherical_equatorial<bg::degree> > SEP;
  typedef boost::geometry::model::polygon<SEP> polygon_type;
  polygon_type testpl;
  SEP testp;
  bool r = false;
  /// wkt for spherical equatorial, format is "longtitude latitude"
  boost::geometry::read_wkt("POLYGON((30 0,30 30,90 30, 90 0, 30 0))", testpl);
  boost::geometry::assign_values(testp, -120,21);
  std::cout << "Result 1: " << boost::geometry::intersects(testpl, testp) << std::endl;
  boost::geometry::read_wkt("POLYGON((0 80,-90 80, -180 80, 90 80, 0 80))", testpl);
  boost::geometry::assign_values(testp, 0,90);
  std::cout << "Result 2: " << boost::geometry::intersects(testpl, testp) << std::endl;
}