PPRep.h

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

#include "mg/NDDArray.h"
#include "mg/Vector.h"

// MGPPRep.h
//

// Forward Declaration
class MGLBRep;
class MGIfstream;
class MGOfstream;


class MG_DLL_DECLR MGPPRep {
 
public:

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


MGPPRep();

MGPPRep(int order, int nbreak, int sdim);

MGPPRep(int order, int sdim, const MGNDDArray& tau);

MGPPRep(const MGLBRep& lbrep);

MGPPRep(int order, const MGPPRep& pprep);

MGPPRep(const MGPPRep& rhs);


~MGPPRep(){if(m_coef) delete[] m_coef;};
                                    

MGPPRep& operator=(const MGPPRep&); 

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

MGVector operator()(int i, int j)const{return coef(i,j);};


const MGNDDArray& break_point()const{return m_break_point;}

double break_point(int i)const{return m_break_point(i);}

double& break_point(int i){return m_break_point(i);}

const double* break_point_data()const{return m_break_point.data();}     

double coef(int i, int j, int k)const;

double& coef(int i, int j, int k);

MGVector coef(int i, int j)const;

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

MGVector eval(      
    double t,   
    int n=0 
)const;

MGVector eval_i(    
    int i,      
    double t,       
    int n       
)const;     

int nbreak()const{return m_nbreak;}

MGPPRep& normalize();

int order()const{return m_order;}

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

void reshape(int nbreak); 

void resize(
    int order,  
    int nbreak, 
    int dim 
);

int sdim()const{return m_sdim;}

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

int dump_size()const;

int dump(MGOfstream& )const;

int restore(MGIfstream& );


private:

    int m_order;        
    int m_nbreak;   
    int m_sdim;     
    MGNDDArray m_break_point;
    double* m_coef;     

};
 // end of GEORelated group
#endif