block.h

//
// Created by milinda on 4/21/17.
#ifndef SFCSORTBENCH_BLOCK_H
#define SFCSORTBENCH_BLOCK_H

#include "dendro.h"
#include "TreeNode.h"
#include <assert.h>
#include <treenode2vtk.h>
#define GHOST_WIDTH 3

namespace ot
{
   enum BlockType{UNSPECIFIED=0, INTRNL_INDEPENDENT, INTRNL_DEPENDENT, FULLY_INDEPENDENT, FULLY_DEPENDENT};

   class Block
   {



   private:
     ot::TreeNode m_uiBlockNode;

     unsigned int m_uiRotID;

     unsigned int m_uiRegGridLev;

     DendroIntL m_uiLocalElementBegin;

     DendroIntL m_uiLocalElementEnd;

     DendroIntL m_uiOffset;

     unsigned int m_uiSize1D;

     unsigned int m_uiPaddingWidth;

     unsigned int m_uiEleOrder;

     unsigned int m_uiSzX;

     unsigned int m_uiSzY;

     unsigned int m_uiSzZ;

     std::vector<unsigned int> m_uiBLK2DIAG;

     std::vector<unsigned int> m_uiBLKVERTX;

     unsigned int m_uiBlkElem_1D;

     bool m_uiIsInternal;
     
     BlockType m_uiBlkType;


   public:
    Block();

    Block(ot::TreeNode pNode, unsigned int rotID ,unsigned int regLev, unsigned int regEleBegin,unsigned int regEleEnd,unsigned int  eleOrder);

    ~Block();

    inline ot::TreeNode getBlockNode()const {return m_uiBlockNode;}

     inline unsigned int getRegularGridLev()const {return m_uiRegGridLev;}

     inline unsigned int getRotationID()const {return m_uiRotID;}

      inline DendroIntL getLocalElementBegin()const {return m_uiLocalElementBegin;}

     inline DendroIntL getLocalElementEnd()const {return m_uiLocalElementEnd;}

     inline unsigned int get1DPadWidth()const {return m_uiPaddingWidth;}

     inline unsigned int getElementOrder() const {return m_uiEleOrder;}

     void setOffset(DendroIntL offset);

     inline void setBlk2DiagMap(unsigned int owner,unsigned int dir, unsigned int id){m_uiBLK2DIAG[dir*(2*m_uiBlkElem_1D)+owner]=id;}

     inline void setBlk2VertexMap(unsigned int dir, unsigned int id){m_uiBLKVERTX[dir]=id;}

     inline void setIsInternal(bool isInternal){m_uiIsInternal=isInternal;}

     inline void setBlkType(BlockType btype) { m_uiBlkType = btype; }

     inline BlockType getBlockType() const {return m_uiBlkType;}

     inline bool isInternal(){return m_uiIsInternal;}

     inline void setBlkNodeFlag(unsigned int flag){m_uiBlockNode.setFlag(flag);};

     inline unsigned int getBlkNodeFlag() const { return (m_uiBlockNode.getFlag()>>NUM_LEVEL_BITS);};

     inline DendroIntL getOffset() const {return m_uiOffset;}

     inline unsigned int get1DArraySize() const {return m_uiSize1D ;}

     inline unsigned int getAllocationSzX() const {return m_uiSzX;}

     inline unsigned int getAllocationSzY() const {return m_uiSzY;}

     inline unsigned int getAllocationSzZ() const {return m_uiSzZ;}

     inline unsigned int getAlignedBlockSz() const
     {
       unsigned int tmp;
      ((m_uiSzX & ((1u<<DENDRO_BLOCK_ALIGN_FACTOR_LOG)-1))==0)? tmp=m_uiSzX : tmp=((m_uiSzX/(1u<<DENDRO_BLOCK_ALIGN_FACTOR_LOG))+1)*(1u<<DENDRO_BLOCK_ALIGN_FACTOR_LOG);
      return tmp*m_uiSzY*m_uiSzZ;
     }


     inline const unsigned int * getBlk2DiagMap() const {return &(*(m_uiBLK2DIAG.begin()));}

     inline const unsigned int * getBlk2VertexMap() const {return &(*(m_uiBLKVERTX.begin()));}

     inline void setAllocationSzX(unsigned int sz) {m_uiSzX=sz;}
     inline void setAllocationSzY(unsigned int sz) {m_uiSzY=sz;}
     inline void setAllocationSzZ(unsigned int sz) {m_uiSzZ=sz;}
     inline void setSiz1D(unsigned int sz){m_uiSize1D=sz;}

     inline unsigned int getElemSz1D() const { return m_uiBlkElem_1D;}


     inline const std::vector<unsigned int >& getBlk2DiagMap_vec() const {return m_uiBLK2DIAG;}
     inline const std::vector<unsigned int >& getBlk2VertexMap_vec() const {return m_uiBLKVERTX;}

     double computeGridDx () const;

     double computeGridDy () const;

     double computeGridDz () const;

     double computeDx(const Point & d_min,const Point & d_max) const ;
     double computeDy(const Point & d_min,const Point & d_max) const ;
     double computeDz(const Point & d_min,const Point & d_max) const ;

     /*** @brief initialize the block diagonal map. */
     void initializeBlkDiagMap(const unsigned int value);

     /*** @brief initialize the block vertex neighbour map. */
     void initializeBlkVertexMap(const unsigned int value);




   };

} // end of namespace ot





#endif //SFCSORTBENCH_BLOCK_H