point_on_surface bug, computed point isn't on the surface of polygon

[petr@…]

The attached test program demonstrates this bug.

[petr@…]

I do not have permission to add attachments. Code is below in this comment.

#include <boost/geometry.hpp>
#include <boost/geometry/geometries/point_xy.hpp>
#include <boost/geometry/geometries/polygon.hpp>

#include <assert.h>

void test_point_on_surface()
{
    typedef boost::geometry::model::point<double, 2, boost::geometry::cs::cartesian> Point;
    boost::geometry::model::polygon<Point> poly;

    boost::geometry::read_wkt("POLYGON((1074699.93 703064.65, 1074702.10 703054.62, 1074704.53 703061.40, 1074703.90 703064.58, 1074699.93 703064.65))",poly);
    boost::geometry::correct(poly);

    Point point;
    boost::geometry::point_on_surface(poly,point);

    //computed point isn't on the surface of polygon
    //[0]	1074839.1013086310	double
		//[1]	703154.77900308778	double

    bool res = boost::geometry::within(point,poly);
    assert(res==true);
    // res is false :-(
}

int main(int argc, char* argv[])
{
  test_point_on_surface();
  return 0;
}