Changeset cc07066 in git for Tst/Short/pfd.tst


Ignore:
Timestamp:
Aug 31, 2020, 11:57:27 AM (4 years ago)
Author:
Janko Boehm <bohem@…>
Branches:
(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
Children:
fdea544a15d2eb16f952ba088de11e5efa5494d9
Parents:
6924d452049dc8bfbb851672c4b3ceac05081b6c
Message:
update
File:
1 edited

Legend:

Unmodified
Added
Removed
  • Tst/Short/pfd.tst

    r6924d4 rcc07066  
    66example pfd;
    77
    8 /* pfd(poly,poly) */
     8// pfd(poly,poly) //////////////////////////////////////////////////////////////
    99ring r1 = 0, x(1..5), dp;
    1010poly f = -2*x(1)*x(3)+3*x(3)*x(4)+x(2)*x(5)+x(3)*x(5)-x(4)*x(5);
     
    3232kill dec;
    3333
     34
     35// pfd(poly,list) (denominator = ideal of factors & intvec of exponents) ///////
    3436ring 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);
     37poly f = x+y+z+1;
     38list g = list(ideal((x^2+y^2+z^2),(x+y^2),(y+z^2),(z+x^2)), intvec(2,1,1,1));
     39list dec = pfd(f,g);
     40displaypfd(dec);
     41checkpfd(list(f,g), dec);
     42checkpfd(list(f,g), dec, 10);
     43kill dec;
     44
     45// different ordering, same polynomials:
     46ring r4 = 5, (x,y,z), lp;
     47poly f = fetch(r3,f);
     48list g = fetch(r3,g);
    3749list dec = pfd(f,g);
    3850displaypfd(dec);
     
    4153kill dec;
    4254
    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;
    5155
    52 /* pfd(list) */
     56// pfd(list) ///////////////////////////////////////////////////////////////////
    5357ring s1 = 0, (x,y,z), dp;
    5458poly f1 = x*y+y*z+z*x-x-y-z+1;
     
    6872kill dec;
    6973
    70 /* pfd(matrix) */
     74// pfd(matrix) /////////////////////////////////////////////////////////////////
    7175ring s2 = 3, (x,y,z), dp;
    7276poly f11 = (x+y+z+1)^3;
     
    8589displaypfd(dec[2][1]);
    8690displaypfd(dec[2][2]);
    87 checkpfd(list(f11,g11), dec[1][1]);
    88 checkpfd(list(f12,g12), dec[1][2]);
    89 checkpfd(list(f21,g21), dec[2][1]);
    90 checkpfd(list(f22,g22), dec[2][2]);
     91checkpfd(list(f1,g1), dec[1][1]);
     92checkpfd(list(f2,g2), dec[1][2]);
     93checkpfd(list(f3,g3), dec[2][1]);
     94checkpfd(list(f4,g4), dec[2][2]);
    9195kill dec;
    9296
Note: See TracChangeset for help on using the changeset viewer.