Changeset ac2c30 in git
- Timestamp:
- Mar 13, 2014, 4:26:09 PM (10 years ago)
- Branches:
- (u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')
- 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
- 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.