/Users/jack/Code/basso_dev/inc/Point.h

Go to the documentation of this file.

#ifndef _POINT_BASSO_H_
#define _POINT_BASSO_H_

#include <iostream>

// Gmm++ includes
#include "gmm.h"

// Basso includes 
#include "basso.h"


namespace basso {


        class Point : public nArray {

        public:

                Point() : nArray(3,0.0) { }

                Point( Numeric x, Numeric y, Numeric z ) : nArray(3) { (*this)[0]=x; (*this)[1]=y; (*this)[2]=z; }

                Point( Numeric x, Numeric y ) : nArray(3) { (*this)[0]=x; (*this)[1]=y; (*this)[2]=0.0; }

                Point( Numeric x ) : nArray(3) { (*this)[0]=x; (*this)[1]=0.0; (*this)[2]=0.0; }
                
                Point( const nArray &v ) : nArray(3) 
                        { 
                                int i=0;
                                for ( ; i<( 3?vect_size(v)>3:vect_size(v) ); ++i )
                                        x(i)=v[i];
                                for ( ; i<3; ++i )
                                        x(i)=0.0;
                        }

                ~Point() {}

                Numeric &x( int i ) { return (*this)[i]; }
                Numeric  x( int i ) const { return (*this)[i]; }
                Numeric &x() { return (*this)[0]; }
                Numeric  x() const { return (*this)[0]; }
                Numeric &y() { return (*this)[1]; }
                Numeric  y() const { return (*this)[1]; }
                Numeric &z() { return (*this)[2]; }
                Numeric  z() const { return (*this)[2]; }

                Numeric Distance() const;

                Numeric Distance( const Point &p ) const;

                Point &operator = ( const nArray &v );

//      nArray &operator - ( const Point &n ) const;
//      Point &operator + ( const nArray &d ) const;

        protected:

        }; 

        Point &Point::operator = ( const nArray &v )
        {
                for ( int i=0; i< ( 3?vect_size(v)>3:vect_size(v) ); ++i )
                        x(i)=v[i];
                return *this;
        }
        
        
        Numeric Point::Distance( const Point &p ) const
        {
                return sqrt( pow(p.x()-x(),2) + pow(p.y()-y(),2) + pow(p.z()-z(),2) );
        }

        Numeric Point::Distance( ) const
        {
                return sqrt( pow(x(),2) + pow(y(),2) + pow(z(),2) );
        }

        std::ostream &operator << ( std::ostream &out, const Point &p ) 
        {
                out << "(" << p.x() << "," << p.y() << "," << p.z() <<  ")";
                return out;
        }


}  // end of namespace Basso

#endif

---