/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
    \\  /    A nd           | www.openfoam.com
     \\/     M anipulation  |
-------------------------------------------------------------------------------
    Copyright (C) 2012-2016 OpenFOAM Foundation
-------------------------------------------------------------------------------
License
    This file is part of OpenFOAM.

    OpenFOAM is free software: you can redistribute it and/or modify it
    under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
    for more details.

    You should have received a copy of the GNU General Public License
    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.

Class
    Foam::pointConversion

Description
    Conversion functions between point (Foam::) and Point (CGAL::)

\*---------------------------------------------------------------------------*/

#ifndef pointConversion_H
#define pointConversion_H

#include "point.H"

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

#ifdef CGAL_INEXACT

    // Define Point to be contiguous for inexact (double storage) kernel
    typedef const Foam::point& pointFromPoint;
    typedef const CGAL::Triangulation_vertex_base_3<K>::Point& PointFrompoint;

#else

    typedef Foam::point pointFromPoint;
    typedef CGAL::Triangulation_vertex_base_3<K>::Point PointFrompoint;

#endif

namespace Foam
{

#ifdef CGAL_INEXACT

    template<class Point>
    inline pointFromPoint topoint(const Point& P)
    {
        return reinterpret_cast<pointFromPoint>(P);
    }

    inline PointFrompoint toPoint(const Foam::point& p)
    {
        return reinterpret_cast<PointFrompoint>(p);
    }

#else

    template<class Point>
    inline pointFromPoint topoint(const Point& P)
    {
        return Foam::point
        (
            CGAL::to_double(P.x()),
            CGAL::to_double(P.y()),
            CGAL::to_double(P.z())
        );
    }

    inline PointFrompoint toPoint(const Foam::point& p)
    {
        return PointFrompoint(p.x(), p.y(), p.z());
    }

#endif

    //- Specialisation for indexedVertex.
    template<>
    inline pointFromPoint topoint<CGAL::indexedVertex<K>>
    (
        const CGAL::indexedVertex<K>& P
    )
    {
        return topoint(P.point());
    }

    //- Specialisation for Foam::point. Used only as a dummy.
    template<>
    inline pointFromPoint topoint<Foam::point>
    (
        const Foam::point& P
    )
    {
        return P;
    }

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

} // End namespace Foam

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

#endif

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