CellNB.h

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

#include <vector>
#include "topo/CellBase.h"
#include "mg/Position.h"

//
//Defines MGCellNB Class.

class MGBox;
class MGComplex;
class MGGeometry;
class MGCellBase;
class MGBVertex;
class MGCellMap;


class MG_DLL_DECLR MGCellNB:public MGCellBase{

public:

typedef std::vector<MGCellBase*>::iterator partnerItr;
typedef std::vector<const MGCellBase*>::const_iterator const_partnerItr;
typedef std::vector<MGCellBase*>::reverse_iterator partnerRItr;
typedef std::vector<const MGCellBase*>::const_reverse_iterator const_partnerRItr;

enum CELL_KIND{
    UNKNOWN=0,
    PCELL=1,
    BCELL=2
};


MGCellNB();

MGCellNB(const MGCellNB& cell);

MGCellNB(const MGGeometry& geo);
explicit MGCellNB(
    MGGeometry* geo
);

MGCellNB(
    MGGeometry* geo,
    MGCellNB* binder    
);

virtual ~MGCellNB();


virtual MGCellNB& operator=(const MGCellNB& gel2);

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


void add_partner(const MGCellBase& partner);

const MGBox& box() const=0;

MGPosition center() const;

virtual MGPosition center_param() const=0;

virtual MGCellNB* clone() const=0;

virtual MGCellNB* clone_without_boundaries() const=0;

virtual MGCellNB* clone_binder(const MGCellBase& c) const=0;

virtual MGUnit_vector direction() const;

virtual void draw3DVertex(mgVBO& vbo)const=0;

const MGGeometry* extent() const{ return m_extent;};
MGGeometry* extent() {return m_extent;};

MGGeometry* free_extent();

MGComplex* free_from_parent();

virtual void free_neighbourhood(int i, int j=0)=0;

void free_partner(const MGCellBase* cellin) const;

virtual long identify_type()const=0;

bool is_bcell() const{return m_partners.size()>=1;};

virtual MGCellNB* make_binder() const=0;

virtual int manifold_dimension() const=0;

const MGCellBase* member_partner(int i)const{return m_partners[i];};

const std::vector<const MGCellBase*>& member_partners()const{return m_partners;};
std::vector<const MGCellBase*>& member_partners(){return m_partners;};

virtual void negate();

int number_of_partner_members() const{return int(m_partners.size());};

virtual std::vector<const MGCellNB*> neighbours() const=0;

virtual double parameter_error()const=0;

const MGComplex* parent_complex() const{return m_parent_complex;};
MGComplex* parent_complex() {return m_parent_complex;};

virtual void set_extent(MGGeometry* extent=0);

const MGCellNB* star() const;
MGCellNB* star();

virtual std::ostream& out(std::ostream&) const;
virtual std::string whoami()const{return "CellNB";};

protected:

    mutable MGComplex* m_parent_complex;
    MGGeometry* m_extent;   
    mutable std::vector<const MGCellBase*> m_partners;

bool bn_binder_tr_necessary()const{ return !parent_complex();}

virtual MGCellNB* clone(MGCellMap& cmap) const;

virtual void compute_box() const=0;

bool is_less_than(const MGCellNB& cell2)const;

virtual void ReadMembers(MGIfstream& buf);


virtual void WriteMembers(MGOfstream& buf) const;

MGCellNB& set_cellnb(const MGCellNB& cell2);

private:

virtual void bn_binder_tr(const MGVector& v)=0;
virtual void bn_binder_tr(double s)=0;
virtual void bn_binder_tr(const MGMatrix& mat)=0;
virtual void bn_binder_tr(const MGTransf& tr)=0;

virtual void set_box_as_null() const=0;

virtual void copy_all_boundaries(const MGCellBase& cell2)=0;

virtual void copy_all_boundaries(const MGCellBase& cell2, MGCellMap& cmap)=0;

virtual void copy_box(const MGCellBase& cell2) const=0;

virtual void copy_perror(const MGCellBase& cell2) const=0;

void free_partner(const_partnerItr itr) const;

virtual void get_all_boundary_binders(std::vector<MGCellNB*>& cvec) const=0;

virtual MGGeometry* make_binder_extent() const=0;

virtual void make_extent() const=0;

virtual void merge_bcell(MGCellNB* bcell2);

virtual void negate_boundary()=0;

friend class MGCellBase;
friend class MGComplex;
friend class MGEdge;

};
 // end of TOPO group
#endif