id,summary,reporter,owner,description,type,status,milestone,component,version,severity,resolution,keywords,cc
13553,intersection gives wrong result,marnix.berg@…,Barend Gehrels,"The following polygons result in a wrong intersection:

using point_type = boost::geometry::model::d2::point_xy<double>;
typedef boost::geometry::model::ring<point_type, false, true> polygon;

polygon op1, op2;

boost::geometry::read_wkt(""POLYGON((7.7058932076134967 -11.523889618379748,8.0348094747518424 0.63492733448631888,7.7720440433489850 0.63492733448631888, 7.7058932076134967 -11.523889618379748))"", op1);
 
boost::geometry::read_wkt(""POLYGON((2.6206910206017642 -32.773696844382265, 5.5835888947200090 -24.273798818378602, 6.7109368565951772 -20.023839227004206, 7.4191426214038723 -15.773870150408516, 7.7058941612236938 -11.523876267001913, -3.1025600674348395 -11.523582486001384, -3.1025610210450365 -32.773541282571529, 2.6206910206017642 -32.773696844382265))"", op2); 
		
std::vector<polygon> result;
boost::geometry::intersection(op1, op2, result);

result is equal to op1, while op1 is mostly outside op2.",Bugs,assigned,Boost 1.68.0,geometry,Boost 1.66.0,Problem,,,henrik@…