vector.h

basic mathematic classes used for math programming

#ifndef VECTOR_H
#define VECTOR_H
#include "mathclass.h"
class unit_vector;
class matrix;
class transf;

class vector
{
  private:
    m_real p[3];

    // negation
    friend vector      operator-( vector const& );
    friend unit_vector operator-( unit_vector const& );

    // addtion
    friend position& operator+=( position&, vector const& );
    friend vector&   operator+=( vector&, vector const& );
    friend vector    operator+( vector const&, vector const& );
    friend position  operator+( position const&, vector const& );
    friend position  operator+( vector const&, position const& );

    // subtraction
    friend vector    operator-( vector const&, vector const& );
    friend position  operator-( position const&, vector const& );
    friend position& operator-=( position&, vector const& );
    friend vector&   operator-=( vector&, vector const& );

    // dot product
    friend m_real    operator%( vector const&, vector const& );
    friend m_real    operator%( position const&, vector const& );
    friend m_real    operator%( vector const&, position const& );

    // cross product
    friend vector    operator*( vector const&, vector const& );
    friend position  operator*( position const&, unit_vector const& );
    friend position  operator*( unit_vector const&, position const& );

    // scalar Multiplication
    friend vector    operator*( vector const&, m_real );
    friend vector    operator*( m_real, vector const& );
    friend vector&   operator*=( vector&, m_real );

    // scalar Division
    friend vector    operator/( vector const&, m_real );
    friend m_real    operator/( vector const&, vector const& );
    friend vector&   operator/=( vector&, m_real );

    // matrix Multiplication
    friend vector    operator*( vector const&, matrix const& );
    friend vector    operator*( matrix const&, vector const& );

    friend vector&   operator*=( vector&, transf const& );
    friend vector    operator*( vector const&, transf const& );

    // functions
    friend m_real       len( vector const& );
    friend unit_vector  normalize( vector const& );

	friend vector       interpolate( m_real, vector const&, vector const& );
	friend vector       slerp_fixangle_interpolate(m_real, vector const&, vector const& );
	friend vector       catmullrom_interpolate(m_real, vector const&, vector const&, vector const&, vector const&);
	friend vector       bezier_interpolate(m_real, vector const&, vector const&, vector const&, vector const&, int);
	friend vector       bezier_interpolatep(m_real, vector const&, vector const&, vector const&, vector const& , int );

    friend m_real       angle( vector const&, vector const& );

    friend position     vector2position( vector const& );
    friend vector       position2vector( position const& );

    // stream
    friend ostream& operator<<( ostream&, vector const& );
    friend istream& operator>>( istream&, vector& );
    friend ostream& operator<<( ostream&, unit_vector const& );
    friend istream& operator>>( istream&, unit_vector& );

  public:
    // constructors
    vector() {p[0] = 0;p[1] = 0;p[2] = 0;}
    vector( m_real x, m_real y, m_real z ) { p[0]=x; p[1]=y; p[2]=z; }
    vector( m_real a[3] ) { p[0]=a[0]; p[1]=a[1]; p[2]=a[2]; }
	vector(const vector& a);
	vector& operator=(const vector& a);
	
	void   clear();

    // inquiry functions
    m_real& operator[](int i) { return p[i]; }

    m_real x() const { return p[0]; };
    m_real y() const { return p[1]; };
    m_real z() const { return p[2]; };
    void   getValue( m_real d[3] ) { d[0]=p[0]; d[1]=p[1]; d[2]=p[2]; }
    void   setValue( m_real d[3] ) { p[0]=d[0]; p[1]=d[1]; p[2]=d[2]; }
	m_real getValue( int n ) const { return p[n]; }
	vector setValue( m_real x, m_real y, m_real z )
								   { p[0]=x, p[1]=y, p[2]=z; return *this; }
	m_real setValue( int n, m_real x )
								   { return p[n]=x; }



	void normalizself();

	m_real length() const;
	m_real length2() const; // return the square of the length
	matrix cross() const;
   
    // change functions
    void set_x( m_real x ) { p[0]=x; };
    void set_y( m_real x ) { p[1]=x; };
    void set_z( m_real x ) { p[2]=x; };
};

#endif