geometry/index/test/algorithms/intersection_content.cpp
2014-06-19 12:36:21 +03:00

71 lines
2.9 KiB
C++

// Boost.Geometry Index
// Unit Test
// Copyright (c) 2011-2013 Adam Wulkiewicz, Lodz, Poland.
// Use, modification and distribution is subject to the Boost Software License,
// Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
// http://www.boost.org/LICENSE_1_0.txt)
#include <algorithms/test_intersection_content.hpp>
#include <boost/geometry/geometries/point_xy.hpp>
#include <boost/geometry/geometries/point.hpp>
#include <boost/geometry/geometries/box.hpp>
//#define BOOST_GEOMETRY_TEST_DEBUG
void test_large_integers()
{
typedef bg::model::point<int, 2, bg::cs::cartesian> int_point_type;
typedef bg::model::point<double, 2, bg::cs::cartesian> double_point_type;
bg::model::box<int_point_type> int_box1, int_box2;
bg::model::box<double_point_type> double_box1, double_box2;
std::string const box_li1 = "POLYGON((1536119 192000, 1872000 528000))";
std::string const box_li2 = "POLYGON((1701234 368250, 2673400 777400))";
bg::read_wkt(box_li1, int_box1);
bg::read_wkt(box_li1, double_box1);
bg::read_wkt(box_li2, int_box2);
bg::read_wkt(box_li2, double_box2);
double int_value = bgi::detail::intersection_content(int_box1, int_box2);
double double_value = bgi::detail::intersection_content(double_box1, double_box2);
BOOST_CHECK_CLOSE(int_value, double_value, 0.0001);
}
int test_main(int, char* [])
{
typedef bg::model::point<int, 2, bg::cs::cartesian> P2ic;
typedef bg::model::point<float, 2, bg::cs::cartesian> P2fc;
typedef bg::model::point<double, 2, bg::cs::cartesian> P2dc;
typedef bg::model::point<int, 3, bg::cs::cartesian> P3ic;
typedef bg::model::point<float, 3, bg::cs::cartesian> P3fc;
typedef bg::model::point<double, 3, bg::cs::cartesian> P3dc;
test_geometry<bg::model::box<P2ic> >("POLYGON((0 1,2 4))", "POLYGON((1 2,3 5))", 2.0);
test_geometry<bg::model::box<P2fc> >("POLYGON((0 1,2 4))", "POLYGON((1 2,3 5))", 2.0);
test_geometry<bg::model::box<P2dc> >("POLYGON((0 1,2 4))", "POLYGON((1 2,3 5))", 2.0);
test_geometry<bg::model::box<P3ic> >("POLYGON((0 1 2,2 4 6))", "POLYGON((1 2 3,3 5 7))", 6.0);
test_geometry<bg::model::box<P3fc> >("POLYGON((0 1 2,2 4 6))", "POLYGON((1 2 3,3 5 7))", 6.0);
test_geometry<bg::model::box<P3dc> >("POLYGON((0 1 2,2 4 6))", "POLYGON((1 2 3,3 5 7))", 6.0);
test_geometry<bg::model::box<P2dc> >("POLYGON((0 1,2 4))", "POLYGON((2 1,3 4))", 0.0);
test_geometry<bg::model::box<P2dc> >("POLYGON((0 1,2 4))", "POLYGON((2 4,3 5))", 0.0);
#ifdef HAVE_TTMATH
typedef bg::model::point<ttmath_big, 2, bg::cs::cartesian> P2ttmc;
typedef bg::model::point<ttmath_big, 3, bg::cs::cartesian> P3ttmc;
test_geometry<bg::model::box<P2ttmc> >("POLYGON((0 1,2 4))", "POLYGON((1 2,3 5))", 2.0);
test_geometry<bg::model::box<P3ttmc> >("POLYGON((0 1 2,2 4 6))", "POLYGON((1 2 3,3 5 7))", 6.0);
#endif
test_large_integers();
return 0;
}