Skip to content

Instantly share code, notes, and snippets.

@sjhardy
sjhardy / boostgeomtestcase.cpp
Created May 17, 2015 12:47
Incorrect result for Boost Geometry line intersection for spherical_equatorial over pole (Boost 1.58.0)
#include <iostream>
#include <boost/geometry.hpp>
int main(int argc, const char *argv[])
{
typedef boost::geometry::model::point<double, 2, boost::geometry::cs::spherical_equatorial<boost::geometry::degree> > spoint;
typedef boost::geometry::model::linestring<spoint> sline;
sline sline1, sline2;