Changeset cec9624 in git for libpolys

Timestamp:

Jul 11, 2019, 2:59:38 PM (5 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

dffcf06de72c5d8f3347abee6740a159f6cabd41

Parents:

18296cd9190bb2f8c652446aa5cf3a783ad2cb66

Message:

add: pp_DivideM, pp_Divide

Location:

libpolys/polys/monomials

Files:

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

libpolys/polys/monomials/p_polys.cc

@@ -1610 +1610 @@
   }
   p_Delete(&b, r);
+  return result;
+}
+
+poly pp_DivideM(poly a, poly b, const ring r)
+{
+  if (a==NULL) { return NULL; }
+  poly result=a;
+
+  if(!p_IsConstant(b,r))
+  {
+    if (rIsLPRing(r))
+    {
+      WerrorS("not implemented for letterplace rings");
+      return NULL;
+    }
+    poly prev=NULL;
+    a=p_Copy(a,r);
+    while (a!=NULL)
+    {
+      if (p_DivisibleBy(b,a,r))
+      {
+        p_ExpVectorSub(a,b,r);
+        prev=a;
+        pIter(a);
+      }
+      else
+      {
+        if (prev==NULL)
+        {
+          p_LmDelete(&result,r);
+          a=result;
+        }
+        else
+        {
+          p_LmDelete(&pNext(prev),r);
+          a=pNext(prev);
+        }
+      }
+    }
+  }
+  if (result!=NULL)
+  {
+    number inv=pGetCoeff(b);
+    //if ((!rField_is_Ring(r)) || n_IsUnit(inv,r->cf))
+    if (rField_is_Zp(r))
+    {
+      inv = n_Invers(inv,r->cf);
+      __p_Mult_nn(result,inv,r);
+      n_Delete(&inv, r->cf);
+    }
+    else
+    {
+      result = p_Div_nn(result,inv,r);
+    }
+  }
   return result;
 }

libpolys/polys/monomials/p_polys.h

@@ -1990 +1990 @@
 poly      p_MDivide(poly a, poly b, const ring r);
 poly      p_DivideM(poly a, poly b, const ring r);
+poly      pp_DivideM(poly a, poly b, const ring r);
 poly      p_Div_nn(poly p, const number n, const ring r);