Changeset c7c7fe4 in git

Timestamp:

Feb 16, 2012, 3:19:10 PM (11 years ago)

Author:

Martin Lee <martinlee84@…>

Branches:

(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'a800fe4b3e9d37a38c5a10cc0ae9dfa0c15a4ee6')

Children:

5b410740267969d96b217ae7a6dc01d921a86161

Parents:

e016ba74132ef655a25823cf0f343ce4b9069f7d

git-author:

Martin Lee <martinlee84@web.de>2012-02-16 15:19:10+01:00

git-committer:

Martin Lee <martinlee84@web.de>2012-04-04 14:42:25+02:00

Message:

chg: added divisibility test for univariate polys

Location:

factory

Files:

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

factory/facMul.cc

@@ -2497 +2497 @@
 }
 
+bool
+uniFdivides (const CanonicalForm& A, const CanonicalForm& B)
+{
+  int p= getCharacteristic();
+  if (p > 0)
+  {
+    zz_p::init (p);
+    Variable alpha;
+    if (hasFirstAlgVar (A, alpha) || hasFirstAlgVar (B, alpha))
+    {
+      zz_pX NTLMipo= convertFacCF2NTLzzpX (getMipo (alpha));
+      zz_pE::init (NTLMipo);
+      zz_pEX NTLA= convertFacCF2NTLzz_pEX (A, NTLMipo);
+      zz_pEX NTLB= convertFacCF2NTLzz_pEX (B, NTLMipo);
+      return divide (NTLB, NTLA);
+    }
+#ifdef HAVE_FLINT
+    nmod_poly_t FLINTA, FLINTB;
+    convertFacCF2nmod_poly_t (FLINTA, A);
+    convertFacCF2nmod_poly_t (FLINTB, B);
+    nmod_poly_rem (FLINTA, FLINTB, FLINTA);
+    bool result= nmod_poly_is_zero (FLINTA);
+    nmod_poly_clear (FLINTA);
+    nmod_poly_clear (FLINTB);
+    return result;
+#else
+    zz_pX NTLA= convertFacCF2NTLzzpX (A);
+    zz_pX NTLB= convertFacCF2NTLzzpX (B);
+    return divide (NTLB, NTLA);
+#endif
+  }
+#ifdef HAVE_FLINT
+  fmpq_poly_t FLINTA,FLINTB;
+  fmpq_poly_init (FLINTA);
+  fmpq_poly_init (FLINTB);
+  convertFacCF2Fmpq_poly_t (FLINTA, A);
+  convertFacCF2Fmpq_poly_t (FLINTB, B);
+  fmpq_poly_rem (FLINTA, FLINTB, FLINTA);
+  bool result= fmpq_poly_is_zero (FLINTA);
+  fmpq_poly_clear (FLINTA);
+  fmpq_poly_clear (FLINTB);
+  return result;
+#else
+  return fdivides (A, B); //maybe NTL?
+#endif
+}
+
 // end division
 

factory/facMul.h

@@ -99 +99 @@
           );
 
+/// divisibility test for univariate polys
+///
+/// @return @a uniFdivides returns true if A divides B
+bool
+uniFdivides (const CanonicalForm& A, ///< [in] univariate poly
+             const CanonicalForm& B  ///< [in] univariate poly
+            );
+
 /// Karatsuba style modular multiplication for bivariate polynomials.
 ///