Changeset 85bcd6 in git for libpolys/polys/monomials/p_polys.h

Timestamp:

Aug 1, 2012, 5:21:12 PM (11 years ago)

Author:

Oleksandr Motsak <motsak@…>

Branches:

(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'c987db42cd2ec943b97ac5746c99892ceddf909c')

Children:

d828d633862fd9b6ec0e9f7b7e95dc3d610527d0fcb5bfd57f91ed4e8f9022cad3b9a32ee509c968

Parents:

6ce030f7be001f1ad6bba7026beec10e5beac410

git-author:

Oleksandr Motsak <motsak@mathematik.uni-kl.de>2012-08-01 17:21:12+02:00

git-committer:

Oleksandr Motsak <motsak@mathematik.uni-kl.de>2012-08-01 19:55:33+02:00

Message:

added p_LmShortDivisibleByNoComp by analogy

add: also needed pDebugLmShortDivisibleByNoComp and p_LmDivisibleByNoComp (with 2 rings) for completeness

File:

• libpolys/polys/monomials/p_polys.h (modified) (3 diffs)

libpolys/polys/monomials/p_polys.h

@@ -1659 +1659 @@
   return _p_LmDivisibleByNoComp(a, b, r);
 }
+
+static inline BOOLEAN p_LmDivisibleByNoComp(poly a, const ring ra, poly b, const ring rb)
+{
+  p_LmCheckPolyRing1(a, ra);
+  p_LmCheckPolyRing1(b, rb);
+  return _p_LmDivisibleByNoComp(a, ra, b, rb);
+}
+
 static inline BOOLEAN p_LmDivisibleBy(poly a, poly b, const ring r)
 {
@@ -1692 +1700 @@
   return _p_LmDivisibleBy(a, r_a, b, r_b);
 }
+
 static inline BOOLEAN p_LmShortDivisibleBy(poly a, unsigned long sev_a,
                                     poly b, unsigned long not_sev_b, const ring r)
@@ -1709 +1718 @@
 #else
   return pDebugLmShortDivisibleBy(a, sev_a, r, b, not_sev_b, r);
+#endif
+}
+
+static inline BOOLEAN p_LmShortDivisibleByNoComp(poly a, unsigned long sev_a,
+                                           poly b, unsigned long not_sev_b, const ring r)
+{
+  p_LmCheckPolyRing1(a, r);
+  p_LmCheckPolyRing1(b, r);
+#ifndef PDIV_DEBUG
+  _pPolyAssume2(p_GetShortExpVector(a, r) == sev_a, a, r);
+  _pPolyAssume2(p_GetShortExpVector(b, r) == ~ not_sev_b, b, r);
+
+  if (sev_a & not_sev_b)
+  {
+    pAssume1(p_LmDivisibleByNoComp(a, b, r) == FALSE);
+    return FALSE;
+  }
+  return p_LmDivisibleByNoComp(a, b, r);
+#else
+  return pDebugLmShortDivisibleByNoComp(a, sev_a, r, b, not_sev_b, r);
 #endif
 }