include/shape/box_face.hpp

Go to the documentation of this file.

#ifndef shape_boxface_hpp
#define shape_boxface_hpp

# include <realroot/Seq.hpp>
# include <shape/vertex.hpp>
# include <shape/bounding_box.hpp>
# include <shape/point_with_idx.hpp>

# define TMPL template<class C, class V>
# define SELF box_face<C,V>
//====================================================================
namespace mmx {
  namespace shape {
    
    TMPL 
    struct box_face {
      typedef vertex<C,3,V>              Point; 
      typedef bounding_box<C,V>        BoundingBox;  
      box_face() {
        m_v<<2;m_s<<0;m_var<<0<<1;
      }
      box_face(int v1, int i1) {  
        m_v<<v1; m_s<<i1;
        switch(v1) {
        case 0:
          m_var<<1<<2;
        case 1:
          m_var<<0<<2;
        default:
          m_var<<0<<1;
        }
      }

      box_face(int v1, int i1, int v2, int i2) {
        m_v<<v1<<v2;
        m_s<<i1<<i2;
        m_var<<(3-v1-v2);
      }
   
      Point* new_point(const BoundingBox & bx, double u) const {
        double pt[3];
        pt[m_var[0]]= u;
        pt[m_v[0]] = (bx)(m_v[0],m_s[0]);
        pt[m_v[1]] = (bx)(m_v[1],m_s[1]);
        return new Point(pt[0],pt[1],pt[2]);
      }

      Point* new_point(const BoundingBox & bx, double u, double v) const {
        double pt[3];
        pt[m_var[0]]= bx(m_var[0],0)+ u*(bx(m_var[0],1)-bx(m_var[0],0));        
        pt[m_var[1]]= bx(m_var[1],0)+ v*(bx(m_var[1],1)-bx(m_var[1],0));        
        pt[m_v[0]] = (bx)(m_v[0],m_s[0]);
        return new Point(pt[0],pt[1],pt[2], m_v[0]);
      }

      bool   is_valid (const BoundingBox & bx, double u, double eps) const {
        if(m_var[0]==0) {
          //      return true;
          if (m_s[0]==0)
            return u-eps>(bx)(m_var[0],0);
          else
            return u+eps<(bx)(m_var[0],1);
        } else {
          if (m_s[0]==1)
            return u-eps>(bx)(m_var[0],0);
          else
            return u+eps<(bx)(m_var[0],1);
        }
      }

      bool   is_valid_scale (double u, double eps) const {
        if(m_var[0]==0)
          return u-eps>0;
        else
          return u+eps<1;
      }
      
      bool is_valid (const Point& p, const BoundingBox & bx, double eps=0.000001) const {

        if (m_s[0] ==0 && p[m_v[0]]<bx(m_v[0],0)-eps) return false;
        if (m_s[0] ==1 && p[m_v[0]]>bx(m_v[0],1)+eps) return false;
        return true;
      }


      int cvar(unsigned i)  const { return m_v[i]; }
      int side(unsigned i)  const { return m_s[i]; }
      int fvar(unsigned i)  const { return m_var[i]; }

      double lower(const BoundingBox & bx) const { 
        return (bx)(m_var[0],0); 
      }
      double upper(const BoundingBox & bx) const { 
        return (bx)(m_var[0],1); 
      }

    public:
      Seq<int> m_v, m_s, m_var;
    };

  
    
  }//namespace shape
}//namespace mmx;
# undef TMPL
# undef SELF
# endif //shape_box_face_hpp

---