Changeset 7600fb1 in git


Ignore:
Timestamp:
Feb 21, 2013, 2:40:12 PM (10 years ago)
Author:
Martin Lee <martinlee84@…>
Branches:
(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'cbf866d4f8c3861325311c042cb2033c053453b2')
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:
1 edited

Legend:

Unmodified
Added
Removed
  • libpolys/polys/ext_fields/algext.cc

    r5523765 r7600fb1  
    926926  assume (nCoeff_is_algExt (dst));
    927927  fraction fa=(fraction)a;
    928   if (!DENIS1 (fa))
    929     WarnS ("ignoring denominators in naCopyTrans2AlgExt");
    930   poly p;
     928  poly p, q;
    931929  if (rSamePolyRep(src->extRing, dst->extRing))
    932930  {
    933931    p = p_Copy(NUM(fa),src->extRing);
     932    if (!DENIS1(fa))
     933    {
     934      q = p_Copy(DEN(fa),src->extRing);
     935      assume (q != NULL);
     936    }
    934937  }
    935938  else
     
    941944    assume (nMap != NULL);
    942945    p= p_PermPoly (NUM (fa), NULL, src->extRing, dst->extRing,nMap, NULL,rVar (src->extRing));
     946    if (!DENIS1(fa))
     947    {
     948      q= p_PermPoly (DEN (fa), NULL, src->extRing, dst->extRing,nMap, NULL,rVar (src->extRing));
     949      assume (q != NULL);
     950    }
    943951  }
    944952  definiteReduce(p, dst->extRing->qideal->m[0], dst);
    945953  assume (p_Test (p, dst->extRing));
    946   return (number)p;
     954  if (!DENIS1(fa))
     955  {
     956    definiteReduce(q, dst->extRing->qideal->m[0], dst);
     957    assume (p_Test (q, dst->extRing));
     958    if (q != NULL)
     959    {
     960      number t= naDiv ((number)p,(number)q, dst);
     961      p_Delete (&p, dst->extRing);
     962      p_Delete (&q, dst->extRing);
     963      return t;
     964    }
     965    WerrorS ("mapping denominator to zero");
     966  }
     967  return (number) p;
    947968}
    948969
Note: See TracChangeset for help on using the changeset viewer.