PbGrid2D.h

/*=======================================================================
 *** THE CONTENT OF THIS WORK IS PROPRIETARY TO FEI S.A.S, (FEI S.A.S.),            ***
 ***              AND IS DISTRIBUTED UNDER A LICENSE AGREEMENT.                     ***
 ***                                                                                ***
 ***  REPRODUCTION, DISCLOSURE,  OR USE,  IN WHOLE OR IN PART,  OTHER THAN AS       ***
 ***  SPECIFIED  IN THE LICENSE ARE  NOT TO BE  UNDERTAKEN  EXCEPT WITH PRIOR       ***
 ***  WRITTEN AUTHORIZATION OF FEI S.A.S.                                           ***
 ***                                                                                ***
 ***                        RESTRICTED RIGHTS LEGEND                                ***
 ***  USE, DUPLICATION, OR DISCLOSURE BY THE GOVERNMENT OF THE CONTENT OF THIS      ***
 ***  WORK OR RELATED DOCUMENTATION IS SUBJECT TO RESTRICTIONS AS SET FORTH IN      ***
 ***  SUBPARAGRAPH (C)(1) OF THE COMMERCIAL COMPUTER SOFTWARE RESTRICTED RIGHT      ***
 ***  CLAUSE  AT FAR 52.227-19  OR SUBPARAGRAPH  (C)(1)(II)  OF  THE RIGHTS IN      ***
 ***  TECHNICAL DATA AND COMPUTER SOFTWARE CLAUSE AT DFARS 52.227-7013.             ***
 ***                                                                                ***
 ***                   COPYRIGHT (C) 1996-2014 BY FEI S.A.S,                        ***
 ***                        BORDEAUX, FRANCE                                        ***
 ***                      ALL RIGHTS RESERVED                                       ***
**=======================================================================*/
/*=======================================================================
** Author      : VSG (MMM YYYY)
**=======================================================================*/
 
#ifndef  _PB_GRID2D_
#define  _PB_GRID2D_
 
#include <MeshViz/3Ddata/PbMesh2D.h>
class  PbGrid2D : public PbMesh2D {
 
 
 public:
 
  PbGrid2D(SbBool is_data_duplicate=TRUE) : PbMesh2D(is_data_duplicate) {};
 
  PbGrid2D(const PbGrid2D &mesh);
 
  ~PbGrid2D();
 
  PbGrid2D& operator=(const PbGrid2D &mesh) ;
 
  void setGeometry(int num_x, int num_y, const float *x, const float *y);
 
  void setGeometry(int num_x, int num_y, const float *x, const float *y, const float *z);
 
  void getGeometry(int &num_x, int &num_y, const float* &x, const float* &y) const;
 
  void getGeometry(int &num_x, int &num_y, const float* &x, const float* &y, const float* &z) const;
 
  virtual SbVec3f getNodeCoord(int nod_index) const=0;
 
  virtual SbVec3f getNodeCoord(int i, int j) const=0;
 
  void getDim(int &num_x, int &num_y) const { num_x = numX; num_y = numY;};
 
  virtual void getNodeOwnerCellsInd(int nod_index, PbArrayOfInt &owner_cells) const;
 
  virtual void getAdjacentCellsIndByNode(int cell_index, PbArrayOfInt &adjacent_cells) const;
 
  virtual void getAdjacentCellsIndByFacet(int cell_index, PbArrayOfInt &adjacent_cells) const;
 
  void  getNodeIndices(int nod_index, int &i, int &j) const;
 
  void  getCellIndices(int cell_index, int &i, int &j) const;
 private:
  int * getCellNodeIds(int i, int j) const;
  int * getCellNodeIds(int cell_index) const;
  virtual SbVec3f getFirstCellNormal() const;
 
  void print(FILE *fp) const;
  virtual void getLimitIndices(PbArrayOfInt &nodesIndex, PbArrayOfInt &limitSizes) const;
 
 private:
  virtual void updateCellBoundingBox() const;
  PbGrid2D(int num_x, int num_y, SbBool isDataDuplicate=TRUE) ;
  PbGrid2D(int num_x, int num_y, const float *z, SbBool isDataDuplicate=TRUE) ;
  void setDim(int num_x, int num_y);
  virtual void setXYGeometry(int num_x, int num_y, const float *x, const float *y)=0;
  virtual void getXYGeometry(const float* &x, const float* &y) const=0;
  void setZCoordMatrix(int num_x, int num_y, const float *z);
  int   numX,numY;
  float **z;
 
  // node ids of the mutable private cell
  // each time the getCell or getTopoCell is called, m_cellNodeIds is updated
  mutable int m_cellNodeIds[4]; 
 
 private:
 
  void copy(const PbGrid2D &mesh) ;
  void destroy();
} ;
 
/*---------------------------------------------------------------------------*/
 
// calculate indices i,j that verify    nod_index = i*numY + j
inline void
PbGrid2D::getNodeIndices(int nod_index, int &i, int &j) const
{
  j = nod_index % numY;
  i = (nod_index-j) / numY;
}
 
// calculate indices i,j that verify    cell_index = i*(numY-1) + j
inline void
PbGrid2D::getCellIndices(int cell_index, int &i, int &j) const
{
  j = cell_index % (numY-1);
  i = (cell_index-j) / (numY-1);
}
 
// set dimensions of the grid
inline void
PbGrid2D::setDim(int num_x, int num_y)
{
  numX = num_x;
  numY = num_y;
  numMeshNodes = numX * numY;
  if (num_x > 1 && num_y > 1) 
    m_numMeshCells = (numX-1) * (numY-1);
  else
    // setDim can be used with (0,0) as argt !
    m_numMeshCells = 0;
}
 
inline int * 
PbGrid2D::getCellNodeIds(int cell_index) const
{
  int i, j;
  getCellIndices(cell_index, i, j);
  return getCellNodeIds(i,j);
}
 
inline int * 
PbGrid2D::getCellNodeIds(int i, int j) const
{
  m_cellNodeIds[0] = i*numY+j;
  m_cellNodeIds[1] = (i+1)*numY+j;
  m_cellNodeIds[2] = (i+1)*numY+j+1;
  m_cellNodeIds[3] = i*numY+j+1;
  return m_cellNodeIds;
}
 
#endif /* _PB_GRID2D_ */