/**********************************************************************

 *

 * GEOS - Geometry Engine Open Source

 * http://geos.osgeo.org

 *

 * Copyright (C) 2001-2002 Vivid Solutions Inc.

 *

 * This is free software; you can redistribute and/or modify it under

 * the terms of the GNU Lesser General Public Licence as published

 * by the Free Software Foundation.

 * See the COPYING file for more information.

 *

 **********************************************************************/

#include <geos/geom/Coordinate.h>

#include <geos/constants.h> // for std::isnan

#include <sstream>

#include <string>

#include <iomanip>

#include <cmath>

namespace geos {

namespace geom { // geos::geom

const CoordinateXY CoordinateXY::_nullCoord = CoordinateXY(DoubleNotANumber, DoubleNotANumber);

const Coordinate Coordinate::_nullCoord = Coordinate(DoubleNotANumber, DoubleNotANumber, DoubleNotANumber);

const CoordinateXYM CoordinateXYM::_nullCoord = CoordinateXYM(DoubleNotANumber, DoubleNotANumber, DoubleNotANumber);

const CoordinateXYZM CoordinateXYZM::_nullCoord = CoordinateXYZM(DoubleNotANumber, DoubleNotANumber, DoubleNotANumber, DoubleNotANumber);

const CoordinateXY&

CoordinateXY::getNull()

{

    return _nullCoord;

}

const Coordinate&

Coordinate::getNull()

{

    return _nullCoord;

}

const CoordinateXYM&

CoordinateXYM::getNull()

{

    return _nullCoord;

}

const CoordinateXYZM&

CoordinateXYZM::getNull()

{

    return _nullCoord;

}

std::string

Coordinate::toString() const

{

    std::ostringstream s;

    s << std::setprecision(17) << *this;

    return s.str();

}

std::string

CoordinateXY::toString() const

{

    std::ostringstream s;

    s << std::setprecision(17) << *this;

    return s.str();

}

std::string

CoordinateXYM::toString() const

{

    std::ostringstream s;

    s << std::setprecision(17) << *this;

    return s.str();

}

std::string

CoordinateXYZM::toString() const

{

    std::ostringstream s;

    s << std::setprecision(17) << *this;

    return s.str();

}

std::ostream&

operator<< (std::ostream& os, const CoordinateType typ)

{

    switch(typ) {

        case CoordinateType::XY: os << "XY"; break;

        case CoordinateType::XYZ: os << "XYZ"; break;

        case CoordinateType::XYM: os << "XYM"; break;

        case CoordinateType::XYZM: os << "XYZM"; break;

    }

    return os;

}

std::ostream&

operator<< (std::ostream& os, const CoordinateXY& c)

{

    os << c.x << " " << c.y;

    return os;

}

std::ostream&

operator<< (std::ostream& os, const Coordinate& c)

{

    os << c.x << " " << c.y;

    if(!std::isnan(c.z)) {

        os << " " << c.z;

    }

    return os;

}

std::ostream&

operator<< (std::ostream& os, const CoordinateXYM& c)

{

    os << c.x << " " << c.y;

    if(!std::isnan(c.m)) {

        os << " " << c.m;

    }

    return os;

}

std::ostream&

operator<< (std::ostream& os, const CoordinateXYZM& c)

{

    os << c.x << " " << c.y;

    if(!std::isnan(c.z) || !std::isnan(c.m)) {

        os << " " << c.z;

    }

    if(!std::isnan(c.m)) {

        os << " " << c.m;

    }

    return os;

}

} // namespace geos::geom

} // namespace geos