BSumSurf.h

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

class MGIfstream;
class MGOfstream;
#include <iosfwd>
#include "mg/MGCL.h"
#include "mg/CSisect_list.h"
#include "mg/Surface.h"

// MGBSumSurf.h


class MG_DLL_DECLR MGBSumSurf: public MGSurface{
 
public:


MGBSumSurf():MGSurface(),m_g1(0), m_g2(0), m_g12(0){;};

MGBSumSurf(
    MGSurface* g1,
    MGSurface* g2,
    MGSurface* g12
);

MGBSumSurf(
    const MGSurface& g1,
    const MGSurface& g2,
    const MGSurface& g12
);

MGBSumSurf(const MGBSumSurf& rhs);


~MGBSumSurf();
                                    

MGBSumSurf& operator=(const MGGel& gel2);
MGBSumSurf& operator=(const MGBSumSurf& gel2);

MGBSumSurf& operator+=(const MGVector& v);
MGBSumSurf& operator-=(const MGVector& v);
MGBSumSurf& operator*=(double scale);
MGBSumSurf& operator*=(const MGMatrix& mat);
MGBSumSurf& operator*=(const MGTransf& tr);

bool operator==(const MGBSumSurf& gel2)const;
bool operator==(const MGGel& gel2)const;
bool operator<(const MGBSumSurf& gel2)const;
bool operator<(const MGGel& gel2)const;
bool operator!=(const MGGel& gel2)const{return !(gel2==(*this));};
bool operator!=(const MGBSumSurf& gel2)const{return !(gel2==(*this));};


int bdim_u() const; 
int bdim_v() const; 

MGBox box_limitted(
    const MGBox& uvrange    
)const;

MGBSumSurf& change_dimension(
    int sdim,   
    int start1=0,
    int start2=0    
);

MGBSumSurf& change_range(
    int is_u,               
    double t1,              
    double t2               
);

MGBSumSurf* clone() const;

MGBSumSurf* copy_change_dimension(
    int sdim,   
    int start1=0,
    int start2=0 
)const;

MGVector eval(
    double u,
    double v 
    , int ndu=0     
    , int ndv=0     
) const;

MGVector eval(
    const MGPosition& uv    
    , int ndu           
    , int ndv           
) const{return eval(uv[0],uv[1],ndu,ndv);}

MGSurface& exchange_uv();

long identify_type() const{return MGBSUMSURF_TID;};

bool in_range(double u, double v) const;

MGCSisect_list isect(const MGCurve& curve)const;
MGCSisect_list isect(const MGStraight& sl)const{return isectSl(sl);};
MGCSisect_list isect(const MGRLBRep& curve)const;
MGCSisect_list isect(const MGEllipse& curve)const;
MGCSisect_list isect(const MGLBRep& curve)const;
MGCSisect_list isect(const MGSurfCurve& curve)const;
MGCSisect_list isect(const MGBSumCurve& curve)const;

MGSSisect_list isect(const MGSurface& srf2)const;
MGSSisect_list isect(const MGPlane& srf2)const;
MGSSisect_list isect(const MGSphere& srf2)const;
MGSSisect_list isect(const MGCylinder& srf2)const;
MGSSisect_list isect(const MGSBRep& srf2)const;
MGSSisect_list isect(const MGRSBRep& srf2)const;
MGSSisect_list isect(const MGBSumSurf& srf2)const;

int isect_order() const;

double knot_u(int i) const;

double knot_v(int i) const;

const MGKnotVector& knot_vector_u() const;
MGKnotVector& knot_vector_u();

const MGKnotVector& knot_vector_v() const;
MGKnotVector& knot_vector_v();

bool less_than(
    int i,  
    const MGPosition& uv1,
    const MGPosition& uv2
)const{return m_g1->less_than(i,uv1,uv2);};

void negate(
    int is_u
);

int order_u() const;

int order_v() const;


std::ostream& out(std::ostream& ostrm) const;

MGCurve* parameter_curve(
    int is_u    
    , double x  
)const;

double param_e_u() const;
double param_e_v() const;

MGBox param_range() const;

double param_s_u() const;
double param_s_v() const;

MGBSumSurf* part(
    const MGBox& bx,
    int multiple=0  
)const;

int perimeter_num() const{return 4;};

int sdim()const{return m_g1->sdim();};

void shrink_to_knot(
    const MGBox& uvbx,
    int multiple=0  
);

MGSURFACE_TYPE type() const{return MGSURFACE_BSUM;}

void ReadMembers(MGIfstream& buf);
    
void WriteMembers(MGOfstream& buf) const;

std::string whoami()const{return "BSumSurf";};


private:

    MGSurface* m_g1;
    MGSurface* m_g2;
    MGSurface* m_g12;

MGBox* compute_box() const;

int intersect_dnum_u() const;
int intersect_dnum_v() const;

MGCurve* isect_incr_pline(
    const MGPosition& uv,   
    int kdt,                
    double du, double dv,
    double& u,              
    double& v,              
    int incr=0          
)const;

MGSBRep* surf1D(const MGPlane& pl)const;

};
 // end of GEO group
#endif