#include <tut/tut.hpp>
// geos
#include <geos_c.h>

#include "capi_test_utils.h"

namespace tut {
//
// Test Group
//

struct test_geosboundary_data : public capitest::utility {};

typedef test_group<test_geosboundary_data> group;
typedef group::object object;

group test_geosboundary("capi::GEOSBoundary");

template<>
template<>
void object::test<1>()
{
    input_ = fromWKT("POLYGON ((0 0, 10 0, 10 10, 0 10, 0 0), (1 1, 1 2, 2 2, 2 1, 1 1))");
    GEOSSetSRID(input_, 3857);

    result_ = GEOSBoundary(input_);
    expected_ = fromWKT("MULTILINESTRING ((0 0, 10 0, 10 10, 0 10, 0 0), (1 1, 1 2, 2 2, 2 1, 1 1))");

    ensure_geometry_equals_exact(result_, expected_, 0.0);
    ensure_equals(GEOSGetSRID(input_), GEOSGetSRID(result_));
}

template<>
template<>
void object::test<2>()
{
    set_test_name("curved inputs");

    input_ = fromWKT("CIRCULARSTRING (0 0, 1 1, 2 0)");
    ensure(input_ != nullptr);

    result_ = GEOSBoundary(input_);
    ensure(result_);

    expected_ = fromWKT("MULTIPOINT (0 0, 2 0)");
    ensure_geometry_equals_identical(result_, expected_);
}

template<>
template<>
void object::test<3>()
{
    input_ = fromWKT("MULTILINESTRING M ((10 11 3, 20 21 4), (20 21 4, 32 21 3))");
    result_ = GEOSBoundary(input_);
    expected_ = fromWKT("MULTIPOINT M ((10 11 3), (32 21 3))");

    ensure_geometry_equals_identical(result_, expected_);
}

template<>
template<>
void object::test<4>()
{
    input_ = fromWKT("POLYGON M ((0 0 0, 1 0 1, 1 1 2, 0 1 3, 0 0 4))");
    result_ = GEOSBoundary(input_);
    expected_ = fromWKT("LINESTRING M (0 0 0, 1 0 1, 1 1 2, 0 1 3, 0 0 4)");

    ensure_geometry_equals_identical(result_, expected_);
}

} // namespace tut

