Changeset 22c7f0 in git for libpolys

Timestamp:

Jul 17, 2013, 6:21:49 PM (11 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', '17f1d200f27c5bd38f5dfc6e8a0879242279d1d8')

Children:

9713eea80d019296cc6a670a2d19b267b228d562

Parents:

891122532e023bc9e189b3dbd0e4d94bec28c798

Message:

fix: p_Content for Q_a

File:

• libpolys/polys/monomials/p_polys.cc (modified) (6 diffs)

libpolys/polys/monomials/p_polys.cc

@@ -2127 +2127 @@
       if (getCoeffType(r->cf)==n_algExt)
       {
-        number hzz = n_Init(1, r->cf->extRing->cf);
+        h = n_Init(1, r->cf->extRing->cf);
         p=ph;
         while (p!=NULL)
@@ -2135 +2135 @@
           while (c_n!=NULL)
           { // each monom: coeff in Q
-            d=n_Lcm(hzz,pGetCoeff(c_n),r->cf->extRing->cf);
-            n_Delete(&hzz,r->cf->extRing->cf);
-            hzz=d;
+            d=n_Lcm(h,pGetCoeff(c_n),r->cf->extRing->cf);
+            n_Delete(&h,r->cf->extRing->cf);
+            h=d;
             pIter(c_n);
           }
           pIter(p);
         }
-        /* hzz contains the 1/lcm of all denominators in c_n_n*/
-        h=n_Invers(hzz,r->cf->extRing->cf);
-        n_Delete(&hzz,r->cf->extRing->cf);
-        n_Normalize(h,r->cf->extRing->cf);
+        /* h contains the 1/lcm of all denominators in c_n_n*/
+        //n_Normalize(h,r->cf->extRing->cf);
         if(!n_IsOne(h,r->cf->extRing->cf))
         {
@@ -2504 +2502 @@
 #endif
 
-  poly start=p;
 
   number d, h;
@@ -2514 +2511 @@
     assume( n_GreaterZero(pGetCoeff(p),C) );
     if(!n_GreaterZero(pGetCoeff(p),C)) p = p_Neg(p,r);
-    return start;
+    return p;
   }
 #endif
@@ -2521 +2518 @@
   {
     if(!n_GreaterZero(pGetCoeff(p),C)) p = p_Neg(p,r);
-    return start;
+    return p;
   }
 
@@ -2547 +2544 @@
     if(!n_GreaterZero(pGetCoeff(p),C)) p = p_Neg(p,r);
 
-    return start;
+    return p;
   }
 
   assume(pNext(p)!=NULL);
+  poly start=p;
 
 #if 0 && CLEARENUMERATORS