include/shape/algebraic_surface_glue.hpp

Go to the documentation of this file.

/******************************************************************************
* MODULE     : implicit_surface_glue.hpp
* DESCRIPTION: Glue for implicit surface
* 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_ALGEBRAIC_SURFACE_GLUE_HPP
#define __MMX_SHAPE_ALGEBRAIC_SURFACE_GLUE_HPP
#include <basix/glue.hpp>
#include <basix/tuple.hpp>
#include <shape/MGXK.hpp>

#include <realroot/ring_sparse_glue.hpp>
#include <shape/axel_glue.hpp>
#include <shape/algebraic_curve_glue.hpp>
#include <shape/algebraic_surface.hpp>

#define shape_algebraic_surface shape::algebraic_surface<rational,shape::MGXK>
#define SURFACE    shape::algebraic_surface<rational,shape::MGXK>
#define CURVE      shape::algebraic_curve<rational,shape::MGXK>
#define Polynomial shape::algebraic_set<shape::MGXK>::Polynomial


namespace mmx {

  namespace shape {
    inline bool operator ==(const SURFACE& v1, const SURFACE& v2) {return v1.equation()==v2.equation();}
    inline bool operator !=(const SURFACE& v1, const SURFACE& v2) {return !(v1==v2);}
    inline bool exact_eq(const SURFACE& v1, const SURFACE& v2) {return v1==v2;}
    inline bool exact_neq(const SURFACE& v1, const SURFACE& v2) {return v1!=v2;}
    
    inline unsigned hash (const SURFACE& 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;
    }
    inline unsigned exact_hash (const SURFACE& m) {return hash(m);}
    inline unsigned soft_hash (const SURFACE& m) {return hash(m);}

    syntactic flatten (const SURFACE& s) 
    {
      using namespace shape;
      
      syntactic res = "AlgebraicSurface";
      
      return apply(res, mmx::flatten(s.equation()));
      
    }
  }
  //--------------------------------------------------------------------  
  CURVE shape_algebraic_surface_intersection (const SURFACE& s1, const SURFACE& s2)
  {
    //    using namespace shape;
    //    BOX *bx= s1.boundingBox() * s2.boundingBox();
    return CURVE(s1.equation(), s2.equation());
  }



//====================================================================
#undef Polynomial
#undef BOX
#undef CURVE
#undef SURFACE
//====================================================================
} // namespace mmx
#endif // __MMX_IMPLICIT_SURFACE_GLUE_HPP

---