Changeset ac2c30 in git
 Timestamp:
 Mar 13, 2014, 4:26:09 PM (10 years ago)
 Branches:
 (u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')
 Children:
 36aaf82c1afa5b8a7881749075e68ae881d2c3e7
 Parents:
 ea0a9dfbfc8e8cda45b109756bdf6380f6d26c3d
 gitauthor:
 Martin Lee <martinlee84@web.de>20140313 16:26:09+01:00
 gitcommitter:
 Hans Schoenemann <hannes@mathematik.unikl.de>20140317 18:31:55+01:00
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

factory/libfac/charset/alg_factor.cc
rea0a9d rac2c30 470 470 CanonicalForm 471 471 QuasiInverse (const CanonicalForm& f, const CanonicalForm& g, 472 CanonicalForm& numt,const Variable& x)472 const Variable& x) 473 473 { 474 474 CanonicalForm pi, pi1, q, t0, t1, Hi, bi, pi2; … … 514 514 if (!isRat) 515 515 Off (SW_RATIONAL); 516 numt= t1; 517 return pi; 516 return t1; 518 517 } 519 518 … … 645 644 CanonicalForm hh= replacevar (h, oldR.mvar(), alpha); 646 645 647 CanonicalForm num t, 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()); 649 648 650 649 Off (SW_RATIONAL); 651 h *= num t;650 h *= numinv; 652 651 h= reduce (h, tmp.getFirst()); 653 den t= LC(h);652 deninv= LC(h); 654 653 655 654 ra= h[0]; 656 denra= gcd (ra, den t);655 denra= gcd (ra, deninv); 657 656 ra /= denra; 658 denra= den t/denra;657 denra= deninv/denra; 659 658 denra= replacevar (denra, ra.mvar(), g.mvar()); 660 659 ra= replacevar(ra, ra.mvar(), g.mvar()); … … 800 799 // make quasi monic 801 800 CFList Rstarlist= CFList (Rstar); 802 CanonicalForm num t, dent;801 CanonicalForm numinv; 803 802 On (SW_RATIONAL); 804 QuasiInverse (Rstar, LC(f), numt, Rstar.mvar());805 806 f *= num t;803 numinv= QuasiInverse (Rstar, LC(f), Rstar.mvar()); 804 805 f *= numinv; 807 806 f= Prem (f, Rstarlist); 808 807 f /= vcontent (f, Rstar.mvar()); … … 858 857 859 858 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; 863 861 h= Prem (h, Rstarlist); 864 862 h /= vcontent (h, Rstar.mvar());
Note: See TracChangeset
for help on using the changeset viewer.