Vector.h

/********************************************************************/
/* Copyright (c) 2015 DG Technologies Inc. and Yuzi Mizuno          */
/* All rights reserved.                                             */
/********************************************************************/
#ifndef _MGVector_HH_
#define _MGVector_HH_

#include <stddef.h>
#include <vector>
#include "mg/MGCL.h"

// MGVector.h
// Header for MGVector.

//Forward Declaration
class MGUnit_vector;
class MGPosition;
class MGIfstream;
class MGOfstream;
class MGIgesOfstream;

class MG_DLL_DECLR MGVector {

public:

MG_DLL_DECLR friend MGVector operator+(const MGVector& vec1,const MGVector& vec2);

MG_DLL_DECLR friend MGVector operator-(const MGVector& vec1,const MGVector& vec2);

MG_DLL_DECLR friend double operator%(const MGVector& vec1,const MGVector& vec2);

MG_DLL_DECLR friend MGVector operator*(const MGVector& vec1,const MGVector& vec2);

MG_DLL_DECLR friend MGVector operator*(const MGVector& vec1,double scale);

MG_DLL_DECLR friend MGVector operator* (double, const MGVector&);

MG_DLL_DECLR friend MGVector operator/(const MGVector& vec1,double scale);

inline friend
bool operator<(const MGVector& v1,const MGVector& v2){return v1.len()<v2.len();};
inline friend
bool operator<=(const MGVector& v1,const MGVector& v2){return v1.len()<=v2.len();};
inline friend
bool operator>(const MGVector& v1,const MGVector& v2){return v1.len()>v2.len();};
inline friend
bool operator>=(const MGVector& v1,const MGVector& v2){return v1.len()>=v2.len();};

MG_DLL_DECLR friend bool operator==(const MGVector& v1,const MGVector& v2);

inline friend bool operator!=(const MGVector& v1,const MGVector& v2){return !(v1==v2);}

MG_DLL_DECLR friend std::ostream& operator<< (std::ostream&, const MGVector&);

MG_DLL_DECLR friend double MGDeterminant(
        const MGVector& v1, const MGVector& v2, const MGVector& v3
);


explicit MGVector(int sdim=0);

MGVector(double x, double y);

MGVector(double x, double y , double z);

MGVector(double x, double y, double z, double w);

MGVector(int sdim, double v);

MGVector(int sdim, const double* v);

MGVector(const MGPosition&);

MGVector(
    const MGVector& dvec,
    const MGVector& svec
);

MGVector(
    int sdim,           
    const MGVector& vec2,
    int start1=0,       
    int start2=0        
);

MGVector(const std::vector<double>& darrays);


MGVector ( const MGVector& );

~MGVector(){if(m_sdim>3) delete[] m_element;}


MGVector& operator =(const MGVector &);

double operator() (int i) const{return ref(i);}  

double operator[] (int i) const{return ref(i);}  

double& operator()(int i);

MGVector& operator=(const double*);

MGVector & operator+= (const MGVector&);

MGVector operator- () const;

MGVector & operator -= ( const MGVector & );

MGVector& operator*= (double scale);

MGVector& operator*= (const MGVector& vec2);

MGVector& operator/= (double scale);


double angle(const MGVector&) const;

double anglepai(const MGVector& v2)const{return angle(v2);};

double angle2pai(const MGVector& v2, const MGVector& N)const;

double cangle(const MGVector&) const;

MGVector& clear(double init=0.0);

const double* data()const{return m_element;};
double* data(){return m_element;};

MGVector interpolate(double t, const MGVector& vec2) const;

MGVector interpolate_by_rotate(
    double t, const MGVector& vec2,
    double* ratio=0
) const;

bool is_collinear(
    const MGVector& v2
)const{ return (*this).parallel(v2);}

bool is_collinear(
    const MGVector& v2,
    const MGVector& v3
)const;

bool is_null()const{return m_sdim==0;}

bool is_unit_vector() const;

bool is_zero_vector() const;

double len() const;

void negate(){operator*=(-1.);};

MGUnit_vector normalize() const;

bool orthogonal(const MGVector& ) const;

int out_to_IGES(
    MGIgesOfstream& igesfile,
    int SubordinateEntitySwitch=0
)const;

MGVector orthogonize(const MGVector& vec2)const;

bool parallel(const MGVector& ) const;

MGVector project(const MGVector& v2) const;

double ref(int i) const{ 
    if(i<sdim()) return m_element[i];
    else         return 0.;
}

void resize(int new_sdim);

double sangle(const MGVector& ) const;

int sdim() const { return m_sdim; };

void set_null();

void set_unit();

void store_at(
    int i,      
    const MGVector& vec2,
    int j=0     
);

void store_at(
    int i,      
    const MGVector& vec2,
    int j,      
    int len     
);

void swap(int i, int j);

int dump_size() const;
int dump(MGOfstream& ) const;

int restore(MGIfstream& );


protected:
    int m_sdim;
    double* m_element;
    double m_data[3];   
    mutable double m_length;

    double& set(int i) {return m_element[i];}

friend class MGLBRep;
friend class MGSBRep;

};

namespace MGCL{


MG_DLL_DECLR MGVector project(const MGVector& V1, const MGVector& V2);

};
 // end of BASE group
#endif