source: git/Tst/New/stdZtests.tst @ 83b334

spielwiese
Last change on this file since 83b334 was 83b334, checked in by Adi Popescu <adi_popescum@…>, 9 years ago
add tests (don't run them yet till i am ready with pshing all the changes into spielwiese)
  • Property mode set to 100644
File size: 18.1 KB
Line 
1LIB "tst.lib";
2tst_init();
3
4LIB("primdec.lib");
5
6
7//Singular Trac #636
8
9ring r=integer,(t,a,b,c,d,e,f,g,h,i,j),ws(1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1);
10ideal inI=
113,
12gh+2fi+ej,
13dh+2ci+bj,
14df+2cg+aj,
15de+2bg+ai,
16ce+2bf+ah;
17std(inI);
18
19//Github Adi #39
20 
21
22ring rng = (integer),(x,y,z),(ls(3),C);
23int trial = 3;
24//option(redSB);
25system("random", 820791554);
26ideal I = 1+7yz3,-4x+2xz+14x2,-9-10z2+8xz;
27ideal J = -2xz3+x2yz,3z-2xy-3xy2;
28// ideal gI = std(I);
29// ideal ggI = std(gI);
30ideal gJ = std(J);
31"Zweite:";
32ideal ggJ = std(gJ);
33//ASSUME(0, 0== size( std(reduce (I,gI ) )) );
34//ASSUME(0, 0== size( std(reduce (J,gJ ) )) );
35//idealsEqual(gI,ggI );
36gJ;ggJ;
37
38
39
40//Github Adi #32
41
42option(redSB);
43
44ring rng = (integer),(x,y,z),dp;
45ideal I = -25*x^3+104*x*y^3-23*y^2*z-128;
46ideal J = -118*x*y^3*z+50*x*y^2*z^2-50*x*y*z-49*y,149*x*y^2*z^2-142*y^2-112;
47ideal gI =  std(I);
48ideal gJ =  std(J);
49reduce (I,gI );
50reduce (J,gJ );
51
52option(redSB);
53
54
55ring rng = (integer),(x,y,z),(dp(3),C);
56int trial = 0;
57system("random", 912758297);
58ideal I = -27*x+25*z,-46*x*z+35*x,21*x-42;
59       ideal gI =  std(I);
60       ideal ggI =  std(gI);
61
62reduce (I,gI );
63
64
65//Github Adi #33
66
67
68ring rng = (integer),(x,y,z),dp;
69option(redSB);
70option(redTail);
71ideal I = -5*z^2-y*z^2, x^3*y*z-z^2+1, y^2*z^2+5*x^2*y*z;
72ideal gI =  std(I);
73ideal ggI =  std(gI);
74gI;
75ggI;
76
77
78
79ring rng = (integer),(x,y,z),(dp(3),C);
80option(redSB);
81option(redTail);
82ideal I = -5*x*z^2-73*y*z^2, 2*x^3*y*z-129*x*z^2+8, 11*y^3*z^2+85*x^2*y*z;
83ideal gI =  std(I);
84ideal ggI =  std(gI);
85size(gI)==size(ggI);
86
87
88ring rng = (integer),(x,y,z),(dp(3),C);
89option(redSB);
90option(redTail);
91ideal I = -5*x*z^2-73*y*z^2, 68*x^3*y*z-129*x*z^2+98, 11*y^3*z^2+85*x^2*y*z;
92ideal gI =  std(I);
93ideal ggI =  std(gI);
94size(gI)==size(ggI);
95
96
97option(redSB);
98
99
100ring rng = (integer),(x,y,z),(dp(3),C);
101int trial = 0;
102system("random", 784341798);
103ideal I = -49*y^2+40*z^2+17,-34*z^2-3,-25*y*z-28*z;
104       ideal gI =  std(I);
105       ideal ggI =  std(gI);
106gI;
107ggI;
108
109//Github Adi #36
110
111ring rng = (0),(x,y,z),(ls(3),C);
112
113system("random", 157161605);
114ideal I =
115(68/23)*z-(142/133)*y*z-(64/37)*x*y+(65/71)*x*y*z^2+(1/3)*x^2*y*z,(-2/7)*x*y*z-(61/91)*x*y^2*z;
116ideal J = (3/53)*x*z+(47/44)*x*y^2*z,(-40/107)*x*z^2+22*x*z^3;
117ideal gI = std(I);
118ideal ggI = std(gI);
119ideal gJ = std(J);
120ideal ggJ = std(gJ);
1210== size( std(reduce (I,gI ) ));
1220== size( std(reduce (J,gJ ) ));
123
124
125//Github Adi #37
126
127ring rng = (integer),(x,y,z),dp;
128ideal i = -5*z^2-y*z^2, x^3*y*z-z^2+1, y^2*z^2+5*x^2*y*z;
129ideal gI =  std(i);
130ideal ggI =  std(gI);
131gI;ggI;
132
133
134
135//Github Adi #39
136//Problem mit ls aber kein problem mit ds,Ds,lp,dp
137//ls solved: in posinL11Ring one should take the Fdeg instead of pLmCmp
138
139option(redSB);
140
141
142ring rng = (integer),(x,y,z),(ls(3),C);
143int trial = 3;
144system("random", 820791554);
145ideal I = 1+7*y*z^3,-4*x+2*x*z+14*x^2,-9-10*z^2+8*x*z;
146ideal J = -14*x*z^3+5*x^2*y*z,15*z-7*x*y+12*x*y*z^2-11*x*y^2;
147       ideal gI =  std(I);
148       ideal ggI =  std(gI);
149       ideal gJ =  std(J);
150       ideal ggJ =  std(gJ);
1510== size( std(reduce (I,gI ) )  );
1520== size( std(reduce (J,gJ ) )  );
153idealsEqual(gI,ggI );
154idealsEqual(gJ,ggJ );
155
156
157
158//Github Adi #39
159//redtail hat Problemen wenn man im posinL11Ring pLmCmp nimmt statt FDeg
160
161
162ring rng = (integer),(x,y,z),(ls(3),C);
163option("redTail");
164option(redSB);
165system("random", 1022116625);
166ideal I = -14*y*z^2+7*x-11*x*y,-15*x*y,-8*y^2*z^2+5*x*y*z,5*y*z^2-15*y^2*z+8*x^2*z^2+3*x^2*y^2;
167ideal gI =  std(I);
168ideal ggI =  std(gI);
169idealsEqual(gI,ggI );
170gI;
171
172
173//Github Adi #40
174option(redSB);
175
176ring rng = (integer),(x,y,z),(ls(3),C);
177int trial = 3;
178system("random", 158291770);
179ideal I = -10*x*y*z,12*z^2-8*y-14*x*z+14*x^2*y*z,-7*y^4+4*x^2+9*x^2*y;
180ideal J = -1+7*y^3*z,-15+14*x+14*x*y*z,15*z+7*x-10*x*z+5*x*y^2;
181       ideal gI =  std(I);
182       
183 ring rng = (integer),(x,y,z),(ls(3),C);
184 intvec op = 67108866,16787537;
185  option() ;
186//options: intStrategy redefine loadLib usage prompt
187  option(set, op) ;
188 int trial = 53;
189 system("random", 1792192072);
190 ideal I = 7*x*z+3*x*z^2-x^2,4*z+6*z^2+2*y;
191 ideal J = 5-5*z^2+3*y^2*z;
192 ideal gI =  std(I);
193       
194
195
196
197//Github Adi #27
198
199ring rng = (integer),(y(1..4),x(1..3)),dp;
200ideal I  = 12*y(3)*y(4)+19*y(1)+6*y(4),2*x(3)*y(4)-2*y(1)+7,2*y(1)*y(3)+y(1)-7*y(3)+6,7*x(3)*y(3)-6*x(3)-6*y(3)-3,y(1)^2+126*y(3)*y(4)+196*y(1)+60*y(4),x(3)*y(1)-3,x(3)*y(3)*y(4)+30*y(3)*y(4)+48*y(1)+15*y(4)+3,48*y(3)*y(4)^2-76*y(1)*y(4)+24*y(4)^2+266*y(4),48*y(3)^2*y(4)+48*y(3)*y(4)+12*y(4),48*y(1)*y(3)*y(4)+24*y(1)*y(4)-168*y(3)*y(4)-84*y(4),48*x(3)*y(3)*y(4)-76*x(3)*y(1)-48*y(1)*y(3)+24*x(3)*y(4)+266*x(3)-24*y(1)+168*y(3)+84,-266*x(3)*y(4)+168*y(3)*y(4)+84*y(4),266*x(3)*y(3)+168*y(3)^2-228*x(3)-60*y(3)-72,336*x(3)*y(3)*y(4)-60*x(3)*y(4)-144*y(3)*y(4)-72*y(4),-266*x(3)^2-168*x(3)*y(3)+144*x(3)+144*y(3)+72,-168*y(3)^2*y(4)+266*y(1)*y(3)+60*y(3)*y(4)-228*y(1)-931*y(3)+72*y(4)+798,-168*y(3)^3-24*y(3)^2+102*y(3)+36,-168*y(1)*y(3)^2+60*y(1)*y(3)+588*y(3)^2+72*y(1)-210*y(3)-252,168*x(3)*y(3)^2+168*x(3)*y(3)-144*y(3)^2+42*x(3)-144*y(3)-36,28*x(3)*y(3)*y(4)-28*y(1)*y(3)+14*x(3)*y(4)-24*y(3)*y(4)+24*y(1)+98*y(3)-12*y(4)-84,-28*x(3)*y(3)^2+10*x(3)*y(3)+12*x(3),-28*x(3)*y(1)*y(3)+24*x(3)*y(1)+98*x(3)*y(3)-84*x(3),28*x(3)^2*y(3)+14*x(3)^2-24*x(3)*y(3)-12*x(3),48*y(1)*y(4)^2-84*y(4)^2,48*y(1)*y(3)*y(4)+24*y(1)*y(4)-84*y(3)*y(4)+72*y(4),72*y(4)^2,48*x(3)*y(1)*y(4)-84*x(3)*y(4)-72*y(4),-48*y(1)^2*y(3)-24*y(1)^2+252*y(1)*y(3)+72*y(3)*y(4)+12*y(1)-294*y(3)+36*y(4)+252,-8*y(1)^2*y(4)+504*y(3)*y(4)^2-756*y(1)*y(4)+252*y(4)^2+2744*y(4),504*y(3)^2*y(4)+492*y(3)*y(4)+120*y(4),504*y(1)*y(3)*y(4)+240*y(1)*y(4)-1764*y(3)*y(4)-840*y(4),-8*x(3)*y(1)^2+504*x(3)*y(3)*y(4)-756*x(3)*y(1)-504*y(1)*y(3)+252*x(3)*y(4)+2744*x(3)-240*y(1)+1764*y(3)+840,168*y(1)*y(3)*y(4)-144*y(1)*y(4)-294*y(3)*y(4)+252*y(4),168*y(1)*y(3)^2-60*y(1)*y(3)-294*y(3)^2-72*y(1)+504*y(3)-216,252*y(3)*y(4)-216*y(4),-168*x(3)*y(1)*y(3)-84*x(3)*y(1)+294*x(3)*y(3)+144*y(1)*y(3)-252*x(3)+72*y(1)-252*y(3)+216,-28*x(3)*y(1)*y(4)-2744*x(3)*y(4)+24*y(1)*y(4)+1764*y(3)*y(4)+840*y(4),28*x(3)*y(1)*y(3)-24*x(3)*y(1)+2744*x(3)*y(3)+1764*y(3)^2-2352*x(3)-672*y(3)-720,3528*x(3)*y(3)*y(4)-672*x(3)*y(4)-1512*y(3)*y(4)-720*y(4),-28*x(3)^2*y(1)-2744*x(3)^2+24*x(3)*y(1)-1764*x(3)*y(3)+1512*x(3)+1512*y(3)+720,28*y(1)^2*y(3)-1764*y(3)^2*y(4)-24*y(1)^2+2646*y(1)*y(3)+630*y(3)*y(4)-2268*y(1)-9604*y(3)+756*y(4)+8232,-1764*y(3)^3-210*y(3)^2+1056*y(3)+360,-1764*y(1)*y(3)^2+672*y(1)*y(3)+6174*y(3)^2+720*y(1)-2352*y(3)-2520,1764*x(3)*y(3)^2+1722*x(3)*y(3)-1512*y(3)^2+420*x(3)-1476*y(3)-360,24*x(3)*y(4)^2+24*y(1)*y(4),24*x(3)*y(3)*y(4)+38*x(3)*y(1)+24*y(1)*y(3)+12*x(3)*y(4)+12*y(1),24*x(3)*y(1)*y(4),24*x(3)^2*y(4),-24*y(1)*y(3)*y(4)+38*y(1)^2-12*y(1)*y(4)-133*y(1),-24*y(1)*y(3)^2-24*y(1)*y(3)-6*y(1),-24*y(1)^2*y(3)-12*y(1)^2+84*y(1)*y(3)+42*y(1),-24*x(3)*y(1)*y(3)-12*x(3)*y(1)+84*x(3)*y(3)+42*x(3),-4*x(3)*y(1)*y(4)-4*y(1)^2+14*x(3)*y(4)+14*y(1),-4*x(3)*y(1)*y(3)-2*x(3)*y(1),-4*x(3)*y(1)^2+14*x(3)*y(1),-4*x(3)^2*y(1)+14*x(3)^2,84*x(3)*y(3)*y(4)+133*x(3)*y(1)-72*x(3)*y(4)-114*y(1),84*x(3)*y(3)^2-30*x(3)*y(3)-36*x(3),-84*x(3)*y(1)*y(3)-42*x(3)*y(1)+72*y(1)*y(3)+36*y(1),-84*x(3)^2*y(3)-42*x(3)^2+72*x(3)*y(3)+36*x(3),-14*x(3)^2*y(4)-14*x(3)*y(1)+12*x(3)*y(4)+12*y(1),14*x(3)^2*y(3)-12*x(3)^2,-14*x(3)^2*y(1)+12*x(3)*y(1),-14*x(3)^3+12*x(3)^2,28*x(3)*y(1)*y(3)-5*x(3)*y(1)-49*x(3)*y(3)-12*y(1)*y(3)+42*x(3)-6*y(1),-24*y(1)^2*y(4)+42*y(1)*y(4),-24*y(1)^2*y(3)-12*y(1)^2+42*y(1)*y(3)-36*y(1),-36*y(1)*y(4),-36*x(3)*y(4),-252*x(3)*y(4)^2-252*y(1)*y(4),-4*x(3)*y(1)^2-252*x(3)*y(3)*y(4)-392*x(3)*y(1)-252*y(1)*y(3)-120*x(3)*y(4)-120*y(1),-252*x(3)*y(1)*y(4),-252*x(3)^2*y(4),-4*y(1)^3+252*y(1)*y(3)*y(4)-378*y(1)^2+126*y(1)*y(4)+1372*y(1),252*y(1)*y(3)^2+246*y(1)*y(3)+60*y(1),252*y(1)^2*y(3)+120*y(1)^2-882*y(1)*y(3)-420*y(1),252*x(3)*y(1)*y(3)+120*x(3)*y(1)-882*x(3)*y(3)-420*x(3),-14*x(3)*y(1)^2-882*x(3)*y(3)*y(4)-1372*x(3)*y(1)+12*y(1)^2+756*x(3)*y(4)+1176*y(1),-882*x(3)*y(3)^2+336*x(3)*y(3)+360*x(3),882*x(3)*y(1)*y(3)+420*x(3)*y(1)-756*y(1)*y(3)-360*y(1),882*x(3)^2*y(3)+420*x(3)^2-756*x(3)*y(3)-360*x(3),-38*x(3)*y(4)^2+24*y(3)*y(4)^2+12*y(4)^2,24*y(3)^2*y(4)+24*y(3)*y(4)+6*y(4),24*x(3)*y(3)*y(4)^2+12*x(3)*y(4)^2,-38*x(3)^2*y(4),-24*y(3)^2*y(4)^2+38*y(1)*y(3)*y(4)-12*y(3)*y(4)^2-133*y(3)*y(4),-24*y(3)^3*y(4)-24*y(3)^2*y(4)-6*y(3)*y(4),-24*y(1)*y(3)^2*y(4)-12*y(1)*y(3)*y(4)+84*y(3)^2*y(4)+42*y(3)*y(4),38*x(3)*y(1)*y(3)+12*x(3)*y(3)*y(4)-133*x(3)*y(3)-12*y(1)*y(3)+6*x(3)*y(4)-6*y(1)+42*y(3)+21,4*x(3)*y(3)*y(4)^2-4*y(1)*y(3)*y(4)+2*x(3)*y(4)^2+120*y(3)*y(4)^2-192*y(1)*y(4)+14*y(3)*y(4)+60*y(4)^2+672*y(4),120*y(3)^2*y(4)+120*y(3)*y(4)+30*y(4),120*y(1)*y(3)*y(4)+60*y(1)*y(4)-420*y(3)*y(4)-210*y(4),4*x(3)^2*y(3)*y(4)-4*x(3)*y(1)*y(3)+2*x(3)^2*y(4)+120*x(3)*y(3)*y(4)-192*x(3)*y(1)+14*x(3)*y(3)-120*y(1)*y(3)+60*x(3)*y(4)+672*x(3)-60*y(1)+420*y(3)+210,114*x(3)*y(4)-72*y(3)*y(4)-36*y(4),-133*x(3)*y(3)^2+114*x(3)*y(3)+42*y(3)^2-15*y(3)-18,-84*x(3)*y(3)^2*y(4)+72*y(3)^2*y(4)-36*x(3)*y(4)+36*y(3)*y(4),133*x(3)^2*y(3)-156*x(3)*y(3)-21*x(3)+36*y(3)+18,-684*x(3)*y(4)+432*y(3)*y(4)+216*y(4),14*x(3)*y(3)^2+660*x(3)*y(3)+420*y(3)^2-576*x(3)-150*y(3)-180,14*x(3)^2*y(3)*y(4)-12*x(3)^2*y(4)+840*x(3)*y(3)*y(4)-150*x(3)*y(4)-360*y(3)*y(4)-180*y(4),-14*x(3)^2*y(3)-672*x(3)^2-408*x(3)*y(3)+366*x(3)+360*y(3)+180,12*x(3)*y(3)*y(4)-432*y(3)^2*y(4)+672*y(1)*y(3)+6*x(3)*y(4)+144*y(3)*y(4)-576*y(1)-2352*y(3)+180*y(4)+2016,-14*x(3)*y(3)^3+5*x(3)*y(3)^2-420*y(3)^3+6*x(3)*y(3)-60*y(3)^2+255*y(3)+90,-14*x(3)*y(1)*y(3)^2+12*x(3)*y(1)*y(3)+49*x(3)*y(3)^2-420*y(1)*y(3)^2-42*x(3)*y(3)+150*y(1)*y(3)+1470*y(3)^2+180*y(1)-525*y(3)-630,14*x(3)^2*y(3)^2+7*x(3)^2*y(3)+408*x(3)*y(3)^2+414*x(3)*y(3)-360*y(3)^2+105*x(3)-360*y(3)-90,-24*y(1)*y(3)*y(4)^2+42*y(3)*y(4)^2,-24*y(1)*y(3)^2*y(4)-12*y(1)*y(3)*y(4)+42*y(3)^2*y(4)-36*y(3)*y(4),-36*y(3)*y(4)^2,12*x(3)*y(1)*y(4)+36*x(3)*y(4)-18*y(4),4*x(3)*y(1)*y(4)^2+392*x(3)*y(4)^2+120*y(1)*y(4)^2-252*y(3)*y(4)^2-336*y(4)^2,120*y(1)*y(3)*y(4)-252*y(3)^2*y(4)+60*y(1)*y(4)-456*y(3)*y(4)+120*y(4),-252*x(3)*y(3)*y(4)^2-120*x(3)*y(4)^2+180*y(4)^2,4*x(3)^2*y(1)*y(4)+392*x(3)^2*y(4)+120*x(3)*y(1)*y(4)-216*x(3)*y(4)-180*y(4),-4*y(1)^2*y(3)*y(4)+252*y(3)^2*y(4)^2-378*y(1)*y(3)*y(4)+126*y(3)*y(4)^2+1372*y(3)*y(4),252*y(3)^3*y(4)+246*y(3)^2*y(4)+60*y(3)*y(4),252*y(1)*y(3)^2*y(4)+120*y(1)*y(3)*y(4)-882*y(3)^2*y(4)-420*y(3)*y(4),-4*x(3)*y(1)^2*y(3)-378*x(3)*y(1)*y(3)-120*y(1)^2*y(3)-120*x(3)*y(3)*y(4)-60*y(1)^2+1372*x(3)*y(3)+756*y(1)*y(3)-60*x(3)*y(4)+180*y(3)*y(4)+90*y(1)-1176*y(3)+90*y(4)+420,-12*x(3)*y(1)*y(4)+432*y(1)*y(3)*y(4)-1176*x(3)*y(4)-360*y(1)*y(4)+1008*y(4),14*x(3)*y(1)*y(3)^2-12*x(3)*y(1)*y(3)+1372*x(3)*y(3)^2+420*y(1)*y(3)^2-1176*x(3)*y(3)-150*y(1)*y(3)-1176*y(3)^2-180*y(1)+1428*y(3)-360,882*x(3)*y(3)^2*y(4)-756*y(3)^2*y(4)+360*x(3)*y(4)+270*y(3)*y(4)-540*y(4),-14*x(3)^2*y(1)*y(3)-1372*x(3)^2*y(3)-408*x(3)*y(1)*y(3)-210*x(3)*y(1)+2352*x(3)*y(3)+360*y(1)*y(3)-420*x(3)+180*y(1)-1008*y(3)+360,-12*x(3)*y(3)*y(4)^2-19*x(3)*y(1)*y(4)+6*y(1)*y(4),-12*x(3)*y(3)^2*y(4)-19*x(3)*y(1)*y(3)-6*x(3)*y(3)*y(4)+6*y(1)*y(3)+3*y(1),6*x(3)*y(1)*y(4),6*x(3)^2*y(4),2*x(3)^2*y(4)^2+2*x(3)*y(1)*y(4)+60*x(3)*y(4)^2+60*y(1)*y(4),2*x(3)*y(1)*y(3)+60*x(3)*y(3)*y(4)+96*x(3)*y(1)+60*y(1)*y(3)+30*x(3)*y(4)+30*y(1),2*x(3)^2*y(1)*y(4)+60*x(3)*y(1)*y(4),2*x(3)^3*y(4)+60*x(3)^2*y(4),-4*x(3)*y(1)*y(3)*y(4)-x(3)*y(1)*y(4)+7*x(3)*y(3)*y(4)-60*y(1)*y(3)*y(4)+96*y(1)^2-30*y(1)*y(4)-336*y(1),-2*x(3)*y(1)*y(3)^2-x(3)*y(1)*y(3)-60*y(1)*y(3)^2-60*y(1)*y(3)-15*y(1),-2*x(3)*y(1)^2*y(3)+7*x(3)*y(1)*y(3)-60*y(1)^2*y(3)-30*y(1)^2+210*y(1)*y(3)+105*y(1),-2*x(3)^2*y(1)*y(3)+7*x(3)^2*y(3)-60*x(3)*y(1)*y(3)-30*x(3)*y(1)+210*x(3)*y(3)+105*x(3),-6*x(3)^2*y(4)+216*x(3)*y(3)*y(4)+336*x(3)*y(1)-180*x(3)*y(4)-288*y(1),7*x(3)^2*y(3)^2-6*x(3)^2*y(3)+210*x(3)*y(3)^2-75*x(3)*y(3)-90*x(3),-7*x(3)^2*y(1)*y(3)-204*x(3)*y(1)*y(3)-105*x(3)*y(1)+180*y(1)*y(3)+90*y(1),-7*x(3)^3*y(3)-204*x(3)^2*y(3)-105*x(3)^2+180*x(3)*y(3)+90*x(3),-2*x(3)*y(1)^2*y(4)-126*x(3)*y(3)*y(4)^2-196*x(3)*y(1)*y(4)-60*y(1)^2*y(4)+168*y(1)*y(4),-2*x(3)*y(1)^2*y(3)-126*x(3)*y(3)^2*y(4)-196*x(3)*y(1)*y(3)-60*y(1)^2*y(3)-60*x(3)*y(3)*y(4)-30*y(1)^2+168*y(1)*y(3)-60*y(1),60*x(3)*y(1)*y(4)-90*y(1)*y(4),60*x(3)^2*y(4)-90*x(3)*y(4);
201ideal J  = std(I);
202
203
204
205
206//Github Adi #20
207ring rng = (integer),(xw,xq,xk),(lp(2),lp(1),C);
208ideal I = -4*xk+10;
209ideal J = -7*xw+4*xk,-13;
210idealsEqual( intersect(I,J) , intersect(std(I),std(J) ) );
211
212
213
214
215
216
217
218//Github Adi #42
219ring rng = (integer),(x,y,z),(dp(1),dp(2),C);
220intvec op = 100663296,16787537;
221 option() ;
222 option(set, op) ;
223int trial = 0;
224system("random", 1310165262);
225ideal I = 20*x*z+7*y^2;
226ideal J = 6*x*y+20*z+6,-5*x-16*y-11,8*x*y+2*y^2;
227       ideal gI =  std(I);
228       ideal ggI =  std(gI);
229       ideal gJ =  std(J);
230       ideal ggJ =  std(gJ);
231std(reduce (I,gI )  );
232std(reduce (J,gJ )  );
233size(gJ)== size(ggJ );
234   
235
236
237
238
239//Github Adi #38
240  ring rng =  (integer),(x(1),x(3),x(4),x(5)),(dp(4),C);
241ideal Q = x(4)^2-x(3)*x(5),
242x(3)*x(4)-x(5),
243x(3)^2-x(4),
244x(1)^2-121*x(5)^2-11*x(5)-1;
245
246Q = std(Q);
247Q;
248qring qr = Q;
249
250 ideal a = x(5)-2;
251 ideal j = 6,x(5)-2,x(4)+2,x(3)-2,x(1)+3;
252 ideal aJJ = quotient(a*j,j); // ?
253 aJJ;
254 std(aJJ);
255 aJJ;
256
257option(redSB);
258
259//segfaults
260
261
262//# segfaults
263
264//ok mit neue posinL11Ring
265
266ring rng = (integer),(x(1),x(2),x(3)),(ls(3),C);
267short = 0 ;
268minpoly = 0;
269intvec op = 67108864,16787537;
270 option(set, op) ;
271 option() ;
272 // option() ;//options: intStrategy redefine loadLib usage prompt warn
273int trial = 96;
274system("random", 970755877);
275ideal I = -3*x(1)-10*x(1)*x(3);
276ideal J = -9*x(2)+9*x(1)*x(2),-3*x(2)*x(3)^2+10*x(1)*x(3)^2+10*x(1)^2*x(2);
277       ideal gI =  std(I);
278       ideal ggI =  std(gI);
279       ideal gJ =  std(J);
280       ideal ggJ =  std(gJ);
281
282
283//ok now
284ring rng = (integer),(x(1),x(2)),(ls(2),C);
285short = 0 ;
286minpoly = 0;
287intvec op = 67108864,16787537;
288 option(set, op) ;
289 option() ;
290 // option() ;//options: intStrategy redefine loadLib usage prompt warn
291int trial = 4;
292system("random", 960625435);
293ideal I = -41*x(1)-10*x(1)^3*x(2),-22*x(2)^2-12*x(1)^3;
294ideal J = 14,35*x(1)*x(2)+12*x(1)*x(2)^3;
295       ideal gI =  std(I);
296       ideal ggI =  std(gI);
297       ideal gJ =  std(J);
298
299//ok now
300
301ring rng = (integer),(x(1),x(2)),(ls(2),C);
302short = 0 ;
303minpoly = 0;
304intvec op = 67108864,16787537;
305 option(set, op) ;
306 option() ;
307 // option() ;//options: intStrategy redefine loadLib usage prompt warn
308int trial = 29;
309system("random", 550741380);
310ideal I = -22-24*x(1)^3,4*x(1)*x(2)^2-19*x(1)^2*x(2);
311ideal J = 16*x(1)*x(2),12-31*x(1)-43*x(1)*x(2);
312       ideal gI =  std(I);
313       ideal ggI =  std(gI);
314       ideal gJ =  std(J);
315
316//in the new posinL11Ring works smooth, in the opther one takes forever
317
318ring rng = (integer),(x(1),x(2)),(ls(2),C);
319short = 0 ;
320minpoly = 0;
321intvec op = 67108864,16787537;
322 option(set, op) ;
323 option() ;
324 // option() ;//options: intStrategy redefine loadLib usage prompt warn
325int trial = 15;
326system("random", 1976068534);
327ideal I = -11*x(2),3*x(2)^2+8*x(1)*x(2)+3*x(1)*x(2)^2+12*x(1)^3,8*x(1)*x(2)^2-15*x(1)^3;
328ideal J = -2*x(1)-11*x(1)^2+10*x(1)^2*x(2),-13*x(1)*x(2)-10*x(1)^3;
329       ideal gI =  std(I);
330
331//ok now
332
333ring rng = (integer),(x(1),x(2)),(ls(2),C);
334short = 0 ;
335minpoly = 0;
336intvec op = 67108864,16787537;
337 option(set, op) ;
338 option() ;
339 // option() ;//options: intStrategy redefine loadLib usage prompt warn
340int trial = 68;
341system("random", 1702067426);
342ideal I = -10,15*x(1)*x(2)+8*x(1)*x(2)^2;
343ideal J = 2*x(1)^2*x(2),-3+2*x(1)*x(2);
344       ideal gI =  std(I);
345
346
347//ok now
348
349ring rng = (integer),(x(1)),(ls(1),C);
350short = 0 ;
351minpoly = 0;
352intvec op = 67108864,16787537;
353 option(set, op) ;
354 option() ;
355 // option() ;//options: intStrategy redefine loadLib usage prompt warn
356int trial = 65;
357system("random", 907801360);
358ideal I = -12+2*x(1)^2;
359ideal J = -6*x(1)^2+11*x(1)^4,-15;
360       ideal gI =  std(I);
361       ideal ggI =  std(gI);
362       ideal gJ =  std(J);
363
364
365
366//Github Adi #47
367
368ring R = integer,x,dp;
369reduce(x,std(2*x));
370
371ring S = integer,x,ls;
372reduce(x,std(2*x));
373
374
375
376//Github Adi #46
377
378ring rng = integer,x,ds;
379short = 0 ;
380ideal J = -7, -3-15*x^2, -6+18*x-x^2;
381ideal gJ =  std(J);
382
383//Github Adi #49
384
385ring rng = (integer),(x,y,z),(dp(3),C);
386ideal Q = x,yz2+5;
387qring QR = std(Q);
388poly a = z;
389ideal J = 5,z,x;
390ideal aJJ =  quotient( a*J, J );
391aJJ;
392
393//Github Adi #50
394
395ring rng = (integer),(x(1)),(lp(1),C);
396ideal I = 4*x(1)^3+35*x(1)^2+24*x(1),45*x(1),35;
397ideal J = 18*x(1)^3+43*x(1)^2,-7*x(1)-10,-12*x(1)^2-20*x(1)+13;
398ideal I_J =  intersect(I,J);
399size(I_J);
400
401//Github Adi #21
402
403ring r = integer,(a,d,g,h,i,j),ws(-1,-1,-1,-1,-1,-1);
404ideal inI= gh+2i, dh+j, d+2g+aj, d+2g+ai;
405ideal g= std(inI);
406///
407ring rng = (integer),(x,y,z),(ws(-1,-1,-1),C);
408ideal I  = 6yz-6,6xy-6yz-8y;
409ideal gI  =std(I);
410///
411ring rng = (integer),(x,y,z),(ws(-1,-1,-1),C);
412ideal I  = -10x+15y,4y-7;
413ideal gI  =std(I);
414///
415ring rng = (integer),(x,y,z),(ws(-1,-1,-1),C);
416ideal I  = 9x-14y,15z-2;
417ideal gI  =std(I);
418///
419ring rng = (integer),(x,y),(ws(-1,-1),C);
420ideal I  = 8x+7,-6x-13y;
421ideal gI  =std(I);
422
423// Singular bug #648
424
425ring r=integer,(t,x(1..4)),ws(1,-1,-11,-3,-19);
426option(prot);
427ideal inJ =
428 3,
429 t*x(3)*x(4)+2*t*x(1)*x(2),
430 t*x(1)*x(2)^2+2*t^3*x(1)*x(2)*x(3),
431 x(1)^2*x(2)*x(3)*x(4)+2*t^2*x(1)^3*x(2)*x(3);
432ideal J =
433 x(1)^2*x(2)*x(3)*x(4)-t*x(1)^2*x(2)^3+5*t*x(1)^2*x(2)*x(3)*x(4)+16*t^3*x(1)*x(2)*x(3)^2*x(4)-2*x(1)^3*x(2)^2-1536*x(1)^3*x(3)*x(4)+32*t^2*x(1)^2*x(3)^2*x(4),
434 t*x(1)*x(2)^2+2*x(1)^2*x(2)+2*t^3*x(1)*x(2)*x(3)-192*x(1)^3+4*t^2*x(1)^2*x(3),
435 t*x(3)*x(4)-3*t^2*x(3)*x(4)-t*x(1)*x(2)-2*x(1)^2,
436 3-t;
437reduce(inJ,J);
438
439//Github Adi #48
440
441ring rng = integer,x,ds;
442ideal I = -8*x+2*x^2, -16+9*x^2-x^3;
443ideal gI =  std(I);
444gI;
445std(gI);
446
447//Github Adi #43
448
449LIB "primdec.lib";
450ring rng = integer,(x),lp;
451short = 0 ;
452ideal J = -10*x^2-9*x,-8*x^2-4;
453       ideal gJ =  std(J);
454       ideal ggJ =  std(gJ);
455ASSUME(0, 0== size( std(reduce (J,gJ ) ))  );
456ASSUME(0, idealsEqual(gJ,ggJ ));
457gJ;
458ggJ;
459
460ring rng = integer,x,lp;
461short = 0 ;
462option(noredSB) ;
463option(noredTail) ;
464ideal I   = -8*x^3+7*x^2-7*x, 8*x^2+10*x-2;
465ideal gI  =  std(I);
466ideal ggI =  std(gI);
467reduce (I,gI ); // <>0 !!
468gI;
469ggI;
470
471LIB("primdec.lib");
472ring rng = integer,x,lp;
473short = 0 ;
474ideal I = 8*x^3-4,-10*x^3+9*x;
475       ideal gI =  std(I);
476       ideal ggI =  std(gI);
477ASSUME(0, 0== size( std(reduce (I,gI ) ))  );
478ASSUME(0, idealsEqual(gI,ggI ));
479size(gI)== size(ggI) ;
480gI;
481ggI;
482
483//Github Adi #25
484
485ring r = integer,(a,d,g,h,i,j),ws(-1,-1,-1,-1,-1,-1);
486ideal inI= gh+2i, dh+j, d+2g+aj, d+2g+ai;
487ideal g= std(inI);
488inI;
489g;
490
491ring rng = (integer),(x,y,z),(ws(-1,-1,-1),C);
492ideal I = 6yz-6,6xy-6yz-8y;
493ideal gI =std(I);
494I;gI;
495
496
497 
498       
499tst_status(1);$;     
500
Note: See TracBrowser for help on using the repository browser.