Changeset c74742 in git

Timestamp:

Jul 31, 2019, 2:24:45 PM (4 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', '0604212ebb110535022efecad887940825b97c3f')

Children:

5cf321599c682cc78b4e3d193a6942c4b89223b2

Parents:

6ca326dd0aa9608cd34d22903d3d6a579fde9c80

Message:

fix: content for gcd

File:

• factory/FLINTconvert.cc (modified) (2 diffs)

factory/FLINTconvert.cc

@@ -752 +752 @@
 }
 
+static CanonicalForm b_content ( const CanonicalForm & f )
+{
+    if ( f.inCoeffDomain() )
+        return f;
+    else
+    {
+        CanonicalForm result = 0;
+        CFIterator i;
+        for ( i = f; i.hasTerms() && (!result.isOne()); i++ )
+            result=bgcd( b_content(i.coeff()) , result );
+        return result;
+    }
+}
+
+
 CanonicalForm gcdFlintMP_QQ(const CanonicalForm& F, const CanonicalForm& G)
 {
@@ -781 +796 @@
     RES=convFlintMPFactoryP(res,ctx,N);
     // gcd(2x,4x) should be 2x, so RES should also have the gcd(lc(F),lc(G))
-    RES*=gcd(F.lc(),G.lc());
+    RES*=bgcd(b_content(F),b_content(G));
   }
   fmpq_mpoly_clear(res,ctx);