256 lines
9.2 KiB
C++
256 lines
9.2 KiB
C++
// Boost.Geometry (aka GGL, Generic Geometry Library)
|
|
// Unit test
|
|
|
|
// Copyright (c) 2015, Oracle and/or its affiliates.
|
|
|
|
// Licensed under the Boost Software License version 1.0.
|
|
// http://www.boost.org/users/license.html
|
|
|
|
// Contributed and/or modified by Menelaos Karavelas, on behalf of Oracle
|
|
|
|
#ifndef BOOST_TEST_MODULE
|
|
#define BOOST_TEST_MODULE test_equals_on_spheroid
|
|
#endif
|
|
|
|
#include <iostream>
|
|
|
|
#include <boost/test/included/unit_test.hpp>
|
|
|
|
#include "test_equals.hpp"
|
|
|
|
#include <boost/geometry/geometries/geometries.hpp>
|
|
|
|
#include <boost/geometry/core/cs.hpp>
|
|
|
|
namespace bgm = bg::model;
|
|
|
|
template <typename P1, typename P2 = P1>
|
|
struct test_point_point
|
|
{
|
|
static inline void apply(std::string const& header)
|
|
{
|
|
std::string const str = header + "-";
|
|
|
|
test_geometry<P1, P2>(str + "pp_01", "POINT(0 0)", "POINT(0 0)", true);
|
|
test_geometry<P1, P2>(str + "pp_02", "POINT(0 0)", "POINT(10 0)", false);
|
|
|
|
// points whose longitudes differ by 360 degrees
|
|
test_geometry<P1, P2>(str + "pp_03", "POINT(0 0)", "POINT(360 0)", true);
|
|
test_geometry<P1, P2>(str + "pp_04", "POINT(10 0)", "POINT(370 0)", true);
|
|
test_geometry<P1, P2>(str + "pp_05", "POINT(10 0)", "POINT(-350 0)", true);
|
|
test_geometry<P1, P2>(str + "pp_06", "POINT(180 10)", "POINT(-180 10)", true);
|
|
test_geometry<P1, P2>(str + "pp_06a", "POINT(540 10)", "POINT(-540 10)", true);
|
|
|
|
#ifdef BOOST_GEOMETRY_NORMALIZE_LATITUDE
|
|
test_geometry<P1, P2>(str + "pp_06b", "POINT(540 370)", "POINT(-540 -350)", true);
|
|
test_geometry<P1, P2>(str + "pp_06c", "POINT(1260 370)", "POINT(-1260 -350)", true);
|
|
test_geometry<P1, P2>(str + "pp_06d", "POINT(2340 370)", "POINT(-2340 -350)", true);
|
|
#endif
|
|
|
|
test_geometry<P1, P2>(str + "pp_06e", "POINT(-180 10)", "POINT(-540 10)", true);
|
|
test_geometry<P1, P2>(str + "pp_06f", "POINT(180 10)", "POINT(-540 10)", true);
|
|
|
|
// north & south pole
|
|
test_geometry<P1, P2>(str + "pp_07", "POINT(0 90)", "POINT(0 90)", true);
|
|
|
|
#ifdef BOOST_GEOMETRY_NORMALIZE_LATITUDE
|
|
test_geometry<P1, P2>(str + "pp_07a", "POINT(0 450)", "POINT(10 -270)", true);
|
|
test_geometry<P1, P2>(str + "pp_07b", "POINT(0 270)", "POINT(10 90)", false);
|
|
test_geometry<P1, P2>(str + "pp_07c", "POINT(0 -450)", "POINT(10 90)", false);
|
|
#endif
|
|
|
|
test_geometry<P1, P2>(str + "pp_08", "POINT(0 90)", "POINT(10 90)", true);
|
|
test_geometry<P1, P2>(str + "pp_09", "POINT(0 90)", "POINT(0 -90)", false);
|
|
test_geometry<P1, P2>(str + "pp_10", "POINT(0 -90)", "POINT(0 -90)", true);
|
|
test_geometry<P1, P2>(str + "pp_11", "POINT(0 -90)", "POINT(10 -90)", true);
|
|
test_geometry<P1, P2>(str + "pp_11a", "POINT(0 -90)", "POINT(10 90)", false);
|
|
test_geometry<P1, P2>(str + "pp_12", "POINT(0 -90)", "POINT(0 -85)", false);
|
|
test_geometry<P1, P2>(str + "pp_13", "POINT(0 90)", "POINT(0 85)", false);
|
|
test_geometry<P1, P2>(str + "pp_14", "POINT(0 90)", "POINT(10 85)", false);
|
|
|
|
// symmetric wrt prime meridian
|
|
test_geometry<P1, P2>(str + "pp_15", "POINT(-10 45)", "POINT(10 45)", false);
|
|
test_geometry<P1, P2>(str + "pp_16", "POINT(-170 45)", "POINT(170 45)", false);
|
|
|
|
// other points
|
|
test_geometry<P1, P2>(str + "pp_17", "POINT(-10 45)", "POINT(10 -45)", false);
|
|
test_geometry<P1, P2>(str + "pp_18", "POINT(-10 -45)", "POINT(10 45)", false);
|
|
test_geometry<P1, P2>(str + "pp_19", "POINT(10 -135)", "POINT(10 45)", false);
|
|
|
|
#ifdef BOOST_GEOMETRY_NORMALIZE_LATITUDE
|
|
test_geometry<P1, P2>(str + "pp_20", "POINT(190 135)", "POINT(10 45)", true);
|
|
test_geometry<P1, P2>(str + "pp_21", "POINT(190 150)", "POINT(10 30)", true);
|
|
test_geometry<P1, P2>(str + "pp_21a", "POINT(-170 150)", "POINT(10 30)", true);
|
|
test_geometry<P1, P2>(str + "pp_22", "POINT(190 -135)", "POINT(10 -45)", true);
|
|
test_geometry<P1, P2>(str + "pp_23", "POINT(190 -150)", "POINT(10 -30)", true);
|
|
test_geometry<P1, P2>(str + "pp_23a", "POINT(-170 -150)", "POINT(10 -30)", true);
|
|
#endif
|
|
}
|
|
};
|
|
|
|
|
|
template <typename P1, typename P2 = P1>
|
|
struct test_point_point_with_height
|
|
{
|
|
static inline void apply(std::string const& header)
|
|
{
|
|
std::string const str = header + "-";
|
|
|
|
test_geometry<P1, P2>(str + "pp_01",
|
|
"POINT(0 0 10)",
|
|
"POINT(0 0 20)",
|
|
true);
|
|
|
|
test_geometry<P1, P2>(str + "pp_02",
|
|
"POINT(0 0 10)",
|
|
"POINT(10 0 10)",
|
|
false);
|
|
|
|
// points whose longitudes differ by 360 degrees
|
|
test_geometry<P1, P2>(str + "pp_03",
|
|
"POINT(0 0 10)",
|
|
"POINT(360 0 10)",
|
|
true);
|
|
|
|
// points whose longitudes differ by 360 degrees
|
|
test_geometry<P1, P2>(str + "pp_04",
|
|
"POINT(10 0 10)",
|
|
"POINT(370 0 10)",
|
|
true);
|
|
|
|
test_geometry<P1, P2>(str + "pp_05",
|
|
"POINT(10 0 10)",
|
|
"POINT(10 0 370)",
|
|
false);
|
|
}
|
|
};
|
|
|
|
|
|
template <typename P>
|
|
void test_segment_segment(std::string const& header)
|
|
{
|
|
typedef bgm::segment<P> seg;
|
|
|
|
std::string const str = header + "-";
|
|
|
|
test_geometry<seg, seg>(str + "ss_01",
|
|
"SEGMENT(10 0,180 0)",
|
|
"SEGMENT(10 0,-180 0)",
|
|
true);
|
|
test_geometry<seg, seg>(str + "ss_02",
|
|
"SEGMENT(0 90,180 0)",
|
|
"SEGMENT(10 90,-180 0)",
|
|
true);
|
|
test_geometry<seg, seg>(str + "ss_03",
|
|
"SEGMENT(0 90,0 -90)",
|
|
"SEGMENT(10 90,20 -90)",
|
|
true);
|
|
test_geometry<seg, seg>(str + "ss_04",
|
|
"SEGMENT(10 80,10 -80)",
|
|
"SEGMENT(10 80,20 -80)",
|
|
false);
|
|
test_geometry<seg, seg>(str + "ss_05",
|
|
"SEGMENT(170 10,-170 10)",
|
|
"SEGMENT(170 10,350 10)",
|
|
false);
|
|
}
|
|
|
|
|
|
BOOST_AUTO_TEST_CASE( equals_point_point_se )
|
|
{
|
|
typedef bg::cs::spherical_equatorial<bg::degree> cs_type;
|
|
|
|
test_point_point<bgm::point<int, 2, cs_type> >::apply("se");
|
|
test_point_point<bgm::point<double, 2, cs_type> >::apply("se");
|
|
test_point_point<bgm::point<long double, 2, cs_type> >::apply("se");
|
|
|
|
// mixed point types
|
|
test_point_point
|
|
<
|
|
bgm::point<double, 2, cs_type>, bgm::point<int, 2, cs_type>
|
|
>::apply("se");
|
|
|
|
test_point_point
|
|
<
|
|
bgm::point<double, 2, cs_type>, bgm::point<long double, 2, cs_type>
|
|
>::apply("se");
|
|
|
|
#if defined(HAVE_TTMATH)
|
|
test_point_point<bgm::point<ttmath_big, 2, cs_type> >::apply("se");
|
|
#endif
|
|
}
|
|
|
|
BOOST_AUTO_TEST_CASE( equals_point_point_with_height_se )
|
|
{
|
|
typedef bg::cs::spherical_equatorial<bg::degree> cs_type;
|
|
|
|
test_point_point<bgm::point<int, 3, cs_type> >::apply("seh");
|
|
test_point_point<bgm::point<double, 3, cs_type> >::apply("seh");
|
|
test_point_point<bgm::point<long double, 3, cs_type> >::apply("seh");
|
|
|
|
// mixed point types
|
|
test_point_point
|
|
<
|
|
bgm::point<double, 3, cs_type>, bgm::point<int, 3, cs_type>
|
|
>::apply("seh");
|
|
|
|
test_point_point
|
|
<
|
|
bgm::point<double, 3, cs_type>, bgm::point<long double, 3, cs_type>
|
|
>::apply("seh");
|
|
|
|
#if defined(HAVE_TTMATH)
|
|
test_point_point<bgm::point<ttmath_big, 3, cs_type> >::apply("seh");
|
|
#endif
|
|
}
|
|
|
|
BOOST_AUTO_TEST_CASE( equals_point_point_geo )
|
|
{
|
|
typedef bg::cs::geographic<bg::degree> cs_type;
|
|
|
|
test_point_point<bgm::point<int, 2, cs_type> >::apply("geo");
|
|
test_point_point<bgm::point<double, 2, cs_type> >::apply("geo");
|
|
test_point_point<bgm::point<long double, 2, cs_type> >::apply("geo");
|
|
|
|
// mixed point types
|
|
test_point_point
|
|
<
|
|
bgm::point<double, 2, cs_type>, bgm::point<int, 2, cs_type>
|
|
>::apply("se");
|
|
|
|
test_point_point
|
|
<
|
|
bgm::point<double, 2, cs_type>, bgm::point<long double, 2, cs_type>
|
|
>::apply("se");
|
|
|
|
#if defined(HAVE_TTMATH)
|
|
test_point_point<bgm::point<ttmath_big, 2, cs_type> >::apply("geo");
|
|
#endif
|
|
}
|
|
|
|
BOOST_AUTO_TEST_CASE( equals_segment_segment_se )
|
|
{
|
|
typedef bg::cs::spherical_equatorial<bg::degree> cs_type;
|
|
|
|
test_segment_segment<bgm::point<int, 2, cs_type> >("se");
|
|
test_segment_segment<bgm::point<double, 2, cs_type> >("se");
|
|
test_segment_segment<bgm::point<long double, 2, cs_type> >("se");
|
|
|
|
#if defined(HAVE_TTMATH)
|
|
test_segment_segment<bgm::point<ttmath_big, 2, cs_type> >("se");
|
|
#endif
|
|
}
|
|
|
|
BOOST_AUTO_TEST_CASE( equals_segment_segment_geo )
|
|
{
|
|
typedef bg::cs::geographic<bg::degree> cs_type;
|
|
|
|
test_segment_segment<bgm::point<int, 2, cs_type> >("geo");
|
|
test_segment_segment<bgm::point<double, 2, cs_type> >("geo");
|
|
test_segment_segment<bgm::point<long double, 2, cs_type> >("geo");
|
|
|
|
#if defined(HAVE_TTMATH)
|
|
test_segment_segment<bgm::point<ttmath_big, 2, cs_type> >("geo");
|
|
#endif
|
|
}
|