[b35b93] | 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 | |
---|
[46976a6] | 94 | LIB "tst.lib";tst_status(1);$ |
---|