RealBox.H

Go to the documentation of this file.

 
#ifndef _REALBOX_H_
#define _REALBOX_H_
 
#include <array>
#include "RealTensor2.H"
#include "Vec3.H"
#include "Types.H"
 
namespace criptic {
 
  class RealBox {
 
  public:
    
    RealBox(const RealVec& lo_ = zeroVec,
        const RealVec& hi_ = zeroVec) : lo(lo_), hi(hi_) { }
 
    RealBox(const RealVec* pts, const IdxType npts) {
      lo = pts[0];
      hi = pts[0];
      for (IdxType n=1; n<npts; n++) {
    lo = min(lo, pts[n]);
    hi = max(hi, pts[n]);
      }   
    }
 
    RealBox(const RealBox& rb1, const RealBox& rb2) {
      lo = min(rb1.lo, rb2.lo);
      hi = max(rb1.hi, rb2.hi);
    }
    
    bool contains(const RealVec& x) const {
      return lo[0] <= x[0] && x[0] <= hi[0] &&
    lo[1] <= x[1] && x[1] <= hi[1] &&
    lo[2] <= x[2] && x[2] <= hi[2];
    }
 
    std::array<RealVec,8> corners() const {
      std::array<RealVec,8> c;
      int ptr = 0;
      for (int i=0; i<2; i++) {
    for (int j=0; j<2; j++) {
      for (int k=0; k<2; k++) {
        c[ptr][0] = (i == 0) ? lo[0] : hi[0];
        c[ptr][1] = (j == 0) ? lo[1] : hi[1];
        c[ptr][2] = (k == 0) ? lo[2] : hi[2];
        ptr++;
      }
    }
      }
      return c;
    }
 
    RealVec size() const { return hi - lo; }
 
    RealBox rotate(const RealTensor2& rot) const {
      std::array<RealVec,8> c = corners();
      for (int i=0; i<8; i++) c[i] = rot * c[i];
      return RealBox(c.data(), 8);
    }
 
    void gaussExtrema(const RealTensor2& eta,
              Real& gMin,
              Real& gMax) const;
 
    // Operations with vectors
 
    RealBox operator+(const RealVec& v) const {
      RealBox rb;
      rb.lo = lo + v;
      rb.hi = hi + v;
      return rb;
    }
 
    RealBox operator-(const RealVec& v) const {
      RealBox rb;
      rb.lo = lo - v;
      rb.hi = hi - v;
      return rb;
    }
    
 
    // Data
    RealVec lo;  
    RealVec hi;  
  };
  
}
 
// Overload the >> and << operators for convenience
std::istream& operator>>(std::istream& is,
             criptic::RealBox& b);
 
std::ostream& operator<<(std::ostream& os,
             const criptic::RealBox& b);
 
 
#endif
// _REALBOX_H_