BaseCell.h

#ifndef __BASECELL__
#define __BASECELL__

#include <Edge.h>
#include <Joint.h>
#include <RefDesc.h>
#include <fstream>

class TBaseCell
{
  protected:
    TRefDesc *RefDesc;

    TJoint **Joints;

    int ClipBoard;

    int Reference_ID;
    
    int Bd_Part;

    int CellIndex;

    int GlobalCellNo;
    
    int *normalOrientation;
    
    int region;

    int LayerCell;
    
#ifdef __3D__

    TEdge **Edges;
#endif

#ifdef  _MPI

  int  ClipBoard_Par; 
    
    int SubDomainNumber;

    int SubDomainLocalCellNo;

    bool OwnCell;

    bool HaloCell;

    bool DependentCell;
    
    bool SubDomainInterfaceCell;

    bool CrossEdgeCell;
    
    bool CrossVertexCell;

    int N_NeibProcesses;

    int *NeibProcessesIds;
 
#endif

    int LocalCellNo;   
    

  public:
    // Constructor
    TBaseCell(TRefDesc *refdesc);

    // Destructor
    virtual ~TBaseCell();

    // Methods
    int SetRefDesc(TRefDesc *newrefdesc)
    {
      if (newrefdesc->GetShapeDesc() == GetShapeDesc())
      {
        RefDesc = newrefdesc;
        return 0;
      }
      else
        if (newrefdesc->GetShapeDesc()->GetN_Vertices() == 
              GetShapeDesc()->GetN_Vertices())
        {
          RefDesc = newrefdesc;
          return 0;
        }
        else
          return -1;
    }

    TRefDesc *GetRefDesc() const
    { return RefDesc; }
    TShapeDesc *GetShapeDesc() const
    { return RefDesc->GetShapeDesc(); }
    Shapes GetType() const
    { return RefDesc->GetShapeDesc()->GetType(); }

    Refinements GetEdgeRef(int i) const
    { return RefDesc->GetEdgeRef(i); }

#ifdef __3D__

    Refinements GetFaceRef(int i)
    { return RefDesc->GetFaceRef(i); }

    int SetEdge(int E_i, TEdge *E)
    {
      Edges[E_i] = E;
      return 0;
    }

    TEdge *GetEdge(int E_i)
    { return Edges[E_i]; }

#endif

    virtual int SetVertex(int Vert_i, TVertex *Vert) = 0;
    virtual TVertex *GetVertex(int Vert_i) = 0;

    int SetJoint(int J_i, TJoint *J)
    {
      Joints[J_i] = J;
      return 0;
    }

    TJoint *GetJoint(int J_i)
    { return Joints[J_i]; }

    int GetN_Vertices()
    { return RefDesc->GetN_OrigVertices(); }
    int GetN_Edges()
    {  return RefDesc->GetN_OrigEdges(); }
    int GetN_Joints()
    {  return RefDesc->GetShapeDesc()->GetN_Joints(); }

    #ifdef __3D__

      int GetN_Faces()
      { return RefDesc->GetN_OrigFaces(); }
    #endif

    virtual int GetN_Children() = 0;
    virtual int GetN_Parents() = 0;

    virtual TBaseCell *GetChild(int C_i) =  0;
    virtual TBaseCell *GetParent() =  0;
    virtual int SetParent(TBaseCell *parent) =  0;
    virtual int GetChildNumber(TBaseCell *Me) =  0;

    virtual int Draw(std::ofstream &dat, double scale, double StartX,
                   double StartY) = 0;
    virtual int PS(std::ofstream &dat, double scale, double StartX,
                   double StartY) = 0;
    virtual int MD_raw(std::ofstream &dat) = 0;

    virtual int Refine(int RefLevel) = 0;

    #ifdef __MORTAR__

      virtual int RefineMortar(int RefLevel) = 0;
    #endif

    virtual int Derefine() = 0;

    virtual int RefDeref() = 0;
    virtual int Gen1RegMarks() = 0;

