GridParams.hpp

#ifndef GRID_PARAMS_HH
#define GRID_PARAMS_HH
 
#include "MCRCProgramming.hpp"
#include "flumy_export.hpp"
 
#include "Point3D.hpp"
#include "Point2D.hpp"
 
class ConvexPolyedra2D;
 
class FLUMY_EXPORT GridParams {
 
public:
  GridParams();
  GridParams(const MP_real & dx,
       const MP_real & dy,
       const MP_int  & nx,
       const MP_int  & ny,
       const MP_real & ox = 0.,
       const MP_real & oy = 0.);
  GridParams(const MP_real & dx,
       const MP_real & dy,
       const MP_real & dz,
       const MP_int  & nx,
       const MP_int  & ny,
       const MP_int  & nz,
       const MP_real & ox = 0.,
       const MP_real & oy = 0.,
       const MP_real & oz = 0.);
  GridParams(const GridParams& grid);
  virtual ~GridParams() {}
  GridParams & operator=(const GridParams &);
  virtual void reset(const GridParams & grid);
  void reset(const MP_real & dx,
             const MP_real & dy,
             const MP_int  & nx,
             const MP_int  & ny,
             const MP_real & ox = 0.,
             const MP_real & oy = 0.);
  void reset(const MP_real & dx,
             const MP_real & dy,
             const MP_real & dz,
             const MP_int  & nx,
             const MP_int  & ny,
             const MP_int  & nz,
             const MP_real & ox = 0.,
             const MP_real & oy = 0.,
             const MP_real & oz = 0.);
  MP_real getDx() const { return _dx; }
  MP_real getDy() const { return _dy; }
  MP_real getDz() const { return _dz; }
  MP_int  getNx() const { return _nx; }
  MP_int  getNy() const { return _ny; }
  MP_int  getNz() const { return _nz; }
  MP_real getOx() const { return _ox; }
  MP_real getOy() const { return _oy; }
  MP_real getOz() const { return _oz; }
  Point2D getOrigin() const;
  MP_int getNbNodes() const { return _nx * _ny * _nz; }
  void getApexRel(Point2D& p1, Point2D& p2, Point2D& p3, Point2D& p4) const;
  void setLocation(const MP_real & x,
                   const MP_real & y,
                   const MP_real & z = 0.);
  void setLocation(const Point3D & location);
  void setLocation(const Point2D & location);
  const Point3D & getLocation() const;
  MP_real getZLoc() const { return _location.get_elevation(); }
  void setRotation(const MP_real& rot);
  // Return the rotation of the grid
  MP_real getRotation() const;
  MP_real getLength() const { return (_nx-1) * _dx; }
  MP_real getWidth()  const { return (_ny-1) * _dy; }
  MP_real getHeight()  const { return is3D() ? (_nz-1) * _dz : MP_MAX_REAL_VALUE; }
  bool is3D() const { return (_nz > 1 && _dz < MP_MAX_REAL_VALUE); }
  bool isCompatibleWith(const GridParams & grid, bool check_geo = false) const;
 
#ifndef SWIG
  friend std::ostream & operator <<(std::ostream & os, const GridParams & grid);
#endif
 
  // 2D Point conversion
  Point2D & rel2Geo (Point2D & point) const;
  Point2D & rel2Grid(Point2D & point, bool round = true) const;
  Point2D & geo2Rel (Point2D & point) const;
  Point2D & geo2Grid(Point2D & point, bool round = true) const;
  Point2D & grid2Rel(Point2D & point) const;
  Point2D   grid2Rel(const MP_int & ix, const MP_int & iy) const;
  Point2D   grid2Geo(const MP_int & ix, const MP_int & iy) const;
 
  // 3D points conversion
  Point3D & rel2Geo(Point3D & point) const;
  Point3D & geo2Rel(Point3D & point) const;
 
  // 2D Polyedra conversion
  ConvexPolyedra2D & rel2Grid(ConvexPolyedra2D & polyedra) const;
  ConvexPolyedra2D & grid2Rel(ConvexPolyedra2D & polyedra) const;
 
  // For all following methods : the given point is assumed to be already put on grid
  bool is_on_grid(const Point2D & point) const;
  bool is_on_grid_nearest(const Point2D & point) const;
 
protected:
 
  MP_real _dx;
  MP_real _dy;
  MP_real _dz;
  MP_int _nx;
  MP_int _ny;
  MP_int _nz;
  MP_real _ox;
  MP_real _oy;
  MP_real _oz;
  Point3D _location;
  MP_real _rotation;
  MP_real _cosrot;
  MP_real _sinrot;
};
 
 
 
#endif /* GRID_PARAMS_HH */