MPolyFunction.hpp

Go to the documentation of this file.

/*********************************************************************
*      This file is part of the source code of BORDERBASIX software. *
*      (C) B. Mourrain,  INRIA                                       *
**********************************************************************
History: 
$Id: MPolyFunction.hpp,v 1.1.1.1 2006/10/06 08:01:40 trebuche Exp $
**********************************************************************/
#ifndef _MPolyFunction_hpp_
#define _MPolyFunction_hpp_
 
//template <class R, class O>int Lvar (const MPoly<R,O> &);

template <class P> inline
P Pdiv (const P & p, const typename P::monom_t & m) {
    P res;
    for (typename P::const_iterator i = p.begin(); i != p.end(); ++i)
        res += div<typename P::monom_t>(*i,m);
    return res;
}

template <class R, class O> inline
int Degree (const MPoly<R,O> & P, int i) {
        int deg = 0;
        for (typename MPoly<R,O>::const_iterator p = P.begin(); p != P.end(); ++p)
                {
                        int d = p->GetDegree(i);
                        if (d > deg)
                                deg = d;
                }
        return deg;
}

template <class R, class O> inline
MPoly<R,O> Diff (const MPoly<R,O> & P, int i) {
        MPoly<R,O> c,res=0;
        int d;
        int deg = 0;
        for (typename MPoly<R,O>::const_iterator p = P.begin(); p != P.end(); ++p)
                {
                        typename MPoly<R,O>::monom_t aux(*p);
                        d = aux.GetDegree(i);
                        c = aux.coeff;
                        //std::cout<<c<<std::endl;
                        if (d > deg)
                        {

                            aux.SetDegree(i,d-1);
                            aux.coeff *=d;
 //                           aux.SetCoeff((double)c);
                            res+=aux;
                            //cout<<"res"<<res<<endl;
                        }
                }
        return res;
}

template <class R, class O> inline
int NbCoef (const MPoly<R,O> & P) {
        int nc = 0;
        for (typename MPoly<R,O>::const_iterator p = P.begin(); p != P.end(); ++p)
          nc++;
        return nc;
}

template <class R, class O> inline
MPoly<R,O> Coeff(const MPoly<R,O> & P, int i) {
   MPoly<R,O> res;
   for (typename MPoly<R,O>::const_iterator j = P.begin(); j != P.end(); ++j)
      if (j->GetDegree(i))
         res += *j;
   return res;
}

template <class R, class O> inline
MPoly<R,O> Coeff(const MPoly<R,O> & P, int i, int d) {
        MPoly<R,O> res;
        for (typename MPoly<R,O>::const_iterator j = P.begin(); j != P.end(); ++j)
                if (j->GetDegree(i) == d) {
                        typename MPoly<R,O>::monom_t aux(*j);
                        aux.SetDegree(i,0);
                        res += aux;
                }
        return res;
}


#endif //_MPolyFunction_hpp_