Changeset e2dd6e in git for libpolys/polys/monomials/p_polys.cc

Timestamp:

Apr 30, 2013, 11:18:19 PM (11 years ago)

Author:

Mathias Schulze <mschulze@…>

Branches:

(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', '38dfc5131670d387a89455159ed1e071997eec94')

Children:

b79fd76a535a05f415a6fd8c00ffbff7f2d2a32d

Parents:

232a34c254032045d3fbc55b3c260a4044ccd608f1d88ec2c56f693af48faa279e80b42e6d11b48b

Message:

Merge branch 'spielwiese' of github.com:mschulze/Sources into spielwiese

File:

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

libpolys/polys/monomials/p_polys.cc

@@ -2120 +2120 @@
     }
     n_Delete(&h,r->cf);
-#ifdef HAVE_FACTORY
-//    if ( (n_GetChar(r) == 1) || (n_GetChar(r) < 0) ) /* Q[a],Q(a),Zp[a],Z/p(a) */
-//    {
-//      singclap_divide_content(ph, r);
-//      if(!n_GreaterZero(pGetCoeff(ph),r->cf)) ph = p_Neg(ph,r);
-//    }
-#endif
     if (rField_is_Q_a(r))
     {
-      // we only need special handling for alg. ext.
+      // special handling for alg. ext.:
       if (getCoeffType(r->cf)==n_algExt)
       {
@@ -2166 +2159 @@
             }
             pIter(p);
+          }
+        }
+        n_Delete(&h,r->cf->extRing->cf);
+      }
+      else
+      {
+      // special handling for rat. functions.:
+        number hzz =NULL;
+        p=ph;
+        while (p!=NULL)
+        { // each monom: coeff in Q_a
+          fraction f=(fraction)pGetCoeff(p); 
+          poly c_n=NUM(f);
+          if (hzz==NULL)
+          {
+            hzz=n_Copy(pGetCoeff(c_n),r->cf->extRing->cf); 
+            pIter(c_n);
+          }
+          while ((c_n!=NULL)&&(!n_IsOne(hzz,r->cf->extRing->cf)))
+          { // each monom: coeff in Q
+            d=n_Gcd(hzz,pGetCoeff(c_n),r->cf->extRing->cf);
+            n_Delete(&hzz,r->cf->extRing->cf);
+            hzz=d;
+            pIter(c_n);
+          }
+          pIter(p);
+        }
+        /* hzz contains the gcd of all numerators in f*/
+        h=n_Invers(hzz,r->cf->extRing->cf);
+        n_Delete(&hzz,r->cf->extRing->cf);
+        n_Normalize(h,r->cf->extRing->cf);
+        if(!n_IsOne(h,r->cf->extRing->cf))
+        {
+          p=ph;
+          while (p!=NULL)
+          { // each monom: coeff in Q_a
+            fraction f=(fraction)pGetCoeff(p);
+            NUM(f)=p_Mult_nn(NUM(f),h,r->cf->extRing);
+            p_Normalize(NUM(f),r->cf->extRing);
+            pIter(p);
           }
         }