1 | //------------------------------------------------- |
---|
2 | //an error of 0.8.5b: donati@hera.unice.fr |
---|
3 | ring r = 0, (x,y), ds; |
---|
4 | poly f = x4+y5+x2y3; |
---|
5 | ideal i = jacob(f); |
---|
6 | ideal j = std(i); |
---|
7 | j; |
---|
8 | //j[1]=2x3+1xy3 |
---|
9 | //j[2]=3x2y2+5y4 |
---|
10 | //j[3]=10xy4 |
---|
11 | //j[4]=1y6 |
---|
12 | |
---|
13 | reduce(x4+y5+x2y3, j); |
---|
14 | //----------------------------------------------- |
---|
15 | //an error of 0.8.6: (ralf) div by zero -> abort |
---|
16 | random(5,5); |
---|
17 | //------------------------------------------------ |
---|
18 | //an error of 0.8.6: (gerik) result was false: |
---|
19 | f=0; |
---|
20 | f==0; |
---|
21 | //------------------------------------------------ |
---|
22 | //an error of 0.8.7: (hgrass) core dump |
---|
23 | matrix a[2][2]=1,2,3,4; |
---|
24 | matrix b[1][0]; |
---|
25 | b; // could not display strange matrices |
---|
26 | a=b; //could not assign strange matrices |
---|
27 | listvar(all); |
---|
28 | //----------------------------------------------- |
---|
29 | //an error of 0.8.7d (hannes) ordering a incorrect |
---|
30 | ring hs=0,(x,y),(c,a(5,3),dp); |
---|
31 | hs; |
---|
32 | poly f=x2+y3; |
---|
33 | f; |
---|
34 | ord(f); |
---|
35 | deg(f); |
---|
36 | //---------------------------------------------- |
---|
37 | //an error of 0.8.7d (gmg) jump to undefined address |
---|
38 | ring r=0,(x,y),a(2,3); |
---|
39 | poly f=x,y; |
---|
40 | f; |
---|
41 | //--------------------------------------------- |
---|
42 | //an error of 0.8.8c (thomas) (last entry missing) |
---|
43 | ring r; |
---|
44 | ideal i1=x,0,y; |
---|
45 | ideal i2; |
---|
46 | i1+i2; |
---|
47 | kill r; |
---|
48 | //--------------------------------------------- |
---|
49 | //an error of 0.8.8c (greuel) result was 0y |
---|
50 | ring r0=0,(x,y),dp; |
---|
51 | poly f=100x+101y; |
---|
52 | subst(f,x,0); |
---|
53 | //--------------------------------------------- |
---|
54 | //an error of 0.8.8d (greuel) result was totally wrong |
---|
55 | // for mixed orderings |
---|
56 | ring r3 = 32003,(t,x,y),(lp(1),ds); |
---|
57 | ideal i = x-t3-t5, y-t10; |
---|
58 | ideal is = std(i); |
---|
59 | is; |
---|
60 | //--------------------------------------------- |
---|
61 | //an error of 0.8.8.f (greuel) segmentation fault |
---|
62 | ring r4 = 31,(t,x,y,z),dp; |
---|
63 | poly tau=(1/2)*(1+6); |
---|
64 | poly alpha=(1/4)*(2*tau+1); |
---|
65 | ideal pencil = (tau^6-1)*x2y2z2-tau^4*(x4y2+y4z2+z4x2)+tau^2*(x2y4+y2z4+z2x4), |
---|
66 | -t2*(x2+y2+z2-t2)*2; |
---|
67 | matrix ppar[2][1] = 1, alpha; |
---|
68 | ideal f = minor(ppar*matrix(pencil),1); |
---|
69 | module jf = module(jacob(f)); |
---|
70 | module kf =std(jf); |
---|
71 | dim(kf); |
---|
72 | mult(kf); |
---|
73 | degree(kf); |
---|
74 | kill r4; |
---|
75 | //--------------------------------------------- |
---|
76 | //an error of 0.9.2d (gmg) input not a resolvente |
---|
77 | LIB "matrix.lib";LIB "homolog.lib"; |
---|
78 | |
---|
79 | |
---|
80 | int n=5; |
---|
81 | int m=3; |
---|
82 | int N=n*(n+1)/2; //Anz. der Variablen |
---|
83 | ring R=32003,x(1..N),dp; |
---|
84 | matrix X=symmat(n); //aus matrix.lib |
---|
85 | ideal J=minor(X,m+1); |
---|
86 | J=std(J); |
---|
87 | module D=J*freemodule(N)+transpose(jacob(J)); // der Differentialmodul |
---|
88 | D = transpose(D); |
---|
89 | ncols(D); nrows(D); |
---|
90 | option(prot); |
---|
91 | |
---|
92 | module E=Ext(2,D,J); |
---|
93 | |
---|
94 | LIB "tst.lib";tst_status(1);$ |
---|