[b1ac99] | 1 | LIB "tst.lib"; |
---|
| 2 | tst_init(); |
---|
| 3 | option(prot); |
---|
| 4 | //option(notRegularity); |
---|
| 5 | //ring r=31991,(t,x,y,z,w),(c,ls); |
---|
| 6 | //ideal i=x2+y2yz,y2+z2zw,t2+x2xy,x3+y3+t3,xyz; |
---|
| 7 | //ideal ii=x3+yzw,y3+zwt,z3+wtx,w3+txy,t2+xy; |
---|
| 8 | //ideal iii=x3+yzw,y3+zwt,z3+wtx,w2+txy,t2+xy; |
---|
| 9 | //ideal tt=t2x2+tx2y+x2yz,t2y2+ty2z+y2zw,t2z2+tz2w+xz2w,t2w2+txw2+xyw2; |
---|
| 10 | //test(17); |
---|
| 11 | //ideal T=std(tt); |
---|
| 12 | //timer = 1; |
---|
| 13 | //res(T,0,t1); |
---|
| 14 | //sres(T,0,t2); |
---|
| 15 | //mres(T,0,t3); |
---|
| 16 | //-------------------------------------------------------------- |
---|
| 17 | ring an=32003,(w,x,y,z),(dp,C); |
---|
| 18 | ideal i= |
---|
| 19 | 1w2xy+1w2xz+1w2yz+1wxyz+1x2yz+1xy2z+1xyz2, |
---|
| 20 | 1w4x+1w4z+1w3yz+1w2xyz+1wx2yz+1x2y2z+1xy2z2, |
---|
| 21 | 1w6+1w5z+1w4xz+1w3xyz+1w2xy2z+1wx2y2z+1x2y2z2; |
---|
| 22 | ideal j=std(i); |
---|
| 23 | list jres=sres(j,0); |
---|
| 24 | list jmin=minres(jres); |
---|
| 25 | print(betti(jmin),"betti"); |
---|
| 26 | list jres1=mres(i,0); |
---|
| 27 | print(betti(jres1),"betti"); |
---|
| 28 | list jres2=res(i,0); |
---|
| 29 | list jmin2=minres(jres2); |
---|
| 30 | print(betti(jmin2),"betti"); |
---|
| 31 | kill an; |
---|
| 32 | //-------------------------------------------------------------- |
---|
| 33 | ring pa2=32003,(w,t,u,x,y,z),(dp,C); |
---|
| 34 | ideal i= |
---|
| 35 | 1w10t2u2-2w5tu2x6-2w5tu2y6+2w5tu2z6+1u2x12+2u2x6y6-2u2x6z6+1u2y12-2u2y6z6+1u2z12, |
---|
| 36 | -1w5t5u-1w5tux4-1w5tuy4-1w5tuz4+1t4ux6+1t4uy6-1t4uz6+1ux10+1ux6y4+1ux6z4 |
---|
| 37 | +1ux4y6-1ux4z6+1uy10+1uy6z4-1uy4z6-1uz10, |
---|
| 38 | 1t8+2t4x4+2t4y4+2t4z4+1x8+2x4y4+2x4z4+1y8+2y4z4+1z8, |
---|
| 39 | 1w5t3u-1w5tux2-1w5tuy2+1w5tuz2-1t2ux6-1t2uy6+1t2uz6+1ux8+1ux6y2-1ux6z2 |
---|
| 40 | +1ux2y6-1ux2z6+1uy8-1uy6z2-1uy2z6+1uz8, |
---|
| 41 | -1t6+1t4x2+1t4y2-1t4z2-1t2x4-1t2y4-1t2z4+1x6+1x4y2-1x4z2+1x2y4+1x2z4+1y6 |
---|
| 42 | -1y4z2+1y2z4-1z6, |
---|
| 43 | 1t4-2t2x2-2t2y2+2t2z2+1x4+2x2y2-2x2z2+1y4-2y2z2+1z4; |
---|
| 44 | ideal j=std(i); |
---|
| 45 | list jres=sres(j,0); |
---|
| 46 | list jmin=minres(jres); |
---|
| 47 | print(betti(jmin),"betti"); |
---|
| 48 | list jres1=mres(i,0); |
---|
| 49 | print(betti(jres1),"betti"); |
---|
| 50 | list jres2=res(i,0); |
---|
| 51 | list jmin2=minres(jres2); |
---|
| 52 | print(betti(jmin2),"betti"); |
---|
| 53 | kill pa2; |
---|
| 54 | //-------------------------------------------------------------- |
---|
| 55 | ring sw5=32003,(a,b,c,d,e,h),(dp,C); |
---|
| 56 | ideal i= |
---|
| 57 | -1ab-1b2-2de-2ch, |
---|
| 58 | -1ac-2bc-1e2-2dh, |
---|
| 59 | -1c2-1ad-2bd-2eh, |
---|
| 60 | -2cd-1ae-2be-1h2, |
---|
| 61 | -1d2-2ce-1ah-2bh; |
---|
| 62 | ideal j=std(i); |
---|
| 63 | list jres=sres(j,0); |
---|
| 64 | list jmin=minres(jres); |
---|
| 65 | print(betti(jmin),"betti"); |
---|
| 66 | list jres1=mres(i,0); |
---|
| 67 | print(betti(jres1),"betti"); |
---|
| 68 | list jres2=res(i,0); |
---|
| 69 | list jmin2=minres(jres2); |
---|
| 70 | print(betti(jmin2),"betti"); |
---|
| 71 | kill sw5; |
---|
| 72 | //-------------------------------------------------------------- |
---|
| 73 | ring r=32003,(t,x,y,z,w),(dp,C); |
---|
| 74 | ideal i= |
---|
| 75 | x2w+y2z, |
---|
| 76 | y2x+z2w, |
---|
| 77 | t2w+x2y, |
---|
| 78 | x3+y3+t3; |
---|
| 79 | ideal j=std(i); |
---|
| 80 | list jres=sres(j,0); |
---|
| 81 | list jmin=minres(jres); |
---|
| 82 | print(betti(jmin),"betti"); |
---|
| 83 | list jres1=mres(i,0); |
---|
| 84 | print(betti(jres1),"betti"); |
---|
| 85 | list jres2=res(i,0); |
---|
| 86 | list jmin2=minres(jres2); |
---|
| 87 | print(betti(jmin2),"betti"); |
---|
| 88 | kill r; |
---|
| 89 | //-------------------------------------------------------------- |
---|
| 90 | |
---|
| 91 | $ |
---|
| 92 | |
---|