include/shape/algebraic_curve_glue.hpp

Go to the documentation of this file.

/******************************************************************************
* MODULE     : implicit_curve_glue.hpp
* DESCRIPTION: Glue for implicit curve
* 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_CURVE_GLUE_HPP
#define __MMX_SHAPE_ALGEBRAIC_CURVE_GLUE_HPP
#include <basix/glue.hpp>
#include <basix/tuple.hpp>
#include <realroot/ring_sparse_glue.hpp>
#include <shape/MGXK.hpp>
#include <shape/axel_glue.hpp>
#include <shape/bounding_box.hpp>
#include <shape/algebraic_curve.hpp>

#define shape_algebraic_curve shape::algebraic_curve<rational,shape::MGXK> 
#define Polynomial  polynomial< rational, with<Sparse,DegRevLex> >
#define CURVE shape_algebraic_curve


namespace mmx {
  namespace shape {

    inline bool operator ==(const CURVE& v1, const CURVE& v2) {return true;}
    inline bool operator !=(const CURVE& v1, const CURVE& v2) {return !(v1==v2);}
    inline bool exact_eq(const CURVE& v1, const CURVE& v2) {return v1==v2;}
    inline bool exact_neq(const CURVE& v1, const CURVE& v2) {return v1!=v2;}
    
    inline unsigned hash (const CURVE& 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 CURVE& m) {return hash(m);}
    inline unsigned soft_hash (const CURVE& m) {return hash(m);}

    inline syntactic flatten (const CURVE& s) 
    {
      using namespace shape;
      
      syntactic res = "AlgebraicCurve";
      if(s.nbequation()>1)
        return apply(res, mmx::flatten(s.equation(0)), mmx::flatten(s.equation(1)));
      else
        return apply(res, mmx::flatten(s.equation(0)));
    }
  }
  
//   inline CURVE shape_algebraic_curve_string (const string& e1)
//   {
//     return CURVE (as_charp(e1));
//   }

//   inline CURVE shape_algebraic_curve_string (const string& e1,const string& e2)
//   {
//     return CURVE (as_charp(e1), as_charp(e2));
//   }

//   inline CURVE algebraic_curve(const Polynomial& e1,const Polynomial& e2)
//   {
//     return CURVE(e1, e2);
//   }
  

 
//====================================================================
#undef CURVE
#undef Polynomial
//====================================================================
} // namespace mmx
#endif // __MMX_ALGEBRAIC_CURVE_GLUE_HPP

---