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 | |
---|