Region.h

Go to the documentation of this file.

 #ifndef REGION_H_
#define REGION_H_
 #ifdef __cplusplus
#ifdef SWIG
#include "Python.h"
#endif
#include "Surface.h"
#include "boundary_type.h"
#include <limits>
#endif

/* Forward declarations to resolve circular dependencies */
class Intersection;
class Union;
class Complement;
class Halfspace;

enum regionType {
  INTERSECTION,
  UNION,
  COMPLEMENT,
  HALFSPACE
};

class Region {

protected:

  regionType _region_type;

  std::vector<Region*> _nodes;

  Region* _parent_region;

public:
  Region();
  virtual ~Region();

  /* Functions for constructing the Region / other Regions */
  virtual void addNode(Region* node, bool clone=true);
  void removeHalfspace(Surface* surface, int halfspace);
  regionType getRegionType();
  void setParentRegion(Region* node);
  Region* getParentRegion();

  /* Getter functions */
  virtual std::vector<Region*> getNodes();
  virtual std::vector<Region*> getAllNodes();
  virtual std::map<int, Halfspace*> getAllSurfaces();

  /* Worker functions */
  virtual double getMinX();
  virtual double getMaxX();
  virtual double getMinY();
  virtual double getMaxY();
  virtual double getMinZ();
  virtual double getMaxZ();
  virtual boundaryType getMinXBoundaryType();
  virtual boundaryType getMaxXBoundaryType();
  virtual boundaryType getMinYBoundaryType();
  virtual boundaryType getMaxYBoundaryType();
  virtual boundaryType getMinZBoundaryType();
  virtual boundaryType getMaxZBoundaryType();
  virtual bool containsPoint(Point* point) =0;
  virtual double minSurfaceDist(LocalCoords* coords);
  virtual double minSurfaceDist(Point* point, double azim, double polar=M_PI_2);
  virtual Region* clone();
};

class Intersection : public Region {

public:
  Intersection();
  bool containsPoint(Point* point);
};

class Union : public Region {

public:
  Union();
  bool containsPoint(Point* point);
};

class Complement : public Region {

public:
  Complement();
  void addNode(Region* node, bool clone=true);
  bool containsPoint(Point* point);
};

class Halfspace : public Region {

public:

  Surface* _surface;

  int _halfspace;

  Halfspace(int halfspace, Surface* surface);
  Halfspace* clone();
  Surface* getSurface();
  int getHalfspace();
  void reverseHalfspace();
  std::map<int, Halfspace*> getAllSurfaces();
  double getMinX();
  double getMaxX();
  double getMinY();
  double getMaxY();
  double getMinZ();
  double getMaxZ();
  boundaryType getMinXBoundaryType();
  boundaryType getMaxXBoundaryType();
  boundaryType getMinYBoundaryType();
  boundaryType getMaxYBoundaryType();
  boundaryType getMinZBoundaryType();
  boundaryType getMaxZBoundaryType();
  bool containsPoint(Point* point);
  double minSurfaceDist(LocalCoords* coords);
  double minSurfaceDist(Point* point, double azim, double polar=M_PI_2);
};

class RectangularPrism : public Intersection {

public:
  RectangularPrism(double width_x, double width_y, double origin_x=0.,
                   double origin_y=0., double width_z=3E300,
                   double origin_z=0.);
  void setBoundaryType(boundaryType boundary_type);
};

#endif /* REGION_H_ */