Changeset 7600fb1 in git for libpolys/polys

Timestamp:

Feb 21, 2013, 2:40:12 PM (11 years ago)

Author:

Martin Lee <martinlee84@…>

Branches:

(u'spielwiese', '17f1d200f27c5bd38f5dfc6e8a0879242279d1d8')

Children:

680ca3e411efc89e2a5f3806c29b9bec84481860

Parents:

5523765c597320ca80eb947391587ecf21dd1539

git-author:

Martin Lee <martinlee84@web.de>2013-02-21 14:40:12+01:00

git-committer:

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

Message:

fix: also map denominators

File:

• libpolys/polys/ext_fields/algext.cc (modified) (2 diffs)

libpolys/polys/ext_fields/algext.cc

@@ -926 +926 @@
   assume (nCoeff_is_algExt (dst));
   fraction fa=(fraction)a;
-  if (!DENIS1 (fa))
-    WarnS ("ignoring denominators in naCopyTrans2AlgExt");
-  poly p;
+  poly p, q;
   if (rSamePolyRep(src->extRing, dst->extRing))
   {
     p = p_Copy(NUM(fa),src->extRing);
+    if (!DENIS1(fa))
+    {
+      q = p_Copy(DEN(fa),src->extRing);
+      assume (q != NULL);
+    }
   }
   else
@@ -941 +944 @@
     assume (nMap != NULL);
     p= p_PermPoly (NUM (fa), NULL, src->extRing, dst->extRing,nMap, NULL,rVar (src->extRing));
+    if (!DENIS1(fa))
+    {
+      q= p_PermPoly (DEN (fa), NULL, src->extRing, dst->extRing,nMap, NULL,rVar (src->extRing));
+      assume (q != NULL);
+    }
   }
   definiteReduce(p, dst->extRing->qideal->m[0], dst);
   assume (p_Test (p, dst->extRing));
-  return (number)p;
+  if (!DENIS1(fa))
+  {
+    definiteReduce(q, dst->extRing->qideal->m[0], dst);
+    assume (p_Test (q, dst->extRing));
+    if (q != NULL)
+    {
+      number t= naDiv ((number)p,(number)q, dst);
+      p_Delete (&p, dst->extRing);
+      p_Delete (&q, dst->extRing);
+      return t;
+    }
+    WerrorS ("mapping denominator to zero");
+  }
+  return (number) p;
 }