Changeset cc07066 in git for Tst/Short/pfd.tst
- Timestamp:
- Aug 31, 2020, 11:57:27 AM (4 years ago)
- Branches:
- (u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
- Children:
- fdea544a15d2eb16f952ba088de11e5efa5494d9
- Parents:
- 6924d452049dc8bfbb851672c4b3ceac05081b6c
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Tst/Short/pfd.tst
r6924d4 rcc07066 6 6 example pfd; 7 7 8 / * pfd(poly,poly) */8 // pfd(poly,poly) ////////////////////////////////////////////////////////////// 9 9 ring r1 = 0, x(1..5), dp; 10 10 poly f = -2*x(1)*x(3)+3*x(3)*x(4)+x(2)*x(5)+x(3)*x(5)-x(4)*x(5); … … 32 32 kill dec; 33 33 34 35 // pfd(poly,list) (denominator = ideal of factors & intvec of exponents) /////// 34 36 ring r3 = 5, (x,y,z), dp; 35 poly f = 2*x^3-x^2*y+x*y^2+y^3-2*x^2*z-2*x*y*z-y^2*z+2*x*z^2+y*z^2-2*z^3; 36 poly g = (x^2+3*y)^2*(x-2*y^2)*(x+y+1)^3*(2*x+3*y+4)*(x^2-x*y+y^2); 37 poly f = x+y+z+1; 38 list g = list(ideal((x^2+y^2+z^2),(x+y^2),(y+z^2),(z+x^2)), intvec(2,1,1,1)); 39 list dec = pfd(f,g); 40 displaypfd(dec); 41 checkpfd(list(f,g), dec); 42 checkpfd(list(f,g), dec, 10); 43 kill dec; 44 45 // different ordering, same polynomials: 46 ring r4 = 5, (x,y,z), lp; 47 poly f = fetch(r3,f); 48 list g = fetch(r3,g); 37 49 list dec = pfd(f,g); 38 50 displaypfd(dec); … … 41 53 kill dec; 42 54 43 ring r4 = 5, (x,y,z), lp;44 poly f = fetch(r3,f);45 poly g = fetch(r3,g);46 list dec = pfd(f,g);47 displaypfd(dec);48 checkpfd(list(f,g), dec);49 checkpfd(list(f,g), dec,10);50 kill dec;51 55 52 / * pfd(list) */56 // pfd(list) /////////////////////////////////////////////////////////////////// 53 57 ring s1 = 0, (x,y,z), dp; 54 58 poly f1 = x*y+y*z+z*x-x-y-z+1; … … 68 72 kill dec; 69 73 70 / * pfd(matrix) */74 // pfd(matrix) ///////////////////////////////////////////////////////////////// 71 75 ring s2 = 3, (x,y,z), dp; 72 76 poly f11 = (x+y+z+1)^3; … … 85 89 displaypfd(dec[2][1]); 86 90 displaypfd(dec[2][2]); 87 checkpfd(list(f1 1,g11), dec[1][1]);88 checkpfd(list(f 12,g12), dec[1][2]);89 checkpfd(list(f 21,g21), dec[2][1]);90 checkpfd(list(f 22,g22), dec[2][2]);91 checkpfd(list(f1,g1), dec[1][1]); 92 checkpfd(list(f2,g2), dec[1][2]); 93 checkpfd(list(f3,g3), dec[2][1]); 94 checkpfd(list(f4,g4), dec[2][2]); 91 95 kill dec; 92 96
Note: See TracChangeset
for help on using the changeset viewer.