Changeset f8fe72 in git for libpolys/polys

Timestamp:

Feb 9, 2022, 10:19:22 PM (2 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', '2a584933abf2a2d3082034c7586d38bb6de1a30a')

Children:

ed6db890624524bb0ee27dd818ebf9dc825cf0df

Parents:

b59b0e73729545f663afa83406856b5d0c5f3cba

Message:

opt: simplified flint-check

File:

• libpolys/polys/monomials/p_polys.cc (modified) (1 diff)

libpolys/polys/monomials/p_polys.cc

@@ -3818 +3818 @@
       pSetCoeff0(p1,c);
       poly h = pNext(p1);
-      while (h!=NULL)
-      {
-        c=n_Div(pGetCoeff(h),k,r->cf);
-        // no need to normalize: Z/p, R
-        // normalize already in nDiv: Q_a, Z/p_a
-        // remains: Q
-        if (rField_is_Q(r) && (!n_IsOne(c,r->cf))) n_Normalize(c,r->cf);
-        p_SetCoeff(h,c,r);
-        pIter(h);
+      if (rField_is_Zp(r))
+      {
+        if (r->cf->ch>32003)
+        {
+          number inv=n_Invers(k,r->cf);
+          while (h!=NULL)
+          {
+            c=n_Mult(pGetCoeff(h),inv,r->cf);
+            // no need to normalize
+            p_SetCoeff(h,c,r);
+            pIter(h);
+          }
+          n_Delete(&inv,r->cf);
+        }
+        else
+        {
+          while (h!=NULL)
+          {
+            c=n_Div(pGetCoeff(h),k,r->cf);
+            // no need to normalize
+            p_SetCoeff(h,c,r);
+            pIter(h);
+          }
+        }
+      }
+      else
+      {
+        while (h!=NULL)
+        {
+          c=n_Div(pGetCoeff(h),k,r->cf);
+          // no need to normalize: Z/p, R
+          // normalize already in nDiv: Q_a, Z/p_a
+          // remains: Q
+          if (rField_is_Q(r) && (!n_IsOne(c,r->cf))) n_Normalize(c,r->cf);
+          p_SetCoeff(h,c,r);
+          pIter(h);
+        }
       }
       n_Delete(&k,r->cf);