// Boost.Geometry (aka GGL, Generic Geometry Library)
// Unit Test

// Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands.
// Copyright (c) 2008-2012 Bruno Lalande, Paris, France.
// Copyright (c) 2009-2012 Mateusz Loskot, London, UK.

// Parts of Boost.Geometry are redesigned from Geodan's Geographic Library
// (geolib/GGL), copyright (c) 1995-2010 Geodan, Amsterdam, the Netherlands.

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


#include <geometry_test_common.hpp>

#include <boost/geometry/extensions/algorithms/parse.hpp>
#include <boost/geometry/geometries/point.hpp>
#include <boost/geometry/strategies/strategies.hpp>

#include <boost/geometry/extensions/gis/geographic/strategies/dms_parser.hpp>

#include <test_common/test_point.hpp>

using namespace boost::geometry;

template <typename P>
void test_2d(double c, bool use_strategy)
{
    // normal order, east=x, north=y
    P p;
    parse(p, std::string("1dE"), std::string("2N"));
    BOOST_CHECK_CLOSE( ((double) bg::get<0>(p)), (double) 1 * c, 1.0e-6);
    BOOST_CHECK_CLOSE( ((double) bg::get<1>(p)), (double) 2 * c, 1.0e-6);

    // reversed order, y,x -> should be interpreted correctly
    parse(p, std::string("1dN"), std::string("2E"));
    BOOST_CHECK_CLOSE( ((double) bg::get<0>(p)), (double) 2 * c, 1.0e-6);
    BOOST_CHECK_CLOSE( ((double) bg::get<1>(p)), (double) 1 * c, 1.0e-6);

    if (use_strategy)
    {
        // DUTCH system NOZW, only for degrees
        bg::strategy::dms_parser<false, 'N', 'O', 'Z', 'W'> strategy;
        parse(p, std::string("1dO"), std::string("2Z"), strategy);
        BOOST_CHECK_CLOSE( ((double) bg::get<0>(p)), (double) 1, 1.0e-6);
        BOOST_CHECK_CLOSE( ((double) bg::get<1>(p)), (double) -2, 1.0e-6);
    }

    // rest of DMS is checked in parse_dms
}

template <typename T, typename P>
void test_3d()
{
}

int test_main(int, char* [])
{
    //test_2d<point<int, 2, cs::geographic<radian> > >();
    //test_2d<point<float, 2, cs::geographic<radian> > >();

    test_2d<bg::model::point<double, 2, bg::cs::geographic<bg::degree> > >(1.0, true);
    test_2d<bg::model::point<double, 2, bg::cs::geographic<bg::radian> > >(bg::math::d2r, false);

    return 0;
}