Changeset 1e4b53 in git for factory/cf_gcd_smallp.cc

Timestamp:

May 3, 2012, 3:53:02 PM (12 years ago)

Author:

Martin Lee <martinlee84@…>

Branches:

(u'spielwiese', 'ec94ef7a30b928574c0c3daf41f6804dff5f6b69')

Children:

60cc1289a79586c1732ba3885997c3cb49bc01d9

Parents:

0d56271c672440c06e14ca2528a5185bb37e1e6f

git-author:

Martin Lee <martinlee84@web.de>2012-05-03 15:53:02+02:00

git-committer:

Martin Lee <martinlee84@web.de>2012-05-07 14:18:40+02:00

Message:

chg: first test LC's before doing a full termination test in GCD

File:

• factory/cf_gcd_smallp.cc (modified) (6 diffs)

factory/cf_gcd_smallp.cc

@@ -51 +51 @@
 TIMING_DEFINE_PRINT(newton_interpolation)
 
+bool
+terminationTest (const CanonicalForm& F, const CanonicalForm& G,
+                 const CanonicalForm& coF, const CanonicalForm& coG,
+                 const CanonicalForm& cand)
+{
+  CanonicalForm LCCand= abs (LC (cand));
+  if (LCCand*abs (LC (coF)) == abs (LC (F)))
+  {
+    if (LCCand*abs (LC (coG)) == abs (LC (G)))
+    {
+      if (abs (cand)*abs (coF) == abs (F))
+      {
+        if (abs (cand)*abs (coG) == abs (G))
+          return true;
+      }
+      return false;
+    }
+    return false;
+  }
+  return false;
+}
+
 static const double log2exp= 1.442695041;
 
@@ -805 +827 @@
         if (((degree (ppCoF,1)+degree (ppH,1) == bound1) &&
              (degree (ppCoG,1)+degree (ppH,1) == bound2) &&
-             (ppA == ppCoF*ppH && ppB == ppCoG*ppH)) || 
+             terminationTest (ppA, ppB, ppCoF, ppCoG, ppH)) || 
              (fdivides (ppH, ppA, ppCoF) && fdivides (ppH, ppB, ppCoG)))
         {
@@ -828 +850 @@
       else if (((degree (ppCoF,1)+degree (ppH,1) == bound1) &&
                 (degree (ppCoG,1)+degree (ppH,1) == bound2) &&
-                (ppA == ppCoF*ppH && ppB == ppCoG*ppH)) ||
+                terminationTest (ppA, ppB, ppCoF, ppCoG, ppH)) ||
                 (fdivides (ppH, ppA, ppCoF) && fdivides (ppH, ppB, ppCoG)))
       {
@@ -1255 +1277 @@
         if (((degree (ppCoF,1)+degree (ppH,1) == bound1) &&
              (degree (ppCoG,1)+degree (ppH,1) == bound2) &&
-             (ppA == ppCoF*ppH && ppB == ppCoG*ppH)) ||
+             terminationTest (ppA, ppB, ppCoF, ppCoG, ppH)) ||
              (fdivides (ppH, ppA, ppCoF) && fdivides (ppH, ppB, ppCoG)))
         {
@@ -1279 +1301 @@
       if (((degree (ppCoF,1)+degree (ppH,1) == bound1) &&
            (degree (ppCoG,1)+degree (ppH,1) == bound2) &&
-           (ppA == ppCoF*ppH && ppB == ppCoG*ppH)) ||
+           terminationTest (ppA, ppB, ppCoF, ppCoG, ppH)) ||
            (fdivides (ppH, ppA, ppCoF) && fdivides (ppH, ppB, ppCoG)))
         {
@@ -1773 +1795 @@
       if (((degree (ppCoF,1)+degree (ppH,1) == bound1) &&
            (degree (ppCoG,1)+degree (ppH,1) == bound2) &&
-           (ppA == ppCoF*ppH && ppB == ppCoG*ppH)) ||
+           terminationTest (ppA, ppB, ppCoF, ppCoG, ppH)) ||
            (fdivides (ppH, ppA, ppCoF) && fdivides (ppH, ppB, ppCoG)))
       {