Rlex.hpp

Go to the documentation of this file.

#ifndef _Rlex_H_
#define _Rlex_H_

template <class M> struct Rlex
{
  static bool less (const M &, const M &);
  bool operator() (const M &, const M &);
};

template <class M>
bool Rlex<M>::less (const M & m1, const M & m2) {
   int n1 = m1.Lvar(), n2 =m2.Lvar();
   if(n1 <0 && n2>0) return false;
   if(n2 <0 && n1>0) return true;

   int dm1=0,dm2=0,firstdiff=0,tmpdeg1,tmpdeg2;
   int n = std::max(n1,n2);
   for(int i=0;i<=n;i++)
   {
       tmpdeg1=m1.GetDegree(i);
       tmpdeg2=m2.GetDegree(i);
       if(tmpdeg1>tmpdeg2 && firstdiff==0) firstdiff=1;
       if(tmpdeg2>tmpdeg1 && firstdiff==0) firstdiff=-1;
       dm1+=tmpdeg1;
       dm2+=tmpdeg2;
   }
   if (dm1 < dm2) return false;
   else if (dm1 > dm2) return true;
   return (firstdiff==1);

}

template <class M>
bool Rlex<M>::operator() (const M & m1, const M & m2) {
        return less(m1,m2);
}

//template <class T1,class T2>
//bool operator< (const Monom<T1,T2> & m1, const Monom<T1,T2> & m2) {
//    return Rlex<Monom<T1,T2> >::less(m1,m2);
//}

#endif // //_Rlex_H_