Changeset ac2c30 in git


Ignore:
Timestamp:
Mar 13, 2014, 4:26:09 PM (10 years ago)
Author:
Martin Lee <martinlee84@…>
Branches:
(u'spielwiese', '17f1d200f27c5bd38f5dfc6e8a0879242279d1d8')
Children:
36aaf82c1afa5b8a7881749075e68ae881d2c3e7
Parents:
ea0a9dfbfc8e8cda45b109756bdf6380f6d26c3d
git-author:
Martin Lee <martinlee84@web.de>2014-03-13 16:26:09+01:00
git-committer:
Hans Schoenemann <hannes@mathematik.uni-kl.de>2014-03-17 18:31:55+01:00
Message:
chg: improve QuasiInverse
File:
1 edited

Legend:

Unmodified
Added
Removed
  • factory/libfac/charset/alg_factor.cc

    rea0a9d rac2c30  
    470470CanonicalForm
    471471QuasiInverse (const CanonicalForm& f, const CanonicalForm& g,
    472       CanonicalForm& numt, const Variable& x)
     472              const Variable& x)
    473473{
    474474  CanonicalForm pi, pi1, q, t0, t1, Hi, bi, pi2;
     
    514514  if (!isRat)
    515515    Off (SW_RATIONAL);
    516   numt= t1;
    517   return pi;
     516  return t1;
    518517}
    519518
     
    645644        CanonicalForm hh= replacevar (h, oldR.mvar(), alpha);
    646645
    647         CanonicalForm numt, dent;
    648         QuasiInverse (tmp.getFirst(), LC (h), numt, tmp.getFirst().mvar());
     646        CanonicalForm numinv, deninv;
     647        numinv= QuasiInverse (tmp.getFirst(), LC (h), tmp.getFirst().mvar());
    649648
    650649        Off (SW_RATIONAL);
    651         h *= numt;
     650        h *= numinv;
    652651        h= reduce (h, tmp.getFirst());
    653         dent= LC(h);
     652        deninv= LC(h);
    654653
    655654        ra= -h[0];
    656         denra= gcd (ra, dent);
     655        denra= gcd (ra, deninv);
    657656        ra /= denra;
    658         denra= dent/denra;
     657        denra= deninv/denra;
    659658        denra= replacevar (denra, ra.mvar(), g.mvar());
    660659        ra= replacevar(ra, ra.mvar(), g.mvar());
     
    800799  // make quasi monic
    801800  CFList Rstarlist= CFList (Rstar);
    802   CanonicalForm numt, dent;
     801  CanonicalForm numinv;
    803802  On (SW_RATIONAL);
    804   QuasiInverse (Rstar, LC(f), numt, Rstar.mvar());
    805 
    806   f *= numt;
     803  numinv= QuasiInverse (Rstar, LC(f), Rstar.mvar());
     804
     805  f *= numinv;
    807806  f= Prem (f, Rstarlist);
    808807  f /= vcontent (f, Rstar.mvar());
     
    858857
    859858      h= alg_gcd (g, fnew, Rstarlist);
    860       QuasiInverse(Rstar, LC(h), numt, Rstar.mvar());
    861       dent= 1;
    862       h *= numt;
     859      numinv= QuasiInverse(Rstar, LC(h), Rstar.mvar());
     860      h *= numinv;
    863861      h= Prem (h, Rstarlist);
    864862      h /= vcontent (h, Rstar.mvar());
Note: See TracChangeset for help on using the changeset viewer.