SQVlex.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: SQVlex.hpp,v 1.1.1.1 2006/10/06 08:01:40 trebuche Exp $
**********************************************************************/
#ifndef _SQVlex_H_
#define _SQVlex_H_

template <class M> struct SQVlex {
  static bool less (const M &, const M &);
  bool operator() (const M &, const M &);
};
 
template <class M>
bool SQVlex<M>::less(const M & m1, const M & m2) {

   int n1 = m1.Lvar(), n2 =m2.Lvar();
   int n = max(n1,n2);
   //cout <<m1<<" "<<m2<<" "<<n1<<" "<<n2 <<endl;
   if(n1 <0 && n2>0) return false;
   if(n2 <0 && n1>0) return true;
   if(n1 <0 && n2<0) return false;

   int d1=0,d2=0, k, c1, c2;
   for(int i =0 ; i<= n1; i++) {c1 = max(d1,(k=m1.GetDegree(i)));d1+=k;}
   for(int i =0 ; i<= n2; i++) {c2 = max(d2,(k=m2.GetDegree(i)));d2+=k;}

   if      (c1 < c2) return false;
   else if (c1 > c2) return true;
   else if (d1 < d2) return false;
   else if (d1 > d2) return true;
   else {
     int i=n;
     while( m1.GetDegree(i) == m2.GetDegree(i) && i>=0) i--;
     if(i <0) 
       return false;
     else 
     return (m1.GetDegree(i)> m2.GetDegree(i));
   }
}

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

template <class M>
bool operator<(const M & m1, const M & m2) {return SQVlex<M>::less(m1,m2);}
#endif // _SQVlex_H_