EQ2EMu/EQ2/source/depends/boost_1_72_0/libs/geometry/test/algorithms/similarity/test_hausdorff_distance.hpp

// Boost.Geometry

// Copyright (c) 2018 Yaghyavardhan Singh Khangarot, Hyderabad, India.
// Contributed and/or modified by Yaghyavardhan Singh Khangarot,
//   as part of Google Summer of Code 2018 program.

// This file was modified by Oracle on 2018.
// Modifications copyright (c) 2018, Oracle and/or its affiliates.
// Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle

// 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)

#ifndef BOOST_GEOMETRY_TEST_HAUSDORFF_DISTANCE_HPP
#define BOOST_GEOMETRY_TEST_HAUSDORFF_DISTANCE_HPP

#include <geometry_test_common.hpp>
#include <boost/geometry/algorithms/discrete_hausdorff_distance.hpp>
#include <boost/geometry/io/wkt/wkt.hpp>
#include <boost/geometry/strategies/strategies.hpp>
#include <boost/variant/variant.hpp>

template <typename Geometry1, typename Geometry2, typename Expected>
void test_hausdorff_distance(Geometry1 const& geometry1, Geometry2 const& geometry2,
                             Expected const& expected_hausdorff_distance)
{
    typedef typename bg::distance_result
        <
            typename bg::point_type<Geometry1>::type,
            typename bg::point_type<Geometry2>::type
        >::type result_type;

    result_type h_distance = bg::discrete_hausdorff_distance(geometry1, geometry2);

#ifdef BOOST_GEOMETRY_TEST_DEBUG
    std::ostringstream out;
    out << typeid(typename bg::coordinate_type<Geometry1>::type).name()
        << std::endl
        << typeid(typename bg::coordinate_type<Geometry2>::type).name()
        << std::endl
        << typeid(h_distance).name()
        << std::endl
        << "hausdorff_distance : " << h_distance
        << std::endl;
    std::cout << out.str();
#endif

    BOOST_CHECK_CLOSE(h_distance, result_type(expected_hausdorff_distance), 0.01);
}



template <typename Geometry1, typename Geometry2, typename Expected>
void test_geometry(std::string const& wkt1, std::string const& wkt2,
                   Expected const& expected_hausdorff_distance)
{
    Geometry1 geometry1;
    bg::read_wkt(wkt1, geometry1);
    Geometry2 geometry2;
    bg::read_wkt(wkt2, geometry2);

    test_hausdorff_distance(geometry1, geometry2, expected_hausdorff_distance);

#if defined(BOOST_GEOMETRY_TEST_DEBUG)
    test_hausdorff_distance(boost::variant<Geometry1>(geometry1),
                            boost::variant<Geometry2>(geometry2),
                            expected_hausdorff_distance);
#endif
}

template <typename Geometry1, typename Geometry2, typename Strategy, typename Expected>
void test_hausdorff_distance(Geometry1 const& geometry1, Geometry2 const& geometry2,
                             Strategy strategy, Expected const& expected_hausdorff_distance)
{
    typedef typename bg::distance_result
        <
            typename bg::point_type<Geometry1>::type,
            typename bg::point_type<Geometry2>::type,
            Strategy
        >::type result_type;

    result_type h_distance = bg::discrete_hausdorff_distance(geometry1, geometry2, strategy);

#ifdef BOOST_GEOMETRY_TEST_DEBUG
    std::ostringstream out;
    out << typeid(typename bg::coordinate_type<Geometry1>::type).name()
        << std::endl
        << typeid(typename bg::coordinate_type<Geometry2>::type).name()
        << std::endl
        << typeid(h_distance).name()
        << std::endl
        << "hausdorff_distance : " << h_distance
        << std::endl;
    std::cout << out.str();
#endif

    BOOST_CHECK_CLOSE(h_distance, result_type(expected_hausdorff_distance), 0.01);
}



template <typename Geometry1, typename Geometry2, typename Strategy, typename Expected>
void test_geometry(std::string const& wkt1, std::string const& wkt2,
                   Strategy strategy, Expected const& expected_hausdorff_distance)
{
    Geometry1 geometry1;
    bg::read_wkt(wkt1, geometry1);
    Geometry2 geometry2;
    bg::read_wkt(wkt2, geometry2);

    test_hausdorff_distance(geometry1, geometry2, strategy, expected_hausdorff_distance);

#if defined(BOOST_GEOMETRY_TEST_DEBUG)
    test_hausdorff_distance(boost::variant<Geometry1>(geometry1),
                            boost::variant<Geometry2>(geometry2),
                            strategy, expected_hausdorff_distance);
#endif
}


template <typename Geometry1,typename Geometry2>
void test_empty_input(Geometry1 const& geometry1, Geometry2 const& geometry2)
{
    try
    {
        bg::discrete_hausdorff_distance(geometry1, geometry2);
    }
    catch(bg::empty_input_exception const& )
    {
        return;
    }
    BOOST_CHECK_MESSAGE(false, "A empty_input_exception should have been thrown");
}


#endif