dtk_pnt.hpp

Go to the documentation of this file.

#ifndef __DTK_PNT_HPP__
#define __DTK_PNT_HPP__
 
#include <fstream>
#include <math.h>
#include "def/define.h" 
#include <cmath>
#include <iosfwd>
 
class Dtk_pnt;
class Dtk_dir;
class Dtk_string;
class Dtk_status;
class CDtk_3dMat;
 
class Dtk_pnt
{
protected:
     double _Coords[ 3 ];
public:
     Dtk_pnt();
     ~Dtk_pnt() {}
     Dtk_pnt( const Dtk_pnt& s );
     Dtk_pnt( double a, double b, double c = 0.0 );
     Dtk_pnt( const double d[] );
     Dtk_pnt& operator = ( const Dtk_pnt& s );
     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 ]; }
 
     inline void GetXYZ( Dtk_Double64& outX, Dtk_Double64& outY, Dtk_Double64& outZ ) const
     {
          outX = _Coords[ 0 ];
          outY = _Coords[ 1 ];
          outZ = _Coords[ 2 ];
     }
 
     inline void SetXYZ( const Dtk_Double64& inX, const Dtk_Double64& inY, const Dtk_Double64& inZ )
     {
          _Coords[ 0 ] = inX;
          _Coords[ 1 ] = inY;
          _Coords[ 2 ] = inZ;
     }
 
     Dtk_pnt& operator*=( double );
     Dtk_pnt& operator/=( double );
     Dtk_pnt operator-();
     Dtk_pnt& operator+=( const Dtk_dir& );
     Dtk_pnt& operator-=( const Dtk_dir& );
     friend Dtk_pnt operator+( const Dtk_pnt&, const Dtk_dir& );  // pour barycentre
     friend Dtk_dir operator-( const Dtk_pnt&, const Dtk_dir& );
 
     inline Dtk_pnt operator* ( double a ) const { return Dtk_pnt( _Coords[ 0 ] * a, _Coords[ 1 ] * a, _Coords[ 2 ] * a ); }
     friend Dtk_pnt operator* ( double, const Dtk_pnt& );
     friend Dtk_pnt operator/ ( const Dtk_pnt&, double );
     friend Dtk_pnt operator/ ( double, const Dtk_pnt& );
 
     friend int operator==( const Dtk_pnt&, const Dtk_pnt& );
     friend int operator!=( const Dtk_pnt&, const Dtk_pnt& );
 
     int IsEqual( const Dtk_pnt& inPointToCompare, const double inTolerance ) const;
     //  Dtk_pnt operator=(double d[]); // inutile
 
     inline static void ProcessMiddlePoint( const Dtk_pnt& inPntA, const Dtk_pnt& inPntB, Dtk_pnt& outMiddlePnt )
     {
          for( Dtk_Size_t i = 0; i < 3; i++ )
               outMiddlePnt[ i ] = ( inPntA[ i ] + inPntB[ i ] ) / 2.0;
     }
 
     inline static Dtk_Double64 ProcessDistance( const Dtk_pnt& inPoint1, const Dtk_pnt& inPoint2 )
     {
          Dtk_Double64 dist = 0.0;
          Dtk_Size_t i;
 
          for( i = 0; i < 3; i++ )
               dist += ( inPoint1[ i ] - inPoint2[ i ] ) * ( inPoint1[ i ] - inPoint2[ i ] );
 
          return static_cast< Dtk_Double64 >( sqrt( dist ) );
     }
 
     /*   void SetInvalid();
          Dtk_bool IsInvalid()const;*/
 
     friend std::ostream& operator<<( std::ostream& o, const Dtk_pnt& d );
 
     friend double Argument( const Dtk_pnt& O, const Dtk_dir& x, const Dtk_dir& y, const Dtk_pnt& P );
     friend double Argument( const Dtk_pnt& O, const Dtk_dir& x, const Dtk_pnt& P );
     friend double Argument( const Dtk_pnt& P );
 
     double Solve( const Dtk_dir& col0, const Dtk_dir& col1, const Dtk_dir& col2, const Dtk_dir& col3 );
     double Solve( const CDtk_3dMat& matrix, const Dtk_dir& col3 );
     void Multiply( const CDtk_3dMat& matrix );
     inline void operator *= ( const CDtk_3dMat& matrix );
     Dtk_pnt Multiplied( const CDtk_3dMat& matrix ) const;
     inline Dtk_pnt operator * ( const CDtk_3dMat& matrix ) const;
     void PostMultiply( const CDtk_3dMat& matrix );
     Dtk_pnt PostMultiplied( const CDtk_3dMat& matrix ) const;
 
     void MultiplyByTransposed( const CDtk_3dMat& matrix );
     Dtk_pnt MultipliedByTransposed( const CDtk_3dMat& matrix ) const;
 
     double Project( const Dtk_pnt& org, const Dtk_dir& zdir, const Dtk_dir& zPrj );
     inline Dtk_pnt Projected( const Dtk_pnt& org, const Dtk_dir& zdir, const Dtk_dir& zPrj ) const;
     double Project( const Dtk_pnt& org, const Dtk_dir& zdir );
     inline Dtk_pnt Projected( const Dtk_pnt& org, const Dtk_dir& zdir ) const;
     double Project( const Dtk_dir& zdir, double d, const Dtk_dir& zPrj );
     inline Dtk_pnt Projected( const Dtk_dir& zdir, double d, const Dtk_dir& zPrj ) const;
     double Project( const Dtk_dir& zdir, double d );
     inline Dtk_pnt Projected( const Dtk_dir& zdir, double d ) const;
     void   ProjectOnLine( const Dtk_pnt& org, const Dtk_dir& dir );
     inline Dtk_pnt ProjectedOnLine( const Dtk_pnt& org, const Dtk_dir& dir ) const;
 
     void    Scale( const Dtk_pnt& p0, double s );
     Dtk_pnt Scaled( const Dtk_pnt& p0, double s ) const;
     void    Translate( const Dtk_dir& dir );
     Dtk_pnt Translated( const Dtk_dir& dir ) const;
     void Translate( const Dtk_pnt& p0, const Dtk_pnt& p1 );
     Dtk_pnt   Translated( const Dtk_pnt& p0, const Dtk_pnt& p1 ) const;
     void    Transform( const CDtk_3dMat& matrix );
     Dtk_pnt Transformed( const CDtk_3dMat& matrix ) const;
     void    Rotate( const Dtk_pnt& p0, const CDtk_3dMat& matrix );
     Dtk_pnt Rotated( const Dtk_pnt& p0, const CDtk_3dMat& matrix ) const;
     void    Rotate( const Dtk_dir& dir, double angle, const Dtk_pnt& center );
     Dtk_pnt   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_pnt Symetrised( const Dtk_pnt& opl, const Dtk_dir& zpl ) const;
     void    Affinity( const Dtk_pnt&, const Dtk_dir& xdir, const Dtk_dir& ydir, double f[ 3 ] );
     inline void  Affinity( double f[ 3 ] ) { _Coords[ 0 ] *= f[ 0 ]; _Coords[ 1 ] *= f[ 1 ];  _Coords[ 2 ] *= f[ 2 ]; }
     Dtk_pnt Affinitised( const Dtk_pnt& p, const Dtk_dir& xdir, const Dtk_dir& ydir, double f[ 3 ] ) const;
};
 
