Changeset 3d5596 in git for libpolys

Timestamp:

Jul 11, 2017, 9:59:44 AM (7 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', '2a584933abf2a2d3082034c7586d38bb6de1a30a')

Children:

0ec309ddd6253287d647226c352fff30a3417b99

Parents:

3092f65270fb472dca33d7d41d19f79dbb1cd638

Message:

add: remainder of polynomial division via factoryi(mod,%), also div

Location:

libpolys/polys

Files:

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

libpolys/polys/clapsing.cc

@@ -575 +575 @@
 }
 
+poly singclap_pmod ( poly f, poly g, const ring r )
+{
+  poly res=NULL;
+  On(SW_RATIONAL);
+  if (rField_is_Zp(r) || rField_is_Q(r))
+  {
+    setCharacteristic( rChar(r) );
+    CanonicalForm F( convSingPFactoryP( f,r ) ), G( convSingPFactoryP( g,r ) );
+    res = convFactoryPSingP( F % G,r );
+  }
+  else if (rField_is_Ring_Z(r))
+  {
+    Off(SW_RATIONAL);
+    setCharacteristic( rChar(r) );
+    CanonicalForm F( convSingPFactoryP( f,r ) ), G( convSingPFactoryP( g,r ) );
+    res = convFactoryPSingP( F % G,r );
+  }
+  else if (r->cf->extRing!=NULL)
+  {
+    if (rField_is_Q_a(r)) setCharacteristic( 0 );
+    else               setCharacteristic( rChar(r) );
+    if (r->cf->extRing->qideal!=NULL)
+    {
+      CanonicalForm mipo=convSingPFactoryP(r->cf->extRing->qideal->m[0],
+                                                 r->cf->extRing);
+      Variable a=rootOf(mipo);
+      CanonicalForm F( convSingAPFactoryAP( f,a,r ) ),
+                    G( convSingAPFactoryAP( g,a,r ) );
+      res= convFactoryAPSingAP(  F % G, r  );
+      prune (a);
+    }
+    else
+    {
+      CanonicalForm F( convSingTrPFactoryP( f,r ) ), G( convSingTrPFactoryP( g,r ) );
+      res= convFactoryPSingTrP(  F % G,r  );
+    }
+  }
+#if 0 // not yet working
+  else if (rField_is_GF())
+  {
+    //Print("GF(%d^%d)\n",nfCharP,nfMinPoly[0]);
+    setCharacteristic( nfCharP,nfMinPoly[0], currRing->parameter[0][0] );
+    CanonicalForm F( convSingGFFactoryGF( f ) ), G( convSingGFFactoryGF( g ) );
+    res = convFactoryGFSingGF( F / G );
+  }
+#endif
+  else
+    WerrorS( feNotImplemented );
+  Off(SW_RATIONAL);
+  return res;
+}
+
 void singclap_divide_content ( poly f, const ring r )
 {

libpolys/polys/clapsing.h

@@ -39 +39 @@
 
 poly singclap_pdivide ( poly f, poly g, const ring r );
+poly singclap_pmod ( poly f, poly g, const ring r );
 
 void singclap_divide_content ( poly f, const ring r);