SPointSeq.h

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

#include "mg/MGCL.h"

// MGSPointSeq.h
//

class MGBox;
class MGNDDArray;
class MGVector;
class MGPosition;
class MGMatrix;
class MGTransf;
class MGBPointSeq;
class MGPlane;
class MGStraight;
class MGIfstream;
class MGOfstream;


class MG_DLL_DECLR MGSPointSeq {

public:

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

MG_DLL_DECLR friend std::ostream& operator<< (std::ostream& out, const MGSPointSeq& sp);


explicit MGSPointSeq(
    int sizeu=0,    
    int sizev=0,    
    int dim=0   
);

MGSPointSeq(
    int dim,                        
    const MGSPointSeq& old_brep,    
    int start1=0,
    int start2=0    
);

MGSPointSeq(const MGSPointSeq& rhs);

~MGSPointSeq(){if(m_spoint) delete[] m_spoint;};

                                          
MGSPointSeq& operator=(const MGSPointSeq&);

double& operator()(int i, int j, int k); 
double operator()(int i, int j, int k) const{return ref(i,j,k);};

MGVector operator()(int i, int j) const;

MGSPointSeq operator+ (const MGVector& ) const;

MGSPointSeq& operator+= (const MGVector& );

MGSPointSeq operator- (const MGVector& ) const;

MGSPointSeq& operator-= (const MGVector& );

MGSPointSeq operator* (double) const;

MGSPointSeq& operator*= (double);

MGSPointSeq operator* (const MGMatrix& ) const;

MGSPointSeq& operator*= (const MGMatrix& );

MGSPointSeq operator* (const MGTransf& ) const;

MGSPointSeq& operator*= (const MGTransf& );

MGSPointSeq operator+ (const MGSPointSeq& sp2) const;
MGSPointSeq& operator+= (const MGSPointSeq& sp2);
MGSPointSeq operator- (const MGSPointSeq& sp2) const;
MGSPointSeq& operator-= (const MGSPointSeq& sp2);

bool operator== (const MGSPointSeq& ) const;

bool operator!= (const MGSPointSeq& brep) const{return !(operator== (brep));}


int average_plane(
    MGPosition& center  
    , MGPlane& plane    
    , MGStraight& line  
    , double& deviation 
)const;                 

MGBox box()const;

const double* data(int i=0, int j=0, int k=0) const;

double* data(int i=0, int j=0, int k=0);

MGSPointSeq& homogeneous_transform(double);
MGSPointSeq& homogeneous_transform(const MGVector&);
MGSPointSeq& homogeneous_transform(const MGMatrix&);
MGSPointSeq& homogeneous_transform(const MGTransf&);

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

void length(int& lengthu, int& lengthv) const
{lengthu=m_lengthu;lengthv=m_lengthv;}

int length_u() const{return m_lengthu;}

int length_v() const{return m_lengthv;}

void make_data_point(MGNDDArray& utau, MGNDDArray& vtau) const;

MGSPointSeq non_homogeneous() const;

double ref(int i, int j, int k=0) const;

void reshape(
      int sizeu     
    , int suzev     
    , int startu=0  
    , int startv=0  
);

void resize(int lenu, int lenv,  int dim);

void reverse(
    int is_u                
);

int sdim() const {return m_sdim;}

void set_null();

void set_length(int lengthu,
                int lengthv );

void capacity(int& capau, int& capav) const{capau=m_capacityu; capav=m_capacityv;}

int capacity_u() const {return m_capacityu;}

int capacity_v() const {return m_capacityv;}

void store_at(
    int i,
    int j,
    const MGVector& vctr,   
    int to=0,   
    int from=0  
);

void store_at(
    int i,
    int j,
    const MGVector& vctr,
    int to, 
    int from,
    int len 
);

void store_at(
    int i,
    int j,
    const double* data, 
    int to=0        
);

void store_BP_along_u_at(
    int j,  
    const MGBPointSeq& bp,  
    int to=0,   
    int from=0  
);

void store_BP_along_v_at(
    int i,  
    const MGBPointSeq& bp,  
    int to=0,   
    int from=0  
);

int dump_size() const;

int dump(MGOfstream& ) const;

int restore(MGIfstream& );

private:

    int m_capacityu;    
    int m_capacityv;    
    int m_sdim;     
    int m_lengthu;  
    int m_lengthv;  
    double* m_spoint;   

MGBox* compute_box()const;

friend class MGSBRep;
friend class MGRSBRep;

};
 // end of BASE group
#endif