include/shape/point_glue.hpp

Go to the documentation of this file.

/******************************************************************************
* MODULE     : point_glue.hpp
* DESCRIPTION: Glue for points
* COPYRIGHT  : (C) 2007 Bernard Mourrain
*******************************************************************************
* This software falls under the GNU general public license and comes WITHOUT
* ANY WARRANTY WHATSOEVER. See the file $TEXMACS_PATH/LICENSE for more details.
* If you don't have this file, write to the Free Software Foundation, Inc.,
* 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
******************************************************************************/
#ifndef __mmx_shape_point_glue_hpp
#define __mmx_shape_point_glue_hpp
#include <basix/glue.hpp>
#include <shape/axel_glue.hpp>
#include <shape/point.hpp>

#define TMPL   template<class C>
#define Point  shape::point<C>
#define Scalar C
#define Vector vector<Scalar>


namespace mmx {
  namespace shape {

    TMPL inline bool operator ==(const Point& v1, const Point& v2) {return true;}
    TMPL inline bool operator !=(const Point& v1, const Point& v2) {return !(v1==v2);}
    TMPL inline bool exact_eq(const Point& v1, const Point& v2) {return v1==v2;}
    TMPL inline bool exact_neq(const Point& v1, const Point& v2) {return v1!=v2;}
    
    TMPL inline unsigned hash (const Point& v) 
    {
      register unsigned i, h= 214365, n=1;
      for (i=0; i<n; i++) h= (h<<1) ^ (h<<5) ^ (h>>27) ;//^ hash(v[i]);
      return h;
    }
    TMPL inline unsigned exact_hash (const Point& m) {return hash(m);}
    TMPL inline unsigned soft_hash (const Point& m) {return hash(m);}

    TMPL inline syntactic flatten (const Point& p) {
      vector<syntactic> v;
      v<<mmx::flatten(p[0]);
      v<<mmx::flatten(p[1]);
      v<<mmx::flatten(p[2]);
      
      return mmx::flatten(v);
    }
  }
  //-------------------------------------------------------------------- 

  TMPL Point shape_point_tuple(const vector<C>& v) {
    Point p;
    for(unsigned i=0;i<N(v) && i<3; i++)
      p[i]=v[i];
    return p;
  }

  TMPL vector<C> operator-(const Point& p1, const Point& p2){
    Vector res((C)0,3);
    for(unsigned i=0;i<3; i++) res[i]=p1[i]-p2[i];
    return res;
  }

  TMPL Point  operator+(const Point& p, const vector<C>& v) {
    Point res;
    for(unsigned i=0;i<3; i++) res[i]=p[i]+v[i];
    return res;
  }
#undef TMPL
#undef Point
//====================================================================
} // namespace mmx
#endif // __mmx_point_glue_hpp

---