| 1 | LIB "primdec.lib"; |
|---|
| 2 | ring r = 32003,(x,y,z),dp; |
|---|
| 3 | poly p = z2+1; |
|---|
| 4 | poly q = z4+2; |
|---|
| 5 | ideal i = p^2*q^3,(y-z3)^3,(x-yz+z4)^4; |
|---|
| 6 | list pr= prepareAss(i); |
|---|
| 7 | pr; |
|---|
| 8 | |
|---|
| 9 | //////////////////////////////////////////////////////////////// |
|---|
| 10 | // test that one more 'is no standard basis' warning is gone away |
|---|
| 11 | system("random",1452807367); |
|---|
| 12 | proc primdecMinpolyExample() |
|---|
| 13 | { |
|---|
| 14 | ring rng = (7,vv),(x,y,z,u),dp; |
|---|
| 15 | minpoly = vv^2+vv+3; |
|---|
| 16 | ideal I = (-3*vv-1)*u^2+2,2*x*z+(2*vv-2)*x+(2*vv+1),(vv+2)*x^2*u+(2*vv+1)*z^2*u+(vv)*z*u^2; |
|---|
| 17 | list result = minAssGTZ (I); |
|---|
| 18 | } |
|---|
| 19 | |
|---|
| 20 | // increase call depth: |
|---|
| 21 | proc testPrimdecMinpolyExample() { return( primdecMinpolyExample() ); } |
|---|
| 22 | testPrimdecMinpolyExample(); |
|---|
| 23 | |
|---|
| 24 | // check that redefinig gnir1 is gone away: |
|---|
| 25 | system("random",1463736066); |
|---|
| 26 | ring rng2 = (0),(x,y,z,u,w),dp; |
|---|
| 27 | ideal I2 = -28*1*x*z*w*1-15*y*z*u*w,-28*1^2*x*u+14*1*x^2*y*w*1+17,12*1*u*w^2*1+18*x*y^2*w; |
|---|
| 28 | list result2 = primdecGTZ (I2); |
|---|
| 29 | |
|---|
| 30 | |
|---|
| 31 | LIB "tst.lib"; tst_status(1); |
|---|
| 32 | $ |
|---|
