monomial_glue.hpp

Go to the documentation of this file.

/******************************************************************************
* MODULE     : monomial_glue.hpp
* DESCRIPTION: Glue for monomials
* COPYRIGHT  : (C) 2008 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_MONOMIAL_GLUE_HPP
#define __MMX_MONOMIAL_GLUE_HPP
#include <basix/category_meta.hpp>
#include <basix/tuple.hpp>
#include <basix/vector.hpp>
#include <realroot/monomial.hpp>

#define TMPL 
//template<typename C>
#define MONOMIAL monom<int>

namespace mmx {
TMPL bool exact_eq (const Monomial& v1, const Monomial& v2) {return v1==v2;}
TMPL bool exact_neq(const Monomial& v1, const Monomial& v2) {return v1!=v2;}

TMPL unsigned hash (const Monomial& v)  {
  register unsigned i, h= 214365^hash(v.coeff()), n= v.size();
  for (i=0; i<n; i++)
      h= (h<<1) ^ (h<<5) ^ (h>>27) ^ hash(v[i]);
  return h;
}
TMPL unsigned exact_hash (const Monomial& m) {return hash(m);}
TMPL unsigned soft_hash (const Monomial& m) {return hash(m);}

TMPL syntactic
flatten (const MONOMIAL& v) {
  syntactic m= flatten(v.coeff());
  for(unsigned i=0;i<v.nbvar();i++)
    m = m*pow(syntactic(MONOMIAL::var[i].data()),v[i]);
  return m;
}


TMPL MONOMIAL 
monomial_from_vector (const vector<int> & t) {
  nat sz=N(t);
  MONOMIAL r(1,sz, AsSize());
  for(nat i=0;i<sz;i++)
    r.set_expt(i,t[i]);
  return r;
}

TMPL vector<int> 
vector_from_monomial (const MONOMIAL& m) {
  vector<int> res;
  for(nat i=0;i<m.size();i++) res<<m[i];
  return res;
}

TMPL inline int  N(const MONOMIAL& v) { return v.size(); }
TMPL inline int size (const MONOMIAL& v) { return N(v); }

template<typename T, typename F, typename TR, typename FR>
struct as_helper<monom<T,TR>, monom<F,FR> >
{
  static inline monom<T,TR>
  cv (const monom<F,FR>& x) {
    monom<T,TR> r(as<T>(x.coeff()),x.size(),AsSize());
    for (unsigned  i=0; i<x.size(); i++) 
      r.set_expt(i,x[i]);
    return r;
  }
};

//====================================================================
} // namespace mmx
#undef TMPL
#undef Monomial
#endif // __MMX_MONOMIAL_GLUE_HPP