/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
    \\  /    A nd           | Copyright (C) 1991-2010 OpenCFD Ltd.
     \\/     M anipulation  |
-------------------------------------------------------------------------------
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::plane

Description
    Geometric class that creates a 2D plane and can return the intersection
     point between a line and the plane.

SourceFiles
    plane.C

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

#ifndef plane_H
#define plane_H

#include <OpenFOAM/point.H>
#include <OpenFOAM/scalarList.H>
#include <OpenFOAM/dictionary.H>
#include <OpenFOAM/line.H>

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

namespace Foam
{

// Forward declaration of friend functions and operators

class plane;
bool operator==(const plane&, const plane&);
bool operator!=(const plane&, const plane&);
Ostream& operator<<(Ostream&, const plane&);


/*---------------------------------------------------------------------------*\
                           Class plane Declaration
\*---------------------------------------------------------------------------*/

class plane
{
public:

        //- A direction and a reference point
        class ray
        {
            point refPoint_;

            vector dir_;

        public:

            ray(const point& refPoint, const vector& dir)
            :
                refPoint_(refPoint),
                dir_(dir)
            {}

            const point& refPoint() const
            {
                return refPoint_;
            }

            const vector& dir() const
            {
                return dir_;
            }
        };


private:

    // Private data

        //- Plane normal
        vector unitVector_;

        //- Base point
        point basePoint_;


    // Private Member Functions

        //- Calculates basePoint and normal vector given plane coefficients
        void calcPntAndVec(const scalarList& C);

        //- Calculates basePoint and normal vector given three points
        //- Normal vector determined using right hand rule
        void calcPntAndVec
        (
            const point& point1,
            const point& point2,
            const point& point3
        );


public:

    // Constructors

        //- Construct from normal vector through the origin
        plane(const vector& normalVector);

        //- Construct from normal vector and point in plane
        plane(const point& basePoint, const vector& normalVector);

        //- Construct from three points in plane
        plane(const point& point1, const point& point2, const point& point3);

        //- Construct from coefficients for the
        //  plane equation: ax + by + cz + d = 0
        plane(const scalarList& C);

        //- Construct from dictionary
        plane(const dictionary& planeDict);

        //- Construct from Istream. Assumes the base + normal notation.
        plane(Istream& is);


    // Member Functions

        //- Return plane normal
        const vector& normal() const;

        //- Return or return plane base point
        const point& refPoint() const;

        //- Return coefficients for the
        //  plane equation: ax + by + cz + d = 0
        FixedList<scalar, 4> planeCoeffs() const;

        //- Return nearest point in the plane for the given point
        point nearestPoint(const point& p) const;

        //- Return distance from the given point to the plane
        scalar distance(const point& p) const;

        //- Return cut coefficient for plane and line defined by
        //  origin and direction
        scalar normalIntersect(const point& pnt0, const vector& dir) const;

        //- Return cut coefficient for plane and ray
        scalar normalIntersect(const ray& r) const
        {
            return normalIntersect(r.refPoint(), r.dir());
        }

        //- Return the cutting point between the plane and
        // a line passing through the supplied points
        template<class Point, class PointRef>
        scalar lineIntersect(const line<Point, PointRef>& l) const
        {
            return normalIntersect(l.start(), l.vec());
        }

        //- Return the cutting line between this plane and another.
        //  Returned as direction vector and point line goes through.
        ray planeIntersect(const plane&) const;

        //- Return the cutting point between this plane and two other planes
        point planePlaneIntersect(const plane&, const plane&) const;

        //- Write to dictionary
        void writeDict(Ostream&) const;


    // friend Operators

        friend bool operator==(const plane&, const plane&);
        friend bool operator!=(const plane&, const plane&);


    // IOstream Operators

        //- Write plane properties
        friend Ostream& operator<<(Ostream&, const plane&);
};


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

} // End namespace Foam

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

#endif

// ************************ vim: set sw=4 sts=4 et: ************************ //
