Changeset 24a96f8 in git for factory/facFqFactorize.cc

Timestamp:

Nov 24, 2014, 10:40:02 PM (9 years ago)

Author:

Martin Lee <martinlee84@…>

Branches:

(u'spielwiese', '2a584933abf2a2d3082034c7586d38bb6de1a30a')

Children:

9905e9eeb3329e6c42fcc9b2c0972cdf89fa3202

Parents:

f4fc9b40c900a9c6dce877b3b60856d2c433e962

git-author:

Martin Lee <martinlee84@web.de>2014-11-24 22:40:02+01:00

git-committer:

Martin Lee <martinlee84@web.de>2014-11-25 20:39:45+01:00

Message:

fix: LCHeuristic

File:

• factory/facFqFactorize.cc (modified) (3 diffs)

factory/facFqFactorize.cc

@@ -2621 +2621 @@
       iter2.getItem() *= power (xx, degree (LC (iter.getItem(),1), xx));
   }
-  CanonicalForm tmp;
+  CanonicalForm tmp, quot1, quot2, quot3;
   iter2= vars1;
   for (iter= leadingCoeffs[lengthAeval-1]; iter.hasItem(); iter++, iter2++)
@@ -2661 +2661 @@
             {
               tmp= ii.getItem().factor();
-              iter2.getItem() /= tmp;
-              CFListIterator iter3= evaluation;
-              for (int jj= A.level(); jj > 2; jj--, iter3++)
-                tmp= tmp (iter3.getItem(), jj);
-              if (!tmp.inCoeffDomain())
+              if (fdivides (tmp, iter2.getItem(), quot1))
               {
-                int index3= 1;
-                for (iter3= biFactors; iter3.hasItem(); iter3++, index3++)
+                CFListIterator iter3= evaluation;
+                for (int jj= A.level(); jj > 2; jj--, iter3++)
+                  tmp= tmp (iter3.getItem(), jj);
+                if (!tmp.inCoeffDomain())
                 {
-                  if (index3 == index2)
+                  int index3= 1;
+                  for (iter3= biFactors; iter3.hasItem(); iter3++, index3++)
                   {
-                    iter3.getItem() /= tmp;
-                    iter3.getItem() /= Lc (iter3.getItem());
-                    break;
+                    if (index3 == index2)
+                    {
+                      if (fdivides (tmp, iter3.getItem(), quot2))
+                      {
+                        if (fdivides (ii.getItem().factor(), A, quot3))
+                        {
+                          A               = quot3;
+                          iter2.getItem() = quot2;
+                          iter3.getItem() = quot3;
+                          iter3.getItem() /= Lc (iter3.getItem());
+                          break;
+                        }
+                      }
+                    }
                   }
                 }
               }
-              A /= ii.getItem().factor();
             }
           }
@@ -2699 +2708 @@
             {
               tmp= power (ii.getItem().factor(), ii.getItem().exp());
-              iter2.getItem() /= tmp;
-              A /= tmp;
-              CFListIterator iter3= evaluation;
-              for (int jj= A.level(); jj > 2; jj--, iter3++)
-                tmp= tmp (iter3.getItem(), jj);
-              if (!tmp.inCoeffDomain())
+              if (fdivides (tmp, A, quot1))
               {
-                int index3= 1;
-                for (iter3= biFactors; iter3.hasItem(); iter3++, index3++)
+                if (fdivides (tmp, iter2.getItem()))
                 {
-                  if (index3 == index2)
+                  CFListIterator iter3= evaluation;
+                  for (int jj= A.level(); jj > 2; jj--, iter3++)
+                    tmp= tmp (iter3.getItem(), jj);
+                  if (!tmp.inCoeffDomain())
                   {
-                    iter3.getItem() /= tmp;
-                    iter3.getItem() /= Lc (iter3.getItem());
-                    break;
+                    int index3= 1;
+                    for (iter3= biFactors; iter3.hasItem(); iter3++, index3++)
+                    {
+                      if (index3 == index2)
+                      {
+                        if (fdivides (tmp, iter3.getItem(), quot3))
+                        {
+                          A               =  quot1;
+                          iter2.getItem() =  quot2;
+                          iter3.getItem() =  quot3;
+                          iter3.getItem() /= Lc (iter3.getItem());
+                          break;
+                        }
+                      }
+                    }
                   }
                 }