| 1 | //124. Beispiel mit Produktordnung, wo Macaulay ca 10 mal schneller ist: |
|---|
| 2 | // ca 43 sec auf PB 520c |
|---|
| 3 | // Macaulay hat 2 sec (ohne set autodegree, autocalc) |
|---|
| 4 | //option(mem); |
|---|
| 5 | //option(prot); |
|---|
| 6 | |
|---|
| 7 | ring P6 = 31991,(s,t,x,y,a,b,c,d,e,f,g),dp; |
|---|
| 8 | ideal i= -6s5x+a, -5s4tx-1s5y+b, -4s3t2x-2s4ty+c, |
|---|
| 9 | -3s2t3x-3s3t2y+d, -2st4x-4s2t3y+e, -1t5x-5st4y+f, -6t5y+g; |
|---|
| 10 | |
|---|
| 11 | ;ideal j= eliminate(i,stxy); |
|---|
| 12 | j; |
|---|
| 13 | |
|---|
| 14 | // andere Ringordnung: |
|---|
| 15 | ring Pa6 = 31991,(s,t,x,y,a,b,c,d,e,f,g),(a(1,1,1,1,6,6,6,6,6,6,6),dp(4),dp(7)) |
|---|
| 16 | ; |
|---|
| 17 | ideal i= -6s5x+a, -5s4tx-1s5y+b, -4s3t2x-2s4ty+c, |
|---|
| 18 | -3s2t3x-3s3t2y+d, -2st4x-4s2t3y+e, -1t5x-5st4y+f, -6t5y+g; |
|---|
| 19 | ;ideal j=std(i); |
|---|
| 20 | LIB "elim.lib"; |
|---|
| 21 | j=nselect(j,1,4); |
|---|
| 22 | j; |
|---|
| 23 | $ |
|---|
| 24 | ============== |
|---|
| 25 | Macaulay input: |
|---|
| 26 | ============== |
|---|
| 27 | ring R |
|---|
| 28 | |
|---|
| 29 | 11 |
|---|
| 30 | stxya-g |
|---|
| 31 | 1:4 6 |
|---|
| 32 | 4 6 |
|---|
| 33 | ideal i |
|---|
| 34 | 7 |
|---|
| 35 | - -6s5x+a |
|---|
| 36 | - -5s4tx-1s5y+b |
|---|
| 37 | - -4s3t2x-2s4ty+c |
|---|
| 38 | - -3s2t3x-3s3t2y+d |
|---|
| 39 | - -2st4x-4s2t3y+e |
|---|
| 40 | - -1t5x-5st4y+f |
|---|
| 41 | - -6t5y+g |
|---|
| 42 | std i j |
|---|
| 43 | - -> 6[126k]7.8.9.10.11.12.13.14.15.16.17.18.19.20.21.22.23.24.25.26.27. |
|---|
| 44 | - -> 28.29.30.31. |
|---|
| 45 | - -> computation complete after degree 31 |
|---|
| 46 | - -> elapsed time : 2 seconds |
|---|
| 47 | |
|---|
| 48 | |
|---|
| 49 | ------- End of Forwarded Message |
|---|
| 50 | |
|---|
| 51 | |
|---|
