Plane.hpp

Go to the documentation of this file.

// This file is part of GranOO, a workbench for DEM simulation.
//
// Author(s)    : - Damien Andre       IRCER/UNILIM, Limoges France
//                  <damien.andre@unilim.fr>
//                - Jean-luc Charles   Arts et Metiers ParisTech, CNRS, I2M, Bordeaux France
//                  <jean-luc.charles@ensam.eu>
//                - Jeremie Girardot   Arts et Metiers ParisTech, CNRS, I2M, Bordeaux France
//                  <jeremie.girardot@ensam.eu>
//                - Cedric Hubert      LAMIH/UPHF, Valenciennes France
//                  <cedric.hubert@uphf.fr>
//                - Ivan Iordanoff     Arts et Metiers ParisTech, CNRS, I2M, Bordeaux France
//                  <ivan.iordanoff@ensam.eu>
//
// Copyright (C) 2008-2019 D. Andre, JL. Charles, J. Girardot, C. Hubert, I. Iordanoff
//
// This program 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.
//
// This program 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 this program.  If not, see <http://www.gnu.org/licenses/>.
 
 
 
#ifndef _GranOO_libGeom_Plane_hpp_
#define _GranOO_libGeom_Plane_hpp_
 
 
#include "GranOO3/Geom/Ray.hpp"
 
 
namespace GranOO3
{
  namespace Geom
  {
 
    class Plane;
    
    bool intersect(const Ray& ray, const Plane& plane, Geom::Point& intersection_pt);
    
    class Plane
    {
 
    public:
      // Constructors and destructor
      Plane (const Geom::Vector& n, double distance);
 
    private:
      Plane() = delete;
 
    public:
      const Geom::Vector& normal;
      const double        distance;
    };
 
 
    inline
    Plane::Plane (const Geom::Vector& n, double d)
      : normal(n),  distance(d) {
      AssertMsg(normal.is_unit() == true, "The plane's normal must be a unit a vector !"); 
    }
 
    inline
    bool intersect(const Ray& ray, const Plane& plane, Geom::Point& intersection_pt) {
      const double t  = -(ray.point.to_vector() * plane.normal + plane.distance) / (ray.direction * plane.normal);
      intersection_pt =  (ray.point.to_vector() + t*ray.direction).to_point();
      return t < 0.;
    }
 
 
  }
}
 
 
#endif