dtk_dir.hpp

Go to the documentation of this file.

#ifndef __DTK_DIR_HPP__
#define __DTK_DIR_HPP__
#include "util/BaseTypes/dtk_pnt.hpp"
#include <def/define.h>
#include <iosfwd>
 
class CDtk_3dMat;
class Dtk_pnt;
class Dtk_dir:public Dtk_pnt
{
public:
     Dtk_dir();
     ~Dtk_dir() {}
     Dtk_dir( const Dtk_dir& s );
     Dtk_dir( const Dtk_pnt& pnt );
     Dtk_dir( double a, double b, double c = 0.0 );
     Dtk_dir( const Dtk_pnt&, const Dtk_pnt& );
     Dtk_dir& operator =( const Dtk_dir& );
     Dtk_dir operator+() const;
     Dtk_dir operator-() const;
     Dtk_dir& operator+=( const Dtk_dir& );
     Dtk_dir& operator-=( const Dtk_dir& );
     Dtk_dir& operator+=( const Dtk_pnt& );
     Dtk_dir& operator-=( const Dtk_pnt& );
     Dtk_dir& operator*=( double );
     Dtk_dir& operator/=( double );
     double norm() const;
     int normalize();
     int rotate( const Dtk_dir& Z, double angle );
     //deprecated
     SetAsDeprecated( "Undefined", "Use ComputeAngleByQuadrant" )
          double angle_with_Dtk_dir( const Dtk_dir& V2, const Dtk_dir& ref_Z );
     Dtk_Double64 ComputeAngleByQuadrant( const Dtk_dir& V2, const Dtk_dir& ref_Z ) const;
 
     void reverse();
 
     inline Dtk_dir operator* ( double a ) const { return Dtk_dir( _Coords[ 0 ] * a, _Coords[ 1 ] * a, _Coords[ 2 ] * a ); }
     friend Dtk_dir operator* ( double, const Dtk_dir& );
     friend Dtk_dir operator/ ( const Dtk_dir&, double );
     friend Dtk_dir operator/ ( double, const Dtk_dir& );
 
     inline Dtk_dir operator+( const Dtk_dir& v ) { return Dtk_dir( _Coords[ 0 ] + v._Coords[ 0 ], _Coords[ 1 ] + v._Coords[ 1 ], _Coords[ 2 ] + v._Coords[ 2 ] ); }
     inline Dtk_dir operator-( const Dtk_dir& v ) { return Dtk_dir( _Coords[ 0 ] - v._Coords[ 0 ], _Coords[ 1 ] - v._Coords[ 1 ], _Coords[ 2 ] - v._Coords[ 2 ] ); }
 
     friend double operator* ( const Dtk_dir&, const Dtk_dir& ); // scalar product - produit scalaire
     friend Dtk_dir operator/ ( const Dtk_dir&, const Dtk_dir& );  // vector product - produit vectoriel
     friend Dtk_dir operator^ ( const Dtk_dir&, const Dtk_dir& );  // vector product - produit vectoriel
     friend double Pdtk_DotProduct( const Dtk_dir&, const Dtk_dir& ); // 
     friend Dtk_dir Pdtk_CrossProduct( const Dtk_dir&, const Dtk_dir& );  // 
 
     friend Dtk_dir Orthogonal( const Dtk_dir& ); // trouve un vecteur orthogonal
     friend double norm( const Dtk_dir& ); // norme
     friend double normC( const Dtk_dir& ); // norme carree
 
     friend std::ostream& operator<<( std::ostream& o, const Dtk_dir& d );
 
     // gestion vecteurs 2D.
     friend Dtk_pnt intersect2d( const Dtk_pnt& p1, const Dtk_dir& d1, const Dtk_pnt& p2, const Dtk_dir& d2 );
     friend Dtk_pnt orthoproj2d( const Dtk_pnt& p1, const Dtk_dir& d1, const Dtk_pnt& p2 );
     friend int onleft( const Dtk_dir&, const Dtk_dir& );
     inline double& operator[] ( Dtk_Size_t i )
     {
          Pdtk_Assert( i < 3 );
          return _Coords[ i ];
     }
     inline const double& operator[] ( Dtk_Size_t i ) const
     {
          Pdtk_Assert( i < 3 );
          return _Coords[ i ];
     }
 
     inline double & x() { return _Coords[ 0 ]; }
 
     inline const double & x() const { return _Coords[ 0 ]; }
 
     inline double & y() { return _Coords[ 1 ]; }
 
     inline const double & y() const { return _Coords[ 1 ]; }
 
     inline double & z() { return _Coords[ 2 ]; }
 
