chg: deleted some unused code

chg: joined precomputeLeadingCoeff from facFactorize.cc and facFqFactorize.cc
Singular · Oct 24, 2012 · b1abf09 · b1abf09
1 parent db9ebe9
commit b1abf09
Show file tree

Hide file tree

Showing 4 changed files with 39 additions and 637 deletions.
diff --git a/factory/facBivar.cc b/factory/facBivar.cc
@@ -184,171 +184,6 @@ CanonicalForm evalPoint (const CanonicalForm& F, int& i)
   } while (1);
 }
 
-CFList
-earlyFactorDetection0 (CanonicalForm& F, CFList& factors,int& adaptedLiftBound,
-                      DegreePattern& degs, bool& success, int deg)
-{
-  DegreePattern bufDegs1= degs;
-  DegreePattern bufDegs2;
-  CFList result;
-  CFList T= factors;
-  CanonicalForm buf= F;
-  CanonicalForm LCBuf= LC (buf, Variable (1));
-  CanonicalForm g, quot;
-  CanonicalForm M= power (F.mvar(), deg);
-  adaptedLiftBound= 0;
-  int d= degree (F) + degree (LCBuf, F.mvar());
-  for (CFListIterator i= factors; i.hasItem(); i++)
-  {
-    if (!bufDegs1.find (degree (i.getItem(), 1)))
-      continue;
-    else
-    {
-      g= i.getItem() (0, 1);
-      g *= LCBuf;
-      g= mod (g, M);
-      if (fdivides (LC (g), LCBuf))
-      {
-        g= mulMod2 (i.getItem(), LCBuf, M);
-        g /= content (g, Variable (1));
-        if (fdivides (g, buf, quot))
-        {
-          result.append (g);
-          buf= quot;
-          d -= degree (g) + degree (LC (g, Variable (1)), F.mvar());
-          LCBuf= LC (buf, Variable (1));
-          T= Difference (T, CFList (i.getItem()));
-
-          // compute new possible degree pattern
-          bufDegs2= DegreePattern (T);
-          bufDegs1.intersect (bufDegs2);
-          bufDegs1.refine ();
-          if (bufDegs1.getLength() <= 1)
-          {
-            result.append (buf);
-            break;
-          }
-        }
-      }
-    }
-  }
-  adaptedLiftBound= d + 1;
-  if (d < deg)
-  {
-    factors= T;
-    degs= bufDegs1;
-    F= buf;
-    success= true;
-  }
-  if (bufDegs1.getLength() <= 1)
-    degs= bufDegs1;
-  return result;
-}
-
-
-CFList
-henselLiftAndEarly0 (CanonicalForm& A, bool& earlySuccess, CFList&
-                    earlyFactors, DegreePattern& degs, int& liftBound,
-                    const CFList& uniFactors, const CanonicalForm& eval)
-{
-  int sizeOfLiftPre;
-  int * liftPre= getLiftPrecisions (A, sizeOfLiftPre, degree (LC (A, 1), 2));
-
-  Variable x= Variable (1);
-  Variable y= Variable (2);
-  CFArray Pi;
-  CFList diophant;
-  CFList bufUniFactors= uniFactors;
-  bufUniFactors.insert (LC (A, x));
-  CFMatrix M= CFMatrix (liftBound, bufUniFactors.length() - 1);
-  earlySuccess= false;
-  int newLiftBound= 0;
-  int smallFactorDeg= tmin (11, liftPre [sizeOfLiftPre- 1] + 1);//this is a tunable parameter
-  if (smallFactorDeg >= liftBound || degree (A,y) <= 4)
-    henselLift12 (A, bufUniFactors, liftBound, Pi, diophant, M);
-  else if (sizeOfLiftPre > 1 && sizeOfLiftPre < 30)
-  {
-    henselLift12 (A, bufUniFactors, smallFactorDeg, Pi, diophant, M);
-    earlyFactors= earlyFactorDetection0 (A, bufUniFactors, newLiftBound,
-                                        degs, earlySuccess,
-                                        smallFactorDeg);
-    if (degs.getLength() > 1 && !earlySuccess &&
-        smallFactorDeg != liftPre [sizeOfLiftPre-1] + 1)
-    {
-      if (newLiftBound > liftPre[sizeOfLiftPre-1]+1)
-      {
-        bufUniFactors.insert (LC (A, x));
-        henselLiftResume12 (A, bufUniFactors, smallFactorDeg,
-                            liftPre[sizeOfLiftPre-1] + 1, Pi, diophant, M);
-        earlyFactors= earlyFactorDetection0 (A, bufUniFactors, newLiftBound,
-                      degs, earlySuccess, liftPre[sizeOfLiftPre-1] + 1);
-      }
-    }
-    else if (earlySuccess)
-      liftBound= newLiftBound;
-    int i= sizeOfLiftPre - 1;
-    while (degs.getLength() > 1 && !earlySuccess && i - 1 >= 0)
-    {
-      if (newLiftBound > liftPre[i] + 1)
-      {
-        bufUniFactors.insert (LC (A, x));
-        henselLiftResume12 (A, bufUniFactors, liftPre[i] + 1,
-                            liftPre[i-1] + 1, Pi, diophant, M);
-        earlyFactors= earlyFactorDetection0 (A, bufUniFactors, newLiftBound,
-                      degs, earlySuccess, liftPre[i-1] + 1);
-      }
-      else
-      {
-        liftBound= newLiftBound;
-        break;
-      }
-      i--;
-    }
-    if (earlySuccess)
-      liftBound= newLiftBound;
-    //after here all factors are lifted to liftPre[sizeOfLiftPre-1]
-  }
-  else
-  {
-    henselLift12 (A, bufUniFactors, smallFactorDeg, Pi, diophant, M);
-    earlyFactors= earlyFactorDetection0 (A, bufUniFactors, newLiftBound,
-                                        degs, earlySuccess,
-                                        smallFactorDeg);
-    int i= 1;
-    while ((degree (A,y)/4)*i + 4 <= smallFactorDeg)
-      i++;
-    if (degs.getLength() > 1 && !earlySuccess)
-    {
-      bufUniFactors.insert (LC (A, x));
-      henselLiftResume12 (A, bufUniFactors, smallFactorDeg,
-                          (degree (A, y)/4)*i + 4, Pi, diophant, M);
-      earlyFactors= earlyFactorDetection0 (A, bufUniFactors, newLiftBound,
-                    degs, earlySuccess, (degree (A, y)/4)*i + 4);
-    }
-    while (degs.getLength() > 1 && !earlySuccess && i < 4)
-    {
-      if (newLiftBound > (degree (A, y)/4)*i + 4)
-      {
-        bufUniFactors.insert (LC (A, x));
-        henselLiftResume12 (A, bufUniFactors, (degree (A,y)/4)*i + 4,
-                            (degree (A, y)/4)*(i+1) + 4, Pi, diophant, M);
-        earlyFactors= earlyFactorDetection0 (A, bufUniFactors, newLiftBound,
-                      degs, earlySuccess, (degree (A, y)/4)*(i+1) + 4);
-      }
-      else
-      {
-        liftBound= newLiftBound;
-        break;
-      }
-      i++;
-    }
-    if (earlySuccess)
-      liftBound= newLiftBound;
-  }
-
-  return bufUniFactors;
-}
-
 CFList biFactorize (const CanonicalForm& F, const Variable& v)
 {
   if (F.inCoeffDomain())