PbRegularCartesianGrid2D.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_REGULARCARTESIANGRID2D_
#define  _PB_REGULARCARTESIANGRID2D_
 
#include <MeshViz/3Ddata/PbGrid2D.h>
 
class PiContouring;
class PoMeshLines;
class PoMeshFilled;
class PoMeshSides;
class PbPixelCell;
class PbQuadrangleCell;
 
class  PbRegularCartesianGrid2D : public PbGrid2D {
 
 
 public:
 
  PbRegularCartesianGrid2D(SbBool isDataDuplicate=TRUE) ;
 
  PbRegularCartesianGrid2D(const PbRegularCartesianGrid2D &mesh);
 
  PbRegularCartesianGrid2D(int num_x,int num_y, float x_min,float y_min, float x_max,float y_max, SbBool isDataDuplicate=TRUE) ;
 
  PbRegularCartesianGrid2D(int num_x,int num_y, float x_min,float y_min, float x_max,float y_max, const float *z, 
                           SbBool isDataDuplicate=TRUE) ;
 
  PbRegularCartesianGrid2D(int num_x,int num_y, 
                           const PbRegularCartesianGrid2D &other_mesh);
 
  ~PbRegularCartesianGrid2D() ;
 
  PbRegularCartesianGrid2D& operator=(const PbRegularCartesianGrid2D &mesh) ;
 
   friend int operator ==(const PbRegularCartesianGrid2D &m1, const PbRegularCartesianGrid2D &m2);
 
   friend int operator !=(const PbRegularCartesianGrid2D &m1, const PbRegularCartesianGrid2D &m2)
  { return !(m1 == m2); }
 
  void setGeometry(int num_x,int num_y, float x_min,float y_min, float x_max,float y_max);
 
  void setGeometry(int num_x,int num_y, float x_min,float y_min, float x_max,float y_max, const float *z);
 
  void getGeometry(int &num_x,int &num_y, float &x_min,float &y_min, float &x_max,float &y_max) const;
 
  void getGeometry(int &num_x,int &num_y, float &x_min,float &y_min, float &x_max,float &y_max, const float *&z) const;
 
  virtual const PbCell* findContainingCell(const SbVec3f &point, float tolerance, SbVec3f &pcoord, const PbCell *adjacent_cell) const;
 
  virtual const PbCell* findContainingCell(const SbVec3f &point, float tolerance, SbVec3f &pcoord) const;
 
  virtual const PbCell *getCell(int cell_index) const;
 
  virtual SbVec3f getNodeCoord(int nod_index) const;
 
  virtual SbVec3f getNodeCoord(int i, int j) const;
 
  virtual float getArea() const;
 
 
 private:
  void meshLines(PoMeshLines *_PoMeshLines) const; 
  void meshFilled(PoMeshFilled *_PoMeshFilled) const; 
  void meshSides(PoMeshSides *_PoMeshSides) const; 
  void meshContouring(PiContouring *_PiContouring, const float *z) const;
  void mesh2DVec(PoMesh2DVec *mesh_vec) const ;
  void getLimits(PoMesh2D *_PoMesh2D, Ppoint_list_list3 &limit) const;
 
  void getNodesCoord(const float *z, SoMFVec3f &vertices) const;
  virtual SbBool hasSidesNormals() const;
  virtual void getSidesNormals(SbVec3f *&normals, 
                               const PbArrayOfInt &limitSizes) const;
 
  void print(FILE *fp) const;
 
 private:
  void constructXYGeometry(float x_min,float y_min, float x_max,float y_max, SbBool is_data_duplicate);
  void setXYGeometry(int numX, int numY, const float *x, const float *y);
  void getXYGeometry(const float* &x, const float* &y) const;
 
 private:
  float m_x[2], m_y[2];
  float m_xmin, m_xmax, m_dx;
  float m_ymin, m_ymax, m_dy;
 
  void updateBoundingBox() const;
  void copy(const PbRegularCartesianGrid2D &mesh) ;
  void destroy();
 
  mutable PbPixelCell      *m_pixel;
  mutable PbQuadrangleCell *m_quadrangle;
} ;
 
/*---------------------------------------------------------------------------*/
 
inline void
PbRegularCartesianGrid2D::getGeometry(int &num_x,int &num_y, float &x_min,float &y_min, float &x_max,float &y_max) const
{
  num_x = numX; x_min = m_xmin; x_max = m_xmax;
  num_y = numY; y_min = m_ymin; y_max = m_ymax;
}
 
inline void
PbRegularCartesianGrid2D::getGeometry(int &num_x,int &num_y, float &x_min,float &y_min, float &x_max,float &y_max, const float *&zn) const
{
  num_x = numX; x_min = m_xmin; x_max = m_xmax;
  num_y = numY; y_min = m_ymin; y_max = m_ymax;
  zn = m_zCoord;
}
 
#endif /* _PB_REGULARCARTESIANGRID2D_ */