     inline const double & z() const { return _Coords[ 2 ]; }
 
 
     void Cross( const Dtk_dir& right );
     Dtk_dir Crossed( const Dtk_dir& right ) const;
     inline void operator ^= ( const Dtk_dir& right );
     void CrossCross( const Dtk_dir& c1, const Dtk_dir& c2 );
     Dtk_dir CrossCrossed( const Dtk_dir& c1, const Dtk_dir& c2 ) const;
 
     double SquareModulus();
     Dtk_dir GetNormal() const;
     double Normalize();
     Dtk_dir Normalized() const;
     Dtk_dir reversed() const;
 
     inline double Dot( const Dtk_dir& other ) const { return ( *this * other ); }
     inline double operator & ( const Dtk_dir& other ) const { return ( *this * other ); }
 
     double DotCross( const Dtk_dir& c1, const Dtk_dir& c2 ) const;
 
     void Multiply( const CDtk_3dMat& matrix );
     inline void operator *= ( const CDtk_3dMat& matrix );
     Dtk_dir Multiplied( const CDtk_3dMat& matrix ) const;
     inline Dtk_dir operator * ( const CDtk_3dMat& matrix ) const;
     void MultiplyNormalize( const CDtk_3dMat& matrix );
     Dtk_dir MultipliedNormalize( const CDtk_3dMat& matrix ) const;
     void PostMultiply( const CDtk_3dMat& matrix );
     Dtk_dir PostMultiplied( const CDtk_3dMat& matrix ) const;
     void MultiplyByTransposed( const CDtk_3dMat& matrix );
     Dtk_dir MultipliedByTransposed( const CDtk_3dMat& matrix ) const;
     double AngleXY( double tolerance = 0.00000001 ) const;
     double Angle( const Dtk_dir& vRef, const Dtk_dir& vAxe, double tolerance = 0.00000001 ) const;
 
     void    Scale( const Dtk_pnt& p0, double s );
     Dtk_dir Scaled( const Dtk_pnt& p0, double s ) const;
     void    Translate( const Dtk_dir& dir );
     Dtk_dir Translated( const Dtk_dir& dir ) const;
     void Translate( const Dtk_pnt& p0, const Dtk_pnt& p1 );
     Dtk_dir   Translated( const Dtk_pnt& p0, const Dtk_pnt& p1 ) const;
     void    Transform( const CDtk_3dMat& matrix );
     Dtk_dir Transformed( const CDtk_3dMat& matrix ) const;
     void    Rotate( const Dtk_pnt& p0, const CDtk_3dMat& matrix );
     Dtk_dir Rotated( const Dtk_pnt& p0, const CDtk_3dMat& matrix ) const;
     void    Rotate( const Dtk_dir& dir, double angle, const Dtk_pnt& center );
     Dtk_dir   Rotated( const Dtk_dir& dir, double angle, const Dtk_pnt& center ) const;
     void    Symetry( const Dtk_dir& zpl );
     void    Symetry( const Dtk_pnt& opl, const Dtk_dir& zpl );
     Dtk_dir Symetrised( const Dtk_pnt& opl, const Dtk_dir& zpl ) const;
     void    Affinity( const Dtk_pnt&, const Dtk_dir& xdir, const Dtk_dir& ydir, double[ 3 ] );
     inline void  Affinity( double f[ 3 ] ) { _Coords[ 0 ] *= f[ 0 ]; _Coords[ 1 ] *= f[ 1 ];  _Coords[ 2 ] *= f[ 2 ]; }
     Dtk_dir Affinitised( const Dtk_pnt& p, const Dtk_dir& xdir, const Dtk_dir& ydir, double f[ 3 ] ) const;
 
 
     Dtk_dir UnitizedDerivative( const Dtk_dir& thisPrime ) const;
     Dtk_dir UnitizedSecondDerivative( const Dtk_dir& thisPrime, const Dtk_dir& thisSecond ) const;
     Dtk_dir DerivativeOfCrossProduct( const Dtk_dir& other, const Dtk_dir& thisPrime, const Dtk_dir& otherPrime ) const;
     Dtk_dir SecondDerivativeOfCrossProduct( const Dtk_dir& other, const Dtk_dir& thisPrime, const Dtk_dir& otherPrime, const Dtk_dir& thisSecond, const Dtk_dir& otherSecond ) const;
     double DerivativeOfNorm( const Dtk_dir& thisPrime ) const;
     double SecondDerivativeOfNorm( const Dtk_dir& thisPrime, const Dtk_dir& thisSecond ) const;
};
 
 
inline void Dtk_dir::operator *= ( const CDtk_3dMat& matrix )
{
     Multiply( matrix );
}
 
inline Dtk_dir Dtk_dir::operator * ( const CDtk_3dMat& matrix ) const
{
     return Multiplied( matrix );
}
 
inline void Dtk_dir::operator ^= ( const Dtk_dir& right )
{
     Cross( right );
}
 
 
#endif