src/monomial_ordering.cpp

Go to the documentation of this file.

#ifndef mmx_realroot_monomial_ordering_cpp
#define mmx_realroot_monomial_ordering_cpp
#include <math.h>
#include <iostream>
#include <realroot/monomial_ordering.hpp>
namespace mmx {

  bool DegRevLex::less(const int *m1, int n1, const int *m2, int n2) const {
    
    int d1=0, d2=0;
    for(int i=0;i<n1;i++) d1+=m1[i];
    for(int i=0;i<n2;i++) d2+=m2[i];
  
    //    std::cout<< "less d1:= "<<d1<<std::endl;
    //    std::cout<< "less d2:= "<<d2<<std::endl;
    if (d1 < d2 ) return true;
    else if (d1 > d2) return false;
    else {
      //     if(n1 <=0 && n2>0) return false;
      //     if(n1 <=0 && n2<=0) return false;
      //     if(n2 <=0 && n1>0) return true;
     int n = (n1>n2?n1:n2); 
     //     std::cout<< "less n:= "<<n<<std::endl;
     int i;
     for(i = (n>0?n:0); i>=0 && (i<n1?m1[i]:0) == (i<n2?m2[i]:0);i--);
     //     std::cout<< "less i:= "<<i<<std::endl;
     if(i < 0) 
       return false;
     else 
       return ((i<n1?m1[i]:0) > (i<n2?m2[i]:0));
   }
  }

  bool LexRevDegree::less(const int *m1, int n1, const int *m2, int n2) const {
    
    int d1=0, d2=0;
    for(int i=0;i<n1;i++) d1+=m1[i];
    for(int i=0;i<n2;i++) d2+=m2[i];
  
    //    std::cout<< "less d1:= "<<d1<<std::endl;
    //    std::cout<< "less d2:= "<<d2<<std::endl;
    if (d1 < d2 ) return false;
    else if (d1 > d2) return true;
    else {
      //     if(n1 <=0 && n2>0) return false;
      //     if(n1 <=0 && n2<=0) return false;
      //     if(n2 <=0 && n1>0) return true;
     int n = (n1>n2?n1:n2); 
     //     std::cout<< "less n:= "<<n<<std::endl;
     int i;
     for(i = (n>0?n:0); i>=0 && (i<n1?m1[i]:0) == (i<n2?m2[i]:0);i--);
     //     std::cout<< "less i:= "<<i<<std::endl;
     if(i < 0) 
       return false;
     else 
       return ((i<n1?m1[i]:0) > (i<n2?m2[i]:0));
   }
  }


  bool Lex::less(const int *m1, int n1, const int *m2, int n2) const {
    int i = 0;
    for(i=0; i< (n1>n2?n1:n2) && (i<n1?m1[i]:0) == (i<n2?m2[i]:0); i++);
    
    return ((i<n1?m1[i]:0)< (i<n2?m2[i]:0));
  }
//====================================================================== 
} //namespace mmx
#endif //mmx_realroot_monomial_ordering_cpp

---