Changeset feed237 in git for libpolys

Timestamp:

Apr 30, 2013, 3:48:57 PM (11 years ago)

Author:

Martin Lee <martinlee84@…>

Branches:

(u'spielwiese', '4a9821a93ffdc22a6696668bd4f6b8c9de3e6c5f')

Children:

81626e384392f28e77ebb8b19445f7959194d982

Parents:

08a955b3ca689a48afd7f7c5849c4179954a1ddb

git-author:

Martin Lee <martinlee84@web.de>2013-04-30 15:48:57+02:00

git-committer:

Martin Lee <martinlee84@web.de>2013-05-02 11:42:46+02:00

Message:

chg: Singular interface to factory absFactorize

Location:

libpolys/polys

Files:

• clapsing.cc (modified) (1 diff)
• clapsing.h (modified) (1 diff)

libpolys/polys/clapsing.cc

@@ -1025 +1025 @@
   return res;
 }
+#ifdef HAVE_NTL
+ideal singclap_absBiFactorize ( poly f, ideal & mipos, intvec ** exps, int & numFactors, const ring r)
+{
+  p_Test(f, r);
+
+  ideal res=NULL;
+
+  int offs = rPar(r);
+  if (f==NULL)
+  {
+    res= idInit (1, 1);
+    mipos= idInit (1, 1);
+    mipos->m[0]= convFactoryPSingTrP (Variable (offs), r); //overkill
+    (*exps)=new intvec (1);
+    (**exps)[0]= 1;
+    numFactors= 0;
+    return res;
+  }
+  CanonicalForm F( convSingTrPFactoryP( f, r) );
+
+  if (getNumVars (F) > 2)
+  {
+    WerrorS( feNotImplemented );
+    return res;
+  }
+  CFAFList absFactors= absFactorize (F);
+
+  int n= absFactors.length();
+  *exps=new intvec (n);
+
+  res= idInit (n, 1);
+
+  mipos= idInit (n, 1);
+
+  Variable x= Variable (offs);
+  Variable alpha;
+  int i= 0;
+  numFactors= 0;
+  int count;
+  CFAFListIterator iter= absFactors;
+  CanonicalForm lead= iter.getItem().factor();
+  if (iter.getItem().factor().inCoeffDomain())
+  {
+    i++;
+    iter++;
+  }
+  for (; iter.hasItem(); iter++, i++)
+  {
+    (**exps)[i]= iter.getItem().exp();
+    alpha= iter.getItem().minpoly().mvar();
+    if (iter.getItem().minpoly().isOne()) //TODO make sure isOn (SW_RATIONAL) == true
+      lead /= bCommonDen (iter.getItem().factor());
+    else
+      lead /= power (bCommonDen (iter.getItem().factor()), degree (iter.getItem().minpoly()));
+    res->m[i]= convFactoryPSingTrP (replacevar (iter.getItem().factor()*bCommonDen (iter.getItem().factor()), alpha, x), r);
+    if (iter.getItem().minpoly().isOne())
+    {
+      count= iter.getItem().exp();
+      mipos->m[i]= convFactoryPSingTrP (x,r);
+    }
+    else
+    {
+      count= iter.getItem().exp()*degree (iter.getItem().minpoly());
+      mipos->m[i]= convFactoryPSingTrP (replacevar (iter.getItem().minpoly(), alpha, x), r);
+    }
+    numFactors += count;
+  }
+
+  (**exps)[0]= 1;
+  res->m[0]= convFactoryPSingTrP (lead, r);
+  mipos->m[0]= convFactoryPSingTrP (x, r);
+  return res;
+}
+#endif
 ideal singclap_sqrfree ( poly f, intvec ** v , int with_exps, const ring r)
 {

libpolys/polys/clapsing.h

@@ -51 +51 @@
 matrix singntl_LLL(matrix A, const ring r);
 intvec* singntl_LLL(intvec* A, const ring r);
+ideal singclap_absBiFactorize ( poly f, ideal & mipos, intvec ** exps, int & n, const ring r);
 #endif
 #endif