inline void Dtk_pnt::operator *= ( const CDtk_3dMat& matrix )
{
     Multiply( matrix );
}
 
inline Dtk_pnt Dtk_pnt::operator * ( const CDtk_3dMat& matrix ) const
{
     return Multiplied( matrix );
}
 
inline Dtk_pnt Dtk_pnt::Projected( const Dtk_pnt& org, const Dtk_dir& zdir, const Dtk_dir& zprj ) const
{
     Dtk_pnt tmp = *this;
     tmp.Project( org, zdir, zprj );
     return tmp;
}
 
inline Dtk_pnt Dtk_pnt::Projected( const Dtk_pnt& org, const Dtk_dir& zdir ) const
{
     Dtk_pnt tmp = *this;
     tmp.Project( org, zdir );
     return tmp;
}
 
inline Dtk_pnt Dtk_pnt::Projected( const Dtk_dir& zdir, double d, const Dtk_dir& zprj ) const
{
     Dtk_pnt tmp = *this;
     tmp.Project( zdir, d, zprj );
     return tmp;
}
 
inline Dtk_pnt Dtk_pnt::Projected( const Dtk_dir& zdir, double d ) const
{
     Dtk_pnt tmp = *this;
     tmp.Project( zdir, d );
     return tmp;
}
 
inline Dtk_pnt Dtk_pnt::ProjectedOnLine( const Dtk_pnt& org, const Dtk_dir& dir ) const
{
     Dtk_pnt tmp = *this;
     tmp.ProjectOnLine( org, dir );
     return tmp;
}
 
#endif