HyperCuboid.h

Go to the documentation of this file.

#ifndef HYPERCUBOID_HH
#define HYPERCUBOID_HH

class HyperCuboid;

// HyperPlot includes
#include "Mint/MessageService.h"
#include "Mint/HyperPoint.h"
#include "Mint/HyperPlane.h"
#include "Mint/HyperLine.h"
#include "Mint/HyperVolume.h"

// Root includes

// std includes
#include <bitset>


class HyperCuboid {

  private:
  
  int _dimension;         
  HyperPoint _lowCorner;  
  HyperPoint _highCorner; 
  mutable std::vector<HyperPlane> _faces; 
  
  void updateFaceCash() const;

  public:
    
  const int& getDimension () const {return _dimension;}
  HyperCuboid(int dimension);
  HyperCuboid(const HyperPoint& lowCorner, const HyperPoint& highCorner);
  HyperCuboid(int dimension, double low, double high);

  void print(std::ostream& os = std::cout, int endline=1) const;

  HyperPoint getCenter() const;
  bool inVolume(const HyperPoint& coords) const;

  const HyperCuboid& inflateCuboid(double percent);

  double volume() const;
  
  double getPositiveIntersectionParameter(const HyperLine& line) const;
  HyperPoint getPositiveIntersectionPoint(const HyperLine& line) const;
  
  double getNegativeIntersectionParameter(const HyperLine& line) const;
  HyperPoint getNegativeIntersectionPoint(const HyperLine& line) const;

  //pick a random point from within the hypercuboid
  HyperPoint getRandomPoint(TRandom* random = gRandom) const;
  HyperPointSet getRandomPoints(int nPoints = 100, TRandom* random = gRandom) const;

  std::vector<HyperPlane> getConnectedHyperPlanes(const HyperPoint& point) const;
  std::vector<HyperPlane> getConnectedHyperPlanes(const HyperPointSet& point) const;
  std::vector<HyperPlane> getBoundaryHyperPlanes() const;

  HyperPoint getOppositePoint(const HyperPoint& point) const;
  HyperPointSet getConnectedVerticies(const HyperPointSet& pointSet) const;
  HyperPointSet getConnectedVerticies(const HyperPoint& point) const;
  HyperPointSet getVertices() const;
  
  HyperPointSet getEdgeCenters() const;


  HyperCuboid project        (std::vector<int> dims) const;
  HyperCuboid projectOpposite(std::vector<int> dims) const;
  bool        inVolume(const HyperPoint& coords, std::vector<int> dims) const;

  const HyperPoint& getLowCorner () const{ return _lowCorner ; }
  const HyperPoint& getHighCorner() const{ return _highCorner; }
  HyperPoint& getLowCorner () { return _lowCorner ; }
  HyperPoint& getHighCorner() { return _highCorner; }
  HyperCuboid splitAbove(int dimension, double fractionalSplitPoint) const;
  HyperCuboid splitBelow(int dimension, double fractionalSplitPoint) const;
  HyperVolume split(int dimension, double fractionalSplitPoint) const;
  void split(int dimension, double fractionalSplitPoint, HyperVolume& set) const;

  HyperPoint getWidth() const;
  double getWidth(int dim) const;

  bool setCorners(const HyperPoint& lowCorner, const HyperPoint& highCorner);

  bool operator ==(const HyperCuboid& other) const;
  bool operator !=(const HyperCuboid& other) const;

  ~HyperCuboid();

};



#endif