Changeset 8222d6 in git

Timestamp:

Sep 9, 2016, 2:42:40 PM (8 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

dd930e594c1cd299f2e2b78fe5444de64aede28a

Parents:

bebd5a0361ed499fa0769cad9424e1bb193acd7c

Message:

add: p_Compare

Location:

libpolys/polys/monomials

Files:

• p_polys.cc (modified) (1 diff)
• p_polys.h (modified) (1 diff)

libpolys/polys/monomials/p_polys.cc

@@ -4708 +4708 @@
 #include <polys/templates/p_Delete__T.cc>
 
+/***************************************************************/
+/*
+* compare the leading terms of a and b
+*/
+static int tCompare(const poly a, const poly b, const ring R)
+{
+  if (b == NULL) return(a != NULL);
+  if (a == NULL) return(-1);
+
+  /* a != NULL && b != NULL */
+  int r = p_LmCmp(a, b,R);
+  if (r != 0) return(r);
+  number h = n_Sub(pGetCoeff(a), pGetCoeff(b),R->cf);
+  r = -1 + n_IsZero(h,R->cf) + 2*n_GreaterZero(h,R->cf);   /* -1: <, 0:==, 1: > */
+  n_Delete(&h,R->cf);
+  return(r);
+}
+
+/*
+* compare a and b
+*/
+int p_Compare(const poly a, const poly b, const ring R)
+{
+  int r = tCompare(a, b,R);
+  if (r != 0) return(r);
+
+  poly aa = a;
+  poly bb = b;
+  while (r == 0 && aa != NULL && bb != NULL)
+  {
+    pIter(aa);
+    pIter(bb);
+    r = tCompare(aa, bb,R);
+  }
+  return(r);
+}
+

libpolys/polys/monomials/p_polys.h

@@ -2021 +2021 @@
 /// shifts components of the vector p by i
 void p_Shift (poly * p,int i, const ring r);
+/*----------------------------------------------------*/
+
+int p_Compare(const poly a, const poly b, const ring R);
 #endif // P_POLYS_H