/Users/jab/Documents/Teaching/02585/GEL2_and_demos/GEL/src/Geometry/AABox.h

#ifndef __GEOMETRY_AABOX__H
#define __GEOMETRY_AABOX__H

#include <iostream>
#include <vector>
#include "Triangle.h"

namespace Geometry
{
  const float DIST_THRESH = 5.0e-4f;

  class AABox
  {
        CGLA::Vec3f pmin, pmax, interior_point;
  public:

        AABox() {}

        AABox(const CGLA::Vec3f& _pmin, const CGLA::Vec3f& _pmax,
                  const CGLA::Vec3f& _interior_point):

        pmin(_pmin), pmax(_pmax), interior_point(_interior_point)
        {
      for(int i=0;i<3;++i)
                if((pmax[i]-pmin[i]) < DIST_THRESH)
                {
                        pmax[i] += DIST_THRESH/2.0f;
                        pmin[i] -= DIST_THRESH/2.0f;
                }
          assert(pmin.all_le(interior_point));
          assert(pmax.all_ge(interior_point));
        }

        const CGLA::Vec3f& get_pmin() const {return pmin;}

        const CGLA::Vec3f& get_pmax() const {return pmax;}

        bool intersect(const CGLA::Vec3f&, const CGLA::Vec3f&) const;

        void minmax_sq_dist(const CGLA::Vec3f& p, float& dmin, float& dmax) const;

        static AABox box_triangle(const Triangle&);

        static AABox box_and_split(const std::vector<Triangle>& invec,
                                                           std::vector<Triangle>& lvec,
                                                           std::vector<Triangle>& rvec);
                                                                                                                 
  };
}

#endif