plane3D.hpp

#ifndef PLANE_3D_HPP
#define PLANE_3D_HPP

#include "jmath/jblas.hpp"
#include <iostream>
#include "model3d/point3D.hpp"

namespace jafar {
  namespace model3d {
    
    class Plane3D
    {
    public:
      Plane3D();
      ~Plane3D();
      
      Plane3D& operator= (const Plane3D& srcPlane);
      
      void setPlane(Point3D& p1, Point3D& p2, Point3D& p3);
      
      void info();
      
      void clear();
      bool isBaseNormal(int& nType);
      void setBaseNormal(int nType);
      bool isVertical();
      
      inline double calculateCosDihedralAngle(Plane3D& plane3D_2)
      {
        return (u[0]*plane3D_2.u[0] + 
                u[1]*plane3D_2.u[1] + 
                u[2]*plane3D_2.u[2]);
      }
      
      
      // signed distance betweed the plane and the point
      inline double signedDistance(const Point3D& p) const
      {
        return (u[0]*p.x()+u[1]*p.y()+u[2]*p.z()+d0);
      }


      bool isCoplanar(Plane3D& plane3D_2, double dist, double cosDihedral);
      
      //data member
      double u[3];
      double d0;
    };
    
    
    
  }// namespace model3d
}// namespace jafar

#endif //PLANE_3D_HPP