    virtual int Gen1RegGrid() = 0;
    virtual int Ref1Reg(int LocJointNum, TBaseCell *&RefCell) = 0;
    virtual int Check1Reg() = 0;
    virtual int SetNoRefinement() = 0;
    virtual int SetRegRefine() = 0;
    virtual int Set1Refine(int i)= 0;    
    virtual int IsToRefine() = 0;
    virtual int ExistChildren() = 0;
    virtual int MakeConfClosure() = 0;

    #ifdef __2D__

      virtual int LineMidXY(int J_i, int P_j, double &X, double &Y) = 0;
      virtual int LineMidT(int J_i, int SJ_j, double &T_0, double &T_1) = 0;
    #else
    #endif

    void SetClipBoard(int value)
    { ClipBoard=value; }
    int GetClipBoard()
    { return ClipBoard; }

    virtual double GetDiameter() = 0;

    virtual double GetShortestEdge() = 0;

    virtual double GetLengthWithReferenceMap() = 0;

    virtual double GetMeasure() = 0;
    
    double Get_hK(int cell_measure);

    virtual bool PointInCell(double X, double Y) = 0;

#ifdef __3D__
    virtual bool PointInCell(double X, double Y, double Z) = 0;

     // added 25.04.2010 for fixing refinement problem
     void CorrectBoundaryVertices(TVertex **NewVertices, TJoint **NewJoints);
#endif

    virtual int GetGeoLevel() = 0;

    virtual int GetSubGridID() = 0;

    void SetReference_ID(int val)
    { Reference_ID = val;}
    
    int GetReference_ID() const
    {return Reference_ID;}
    
    void SetBd_Part(int val)
       {Bd_Part = val;}

    int GetBd_Part() const
       {return Bd_Part;}

    void SetCellIndex(int val)
       {CellIndex = val;}

    int GetCellIndex() const
       {return CellIndex;}

    void SetGlobalCellNo(int val)
       {GlobalCellNo = val;}

    int GetGlobalCellNo() const
       {return GlobalCellNo;}

    void SetLocalCellNo(int val)
       {LocalCellNo = val;}

    int GetLocalCellNo() const 
       {return LocalCellNo;}

    void SetRegionID(int val)
       {region = val;}

    int GetRegionID() const
       {return region;}       

    void SetAsLayerCell(int val)
       {LayerCell = val;}

    int IsLayerCell() const
       {return LayerCell;}          

    void SetNormalOrientation();
    
    int GetNormalOrientation(int i) const
    { return normalOrientation[i]; }

       
#ifdef  _MPI

    void SetClipBoard_Par(int value)
    { ClipBoard_Par=value; }
    
    int GetClipBoard_Par()
    { return ClipBoard_Par; }
    
    void SetSubDomainNo(int val)
       {SubDomainNumber = val;}

    int GetSubDomainNo() const
       {return SubDomainNumber;}

    void SetAsOwnCell()
     { OwnCell=TRUE; }

    void SetAsDependentCell()
     { DependentCell=TRUE; }

    void SetAsHaloCell()
     { 
      OwnCell=FALSE;
      HaloCell=TRUE;
     }

    bool IsHaloCell() const
     {  return HaloCell; }

    void SetAsSubDomainInterfaceCell()
     { SubDomainInterfaceCell=TRUE; }

    bool IsSubDomainInterfaceCell() const
     {  return SubDomainInterfaceCell; }

    void SetAsCrossEdgeCell()
     { CrossEdgeCell=TRUE; }

    bool IsCrossEdgeCell() const
     {  return CrossEdgeCell; }


    void SetAsCrossVertexCell()
     { CrossVertexCell=TRUE; }

    bool IsCrossVertexCell() const
     {  return CrossVertexCell; }

    bool IsDependentCell() const
     {  return DependentCell; }

    void SetN_NeibProcesses(int n)
     { N_NeibProcesses = n; }

    int GetN_NeibProcesses() const
     { return N_NeibProcesses; }

    void SetNeibProcessesIds(int *Neiblist);

    int *GetNeibProcessesIds() const
     { return NeibProcessesIds; } 
     
#endif
};

#endif