Changeset 6c4143 in git

Timestamp:

Oct 18, 2018, 12:58:43 PM (6 years ago)

Author:

Christian Eder

Branches:

(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'b4f17ed1d25f93d46dbe29e4b499baecc2fd51bb')

Children:

ab6ae3e4daabcf18b48dab4c572a7e169b83bbd1

Parents:

b2bceb6030806c7016ffc929e51bb5db6615e546

Message:

adds special handling of one-term-reductions over Z

File:

• kernel/GBEngine/kspoly.cc (modified) (1 diff)

kernel/GBEngine/kspoly.cc

@@ -95 +95 @@
 #endif
 
-  /* if (t2==NULL)           // Divisor is just one term, therefore it will
-   * {                       // just cancel the leading term
-   *   PR->LmDeleteAndIter();
-   *   if (coef != NULL) *coef = n_Init(1, tailRing->cf);
-   *   return 0;
-   * } */
+  if (t2==NULL)           // Divisor is just one term, therefore it will
+  {                       // just cancel the leading term
+    // adjust lead coefficient if needed
+    if (! n_IsOne(pGetCoeff(p2), tailRing->cf))
+    {
+      number bn = pGetCoeff(lm);
+      number an = pGetCoeff(p2);
+      int ct = ksCheckCoeff(&an, &bn, tailRing->cf);    // Calculate special LC
+      p_SetCoeff(lm, bn, tailRing);
+      if ((ct == 0) || (ct == 2))
+      PR->Tail_Mult_nn(an);
+      if (coef != NULL) *coef = an;
+      else n_Delete(&an, tailRing->cf);
+    }
+    PR->LmDeleteAndIter();
+    if (coef != NULL) *coef = n_Init(1, tailRing->cf);
+    return 0;
+  }
 
   p_ExpVectorSub(lm, p2, tailRing); // Calculate the Monomial we must multiply to p2