Changeset 4609b7 in git for libpolys/polys

Timestamp:

Oct 7, 2013, 3:56:41 PM (11 years ago)

Author:

Martin Lee <martinlee84@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

0b4a87ab19d7bd751e10156552aa4184fbff4ea7

Parents:

9030ce67e86f3969436b1e54affb5d969d7143bc

git-author:

Martin Lee <martinlee84@web.de>2013-10-07 15:56:41+02:00

git-committer:

Martin Lee <martinlee84@web.de>2013-10-07 18:15:10+02:00

Message:

chg: extra treatment of transcendental extensions in p_Content

File:

• libpolys/polys/monomials/p_polys.cc (modified) (1 diff)

libpolys/polys/monomials/p_polys.cc

@@ -2212 +2212 @@
         n_Delete(&h,r->cf->extRing->cf);
       }*/
+    }
+
+    /* normalize ph over a transcendental extension s.t.
+       lead (ph) is > 0 if extRing->cf == Q
+       or lead (ph) is monic if extRing->cf == Zp*/
+    const coeffs C = r->cf;
+    if (nCoeff_is_transExt(C))
+    {
+      p= ph;
+      h= p_GetCoeff (p, C);
+      fraction f = (fraction) h;
+      poly g = NUM(f);
+      number n=p_GetCoeff (g,C->extRing->cf);
+      if (rField_is_Q (C->extRing))
+      {
+        if (!n_GreaterZero(n,C->extRing->cf))
+        {
+          p=p_Neg (p,r);
+        }
+      }
+      else if (rField_is_Zp(C->extRing))
+      {
+        if (!n_IsOne (n, C->extRing->cf))
+        {
+          n=n_Invers (n,C->extRing->cf);
+          nMapFunc nMap;
+          nMap= n_SetMap (C->extRing->cf, C);
+          number ninv= nMap (n,C->extRing->cf, C);
+          p=p_Mult_nn (p, ninv, r);
+          n_Delete (&ninv, C);
+          n_Delete (&n, C->extRing->cf);
+        }
+      }
+      p= ph;
     }
   }