Changeset cd964c in git

Timestamp:

Feb 22, 2013, 11:24:13 AM (11 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', '4a9821a93ffdc22a6696668bd4f6b8c9de3e6c5f')

Children:

088e67fc66fb195ab2e00e54b3494d19574a1092

Parents:

e070751c183f5cec7c3ec5339a0b53ff93302566

git-author:

Hans Schoenemann <hannes@mathematik.uni-kl.de>2013-02-22 11:24:13+01:00

git-committer:

Hans Schoenemann <hannes@mathematik.uni-kl.de>2013-02-22 11:56:22+01:00

Message:

add: new function singclap_gcd_and_divide

Location:

libpolys/polys

Files:

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

libpolys/polys/clapsing.cc

@@ -99 +99 @@
 }
 
+void singclap_gcd_and_divide ( poly& f, poly& g, const ring r)
+{
+  poly res=NULL;
+
+  if (f!=NULL) p_Cleardenom(f, r);
+  if (g!=NULL) p_Cleardenom(g, r);
+  else         return; // g==0 => but do a p_Cleardenom(f)
+  if (f==NULL) return; // f==0 => but do a p_Cleardenom(g)
+
+  Off(SW_RATIONAL);
+  CanonicalForm F,G,GCD;
+  if (rField_is_Q(r) || (rField_is_Zp(r)))
+  {
+    bool b1=isOn(SW_USE_EZGCD_P);
+    Off (SW_USE_NTL_GCD_P);
+    setCharacteristic( rChar(r) );
+    F=convSingPFactoryP( f,r );
+    G=convSingPFactoryP( g,r );
+    GCD=gcd(F,G);
+    if (!GCD.isOne())
+    {
+      p_Delete(&f,r);
+      p_Delete(&g,r);
+      f=convFactoryPSingP( F/GCD, r);
+      g=convFactoryPSingP( G/GCD, r);
+    }
+    if (!b1) Off (SW_USE_EZGCD_P);
+  }
+  // and over Q(a) / Fp(a)
+  else if ( rField_is_Extension(r))
+  {
+    if ( rField_is_Q_a(r)) setCharacteristic( 0 );
+    else                   setCharacteristic( rChar(r) );
+    if (r->cf->extRing->qideal!=NULL)
+    {
+      bool b1=isOn(SW_USE_QGCD);
+      if ( rField_is_Q_a(r) ) On(SW_USE_QGCD);
+      CanonicalForm mipo=convSingPFactoryP(r->cf->extRing->qideal->m[0],
+                                           r->cf->extRing);
+      Variable a=rootOf(mipo);
+      F=( convSingAPFactoryAP( f,a,r ) );
+      G=( convSingAPFactoryAP( g,a,r ) );
+      GCD=gcd(F,G);
+      if (!GCD.isOne())
+      {
+        p_Delete(&f,r);
+        p_Delete(&g,r);
+        f= convFactoryAPSingAP( F/GCD,r );
+        g= convFactoryAPSingAP( G/GCD,r );
+      }
+      if (!b1) Off(SW_USE_QGCD);
+    }
+    else
+    {
+      F=( convSingTrPFactoryP( f,r ) );
+      G=( convSingTrPFactoryP( g,r ) );
+      GCD=gcd(F,G);
+      if (!GCD.isOne())
+      {
+        p_Delete(&f,r);
+        p_Delete(&g,r);
+        f= convFactoryPSingTrP( F/GCD,r );
+        g= convFactoryPSingTrP( G/GCD,r );
+      }
+    }
+  }
+  else
+    WerrorS( feNotImplemented );
+  Off(SW_RATIONAL);
+}
+
 poly singclap_gcd ( poly f, poly g, const ring r)
 {
   poly res=NULL;
 
-  if (f!=NULL) p_Cleardenom(f, r);
+  if (f!=NULL) p_Cleardenom(f, r);             
   if (g!=NULL) p_Cleardenom(g, r);
   else         return f; // g==0 => gcd=f (but do a p_Cleardenom)

libpolys/polys/clapsing.h

@@ -26 +26 @@
 poly singclap_gcd ( poly f, poly g, const ring r );
 
+/// clears denominators of f and g, divides by gcd(f,g)
+void singclap_gcd_and_divide ( poly& f, poly& g, const ring r);
 
 // commented out!