source: git/Tst/New/stdZtests.tst

spielwiese
Last change on this file was d4a6939, checked in by Hans Schoenemann <hannes@…>, 5 months ago
prepare for Singular 4.4.0
  • Property mode set to 100644
File size: 20.4 KB
Line 
1LIB "tst.lib";
2tst_init();
3
4LIB "primdec.lib";
5
6//Singular Trac #636
7
8ring 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);
9ideal inI=
103,
11gh+2fi+ej,
12dh+2ci+bj,
13df+2cg+aj,
14de+2bg+ai,
15ce+2bf+ah;
16lead(std(inI));
17kill r;
18
19//Github Adi #39
20
21ring rng = (integer),(x,y,z),(ls(3),C);
22ideal J = -2xz3+x2yz,3z-2xy-3xy2;
23ideal gJ = std(J);
24"Zweite:";
25ideal ggJ = std(gJ);
26lead(gJ);lead(ggJ);
27reduce(J,gJ);
28kill rng;
29
30//Github Adi #32
31
32ring rng = (integer),(x,y,z),dp;
33ideal I = -25*x^3+104*x*y^3-23*y^2*z-128;
34ideal 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;
35ideal gI =  std(I);
36ideal gJ =  std(J);
37reduce (I,gI );
38reduce (J,gJ );
39kill rng;
40
41
42option(redSB);
43
44
45ring rng = (integer),(x,y,z),(dp(3),C);
46ideal I = -27*x+25*z,-46*x*z+35*x,21*x-42;
47ideal gI =  std(I);
48ideal ggI =  std(gI);
49reduce (I,gI );
50kill rng;
51
52//Github Adi #33
53
54ring rng = (integer),(x,y,z),dp;
55option(redSB);
56option(redTail);
57ideal I = -5*z^2-y*z^2, x^3*y*z-z^2+1, y^2*z^2+5*x^2*y*z;
58ideal gI =  std(I);
59ideal ggI =  std(gI);
60lead(gI);
61lead(ggI);
62reduce(I,gI);
63kill rng;
64
65
66ring rng = (integer),(x,y,z),(dp(3),C);
67option(redSB);
68option(redTail);
69ideal 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;
70ideal gI =  std(I);
71ideal ggI =  std(gI);
72lead(gI);
73lead(ggI);
74reduce(I,gI);
75kill rng;
76
77option(redSB);
78
79ring rng = (integer),(x,y,z),(dp(3),C);
80ideal I = -49*y^2+40*z^2+17,-34*z^2-3,-25*y*z-28*z;
81ideal gI =  std(I);
82ideal ggI =  std(gI);
83lead(gI);
84lead(ggI);
85reduce(I,gI);
86kill rng;
87
88//Github Adi #36
89
90ring rng = (0),(x,y,z),(ls(3),C);
91ideal I =
92(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;
93ideal J = (3/53)*x*z+(47/44)*x*y^2*z,(-40/107)*x*z^2+22*x*z^3;
94ideal gI = std(I);
95ideal ggI = std(gI);
96ideal gJ = std(J);
97ideal ggJ = std(gJ);
98lead(gI);lead(ggI);lead(gJ);lead(ggJ);
99reduce (I,gI );
100reduce (J,gJ );
101kill rng;
102
103//Github Adi #37
104
105ring rng = (integer),(x,y,z),dp;
106ideal i = -5*z^2-y*z^2, x^3*y*z-z^2+1, y^2*z^2+5*x^2*y*z;
107ideal gI =  std(i);
108ideal ggI =  std(gI);
109lead(gI);lead(ggI);
110reduce(i,gI);
111kill rng;
112
113//Github Adi #39
114//Problem mit ls aber kein problem mit ds,Ds,lp,dp
115//ls solved: in posinL11Ring one should take the Fdeg instead of pLmCmp
116
117option(redSB);
118
119
120ring rng = (integer),(x,y,z),(ls(3),C);
121ideal I = 1+7*y*z^3,-4*x+2*x*z+14*x^2,-9-10*z^2+8*x*z;
122ideal J = -14*x*z^3+5*x^2*y*z,15*z-7*x*y+12*x*y*z^2-11*x*y^2;
123ideal gI =  std(I);
124ideal ggI =  std(gI);
125ideal gJ =  std(J);
126ideal ggJ =  std(gJ);
127reduce (I,gI );
128reduce (J,gJ );
129idealsEqual(gI,ggI );
130idealsEqual(gJ,ggJ );
131kill rng;
132
133//Github Adi #39
134//redtail hat Problemen wenn man im posinL11Ring pLmCmp nimmt statt FDeg
135
136
137ring rng = (integer),(x,y,z),(ls(3),C);
138option("redTail");
139option(redSB);
140ideal 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;
141ideal gI =  std(I);
142ideal ggI =  std(gI);
143reduce(I,gI);
144idealsEqual(gI,ggI );
145lead(gI);
146kill rng;
147
148
149//Github Adi #40
150
151ring rng = (integer),(x,y,z),(ls(3),C);
152ideal 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;
153ideal J = -1+7*y^3*z,-15+14*x+14*x*y*z,15*z+7*x-10*x*z+5*x*y^2;
154ideal gI =  std(I);
155lead(gI);
156reduce(I,gI);
157ideal gJ =  std(J);
158lead(gJ);
159reduce(J,gJ);
160kill rng;
161
162ring rng = (integer),(x,y,z),(ls(3),C);
163ideal I = 7*x*z+3*x*z^2-x^2,4*z+6*z^2+2*y;
164ideal J = 5-5*z^2+3*y^2*z;
165ideal gI =  std(I);
166lead(gI);
167reduce(I,gI);
168ideal gJ =  std(J);
169lead(gJ);
170reduce(J,gJ);
171kill rng;
172
173//Github Adi #27
174
175ring rng = (integer),(y(1..4),x(1..3)),dp;
176ideal 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);
177ideal J  = std(I);
178lead(J);
179reduce(I,J);
180kill rng;
181
182//Github Adi #20
183
184ring rng = (integer),(xw,xq,xk),(lp(2),lp(1),C);
185ideal I = -4*xk+10;
186ideal J = -7*xw+4*xk,-13;
187ideal IJ = intersect(I,J);
188ideal stdIstdJ = intersect(std(I),std(J) );
189reduce(IJ,stdIstdJ);
190reduce(stdIstdJ,IJ);
191kill rng;
192
193//Github Adi #42
194
195ring rng = (integer),(x,y,z),(dp(1),dp(2),C);
196ideal I = 20*x*z+7*y^2;
197ideal J = 6*x*y+20*z+6,-5*x-16*y-11,8*x*y+2*y^2;
198ideal gI =  std(I);
199lead(gI);
200ideal ggI =  std(gI);
201ideal gJ =  std(J);
202lead(gJ);
203ideal ggJ =  std(gJ);
204std(reduce (I,gI )  );
205std(reduce (J,gJ )  );
206size(gJ)== size(ggJ );
207kill rng;
208
209//Github Adi #38
210
211ring rng =  (integer),(x(1),x(3),x(4),x(5)),(dp(4),C);
212ideal Q = x(4)^2-x(3)*x(5),
213x(3)*x(4)-x(5),
214x(3)^2-x(4),
215x(1)^2-121*x(5)^2-11*x(5)-1;
216Q = std(Q);
217qring qr = Q;
218ideal a = x(5)-2;
219ideal j = 6,x(5)-2,x(4)+2,x(3)-2,x(1)+3;
220ideal aJJ = quotient(a*j,j); // ?
221aJJ;
222std(aJJ);
223kill rng;
224
225//segfaults
226//ok mit neue posinL11Ring
227
228ring rng = (integer),(x(1),x(2),x(3)),(ls(3),C);
229ideal I = -3*x(1)-10*x(1)*x(3);
230ideal 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);
231ideal gI =  std(I);
232ideal ggI =  std(gI);
233ideal gJ =  std(J);
234ideal ggJ =  std(gJ);
235lead(gI);
236lead(gJ);
237std(reduce (I,gI )  );
238std(reduce (J,gJ )  );
239kill rng;
240
241//ok now
242
243ring rng = (integer),(x(1),x(2)),(ls(2),C);
244ideal I = -41*x(1)-10*x(1)^3*x(2),-22*x(2)^2-12*x(1)^3;
245ideal J = 14,35*x(1)*x(2)+12*x(1)*x(2)^3;
246ideal gI =  std(I);
247ideal ggI =  std(gI);
248ideal gJ =  std(J);
249lead(gI);
250lead(gJ);
251std(reduce (I,gI )  );
252std(reduce (J,gJ )  );
253kill rng;
254
255//ok now
256
257ring rng = (integer),(x(1),x(2)),(ls(2),C);
258ideal I = -22-24*x(1)^3,4*x(1)*x(2)^2-19*x(1)^2*x(2);
259ideal J = 16*x(1)*x(2),12-31*x(1)-43*x(1)*x(2);
260ideal gI =  std(I);
261ideal ggI =  std(gI);
262ideal gJ =  std(J);
263lead(gI);
264lead(gJ);
265std(reduce (I,gI )  );
266std(reduce (J,gJ )  );
267kill rng;
268
269
270//in the new posinL11Ring works smooth, in the other one takes forever
271
272ring rng = (integer),(x(1),x(2)),(ls(2),C);
273ideal 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;
274ideal J = -2*x(1)-11*x(1)^2+10*x(1)^2*x(2),-13*x(1)*x(2)-10*x(1)^3;
275ideal gI =  std(I);
276ideal ggI =  std(gI);
277ideal gJ =  std(J);
278lead(gI);
279lead(gJ);
280std(reduce (I,gI )  );
281std(reduce (J,gJ )  );
282kill rng;
283
284//ok now
285
286ring rng = (integer),(x(1),x(2)),(ls(2),C);
287ideal I = -10,15*x(1)*x(2)+8*x(1)*x(2)^2;
288ideal J = 2*x(1)^2*x(2),-3+2*x(1)*x(2);
289ideal gI =  std(I);
290ideal gJ =  std(J);
291lead(gI);
292lead(gJ);
293std(reduce (I,gI )  );
294std(reduce (J,gJ )  );
295kill rng;
296
297//ok now
298
299ring rng = (integer),(x(1)),(ls(1),C);
300ideal I = -12+2*x(1)^2;
301ideal J = -6*x(1)^2+11*x(1)^4,-15;
302ideal gI =  std(I);
303ideal gJ =  std(J);
304lead(gI);
305lead(gJ);
306std(reduce (I,gI )  );
307std(reduce (J,gJ )  );
308kill rng;
309
310//Github Adi #47
311
312ring R = integer,x,dp;
313reduce(x,std(2*x));
314
315ring S = integer,x,ls;
316reduce(x,std(2*x));
317kill R,S;
318
319//Github Adi #46
320
321ring rng = integer,x,ds;
322ideal J = -7, -3-15*x^2, -6+18*x-x^2;
323ideal gJ =  std(J);
324reduce(J,gJ);
325kill rng;
326
327//Github Adi #49
328
329ring rng = (integer),(x,y,z),(dp(3),C);
330ideal Q = x,yz2+5;
331qring QR = std(Q);
332poly a = z;
333std(a);
334ideal J = 5,z,x;
335ideal aJJ =  quotient( a*J, J );
336aJJ;
337kill rng;
338kill QR;
339
340//Github Adi #50
341
342ring rng = (integer),(x(1)),(lp(1),C);
343ideal I = 4*x(1)^3+35*x(1)^2+24*x(1),45*x(1),35;
344ideal J = 18*x(1)^3+43*x(1)^2,-7*x(1)-10,-12*x(1)^2-20*x(1)+13;
345ideal I_J =  intersect(I,J);
346size(I_J);
347kill rng;
348
349//Github Adi #21
350
351ring r = integer,(a,d,g,h,i,j),ws(-1,-1,-1,-1,-1,-1);
352ideal inI= gh+2i, dh+j, d+2g+aj, d+2g+ai;
353ideal g= std(inI);
354reduce(inI,g);
355///
356ring rng = (integer),(x,y,z),(ws(-1,-1,-1),C);
357ideal I  = 6yz-6,6xy-6yz-8y;
358ideal gI  =std(I);
359reduce(I,gI);
360kill rng;
361///
362ring rng = (integer),(x,y,z),(ws(-1,-1,-1),C);
363ideal I  = -10x+15y,4y-7;
364ideal gI  =std(I);
365reduce(I,gI);
366kill rng;
367///
368ring rng = (integer),(x,y,z),(ws(-1,-1,-1),C);
369ideal I  = 9x-14y,15z-2;
370ideal gI  =std(I);
371reduce(I,gI);
372kill rng;
373///
374ring rng = (integer),(x,y),(ws(-1,-1),C);
375ideal I  = 8x+7,-6x-13y;
376ideal gI  =std(I);
377reduce(I,gI);
378kill rng;
379kill r;
380
381// Singular bug #648
382
383ring r=integer,(t,x(1..4)),ws(1,-1,-11,-3,-19);
384ideal inJ =
385 3,
386 t*x(3)*x(4)+2*t*x(1)*x(2),
387 t*x(1)*x(2)^2+2*t^3*x(1)*x(2)*x(3),
388 x(1)^2*x(2)*x(3)*x(4)+2*t^2*x(1)^3*x(2)*x(3);
389ideal J =
390 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),
391 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),
392 t*x(3)*x(4)-3*t^2*x(3)*x(4)-t*x(1)*x(2)-2*x(1)^2,
393 3-t;
394reduce(inJ,J);
395kill r;
396
397//Github Adi #48
398
399ring rng = integer,x,ds;
400ideal I = -8*x+2*x^2, -16+9*x^2-x^3;
401ideal gI =  std(I);
402lead(gI);
403lead(std(gI));
404reduce(I,gI);
405kill rng;
406
407//Github Adi #43
408
409ring rng = integer,(x),lp;
410short = 0 ;
411ideal J = -10*x^2-9*x,-8*x^2-4;
412       ideal gJ =  std(J);
413       ideal ggJ =  std(gJ);
414reduce (J,gJ );
415kill rng;
416
417ring rng = integer,x,lp;
418option(noredSB) ;
419option(noredTail) ;
420ideal I   = -8*x^3+7*x^2-7*x, 8*x^2+10*x-2;
421ideal gI  =  std(I);
422ideal ggI =  std(gI);
423reduce (I,gI ); // <>0 !!
424lead(gI);
425lead(ggI);
426kill rng;
427
428ring rng = integer,x,lp;
429short = 0 ;
430ideal I = 8*x^3-4,-10*x^3+9*x;
431       ideal gI =  std(I);
432       ideal ggI =  std(gI);
433reduce (I,gI );
434size(gI)== size(ggI) ;
435lead(gI);
436lead(ggI);
437kill rng;
438
439//Github Adi #25
440
441ring r = integer,(a,d,g,h,i,j),ws(-1,-1,-1,-1,-1,-1);
442ideal inI= gh+2i, dh+j, d+2g+aj, d+2g+ai;
443ideal g= std(inI);
444lead(g);
445reduce(inI,g);
446kill r;
447
448ring rng = (integer),(x,y,z),(ws(-1,-1,-1),C);
449ideal I = 6yz-6,6xy-6yz-8y;
450ideal gI =std(I);
451lead(I);lead(gI);
452reduce(I,gI);
453kill rng;
454
455//  Yue's Bug
456
457ring r = 0,(t,x(1..4)),ws(1,-1,-1,-1,-1);
458ideal I =
4592*x(1)^2+3*x(1)*x(2)+24*x(3)*x(4),
4608*x(1)^3+x(2)*x(3)*x(4)+18*x(3)^2*x(4);
461I = std(I);
462ring s = integer,(t,x(1..4)),ws(1,-1,-1,-1,-1);
463ideal I = fetch(r,I);
464ideal J = I,2-t;
465ideal J2 = std(J);
466lead(J2);
467reduce(J,J2);
468lead(std(J2));
469kill r;
470
471//  Infinite reduction small example
472
473ring r = integer,(t,x,y),(dp(1),ds);
474ideal i = 6+y+x2,4+x;
475ideal gi = std(i);
476lead(gi);
477reduce(i,gi);
478kill r;
479
480//  Intersect with the elimination takes for ever. With the module is immediate
481
482ring rng = integer,(x),dp;
483ideal I = -3*x^2-8;
484ideal J = 9,5*x^2+8;
485option(intersectSyz);
486ideal I_J =  intersect(I,J);
487lead(I_J);
488kill rng;
489
490//  Github Adi #61
491
492ring rng = integer,x,lp;
493ideal  QQQ = 42*x^2;
494qring QR = QQQ;
495ideal I = 41,-36*x1;
496module syzI = syz (I);
497print(syzI);
498kill rng,QR;
499
500//  Github Adi #60
501
502ring rng = integer,x,dp;
503ideal  QQQ = 4*x;
504qring QR = QQQ;
505ideal I = 13;
506ideal J = -7;
507ideal gJ = std(J);
508ideal I_J =  intersect(I,J);
509ideal IGJ = I+gJ;
510intersect(I_J, IGJ) ;
511kill rng,QR;
512
513//  Github Adi #59
514
515ring rng = integer,(x,y),lp;
516ideal I = -2*y^2+11,14*x^2+5*y;
517ideal J = x^2-3*x,9*x*y+2*y^2;
518ideal I_J =  intersect(I,J);
519lead(I_J);
520reduce(I_J,std(I));
521reduce(I_J,std(J));
522kill rng;
523
524//  Github Adi #58
525
526ring rng = (integer),(x(1)),(lp(1),C);
527short = 0 ;
528ideal  QQQ = 4*x(1);
529qring QR = QQQ;
530ideal I = 13;
531ideal J = -7;
532ideal gI = std(I);
533ideal gJ = std(J);
534ideal I_J =  intersect(I,J);
535intersect(I_J,  I+gJ);
536kill rng;
537
538//  Github Adi #56
539
540ring rng = (integer),(x,y),(ls(2),C);
541ideal I = -8-y^2+13*x^2*y,-12+2*x-9*x^2*y;
542ideal gI =  std(I);
543lead(gI);
544reduce(I,gI);
545kill rng;
546
547ring rng = (integer),(x,y,z),(dp(3),C);
548ideal I = -10*y^2*z^2-7*x*y-2*x+2, 7*x*y^3+11, 15*x*y^2*z-8*y*z-11;
549ideal gI =  std(I);
550reduce(I,gI);
551I = 8*y^2*z+3*x*z, -3*x^2*z-9*y^2*z+3, 5*x*z^2-9*y*z+2*z;
552gI =  std(I);
553reduce(I,gI);
554kill rng;
555
556// Github Adi #51
557
558ring rng = (integer),(x(1),x(2)),(ls(1),ls(1),C);
559ideal I = -36+26*x(1)-29*x(1)^3,26+6*x(2)-13*x(1);
560ideal J = 44+28*x(1)*x(2);
561ideal gI =  std(I);
562ideal gJ =  std(J);
563lead(gI);
564lead(gJ);
565std(gI);
566reduce(I,gI);
567reduce(J,gJ);
568kill rng;
569
570ring rng = (integer),(x,y,z),(Ds(1),Ds(2),C);
571option(redSB);
572ideal J = -3+4*z-x,-2*z+y^2-3*z^3;
573ideal gJ =  std(J);
574lead(gJ);
575reduce (J,gJ ) ;
576kill rng;
577
578ring rng = (integer),(x(1),x(2),x(3)),(is(2),is(1),C);
579ideal I = 2*x(3)+15*x(3)^3+x(1),-15+13*x(3)+6*x(1);
580ideal gI =  std(I);
581lead(gI);
582reduce (I,gI );
583kill rng;
584
585//  Singular Trac: #719
586
587ring r = integer,(t,x(1..5)),ws(1,-1,-1,-1,-1,-1);
588poly h = 3*t*x(2)*x(4)-21*t*x(1)*x(5)-12*t^2*x(1)*x(2)-2*t^2*x(1)*x(3)-6*t^2*x(1)*x(4);
589ideal G =
590x(2)+6*x(4)+7*x(5)+4*t*x(2),
591x(2)+8*x(3)+7*x(5)+t*x(2),
592t*x(1)+3*t*x(2)+t^2*x(1),
593t*x(4)+21*t*x(5)+12*t^2*x(2)+2*t^2*x(3)+7*t^2*x(4),
5942-t;
595ideal Gstd = std(G);
596lead(Gstd);
597reduce(h,Gstd);
598reduce(G,Gstd);
599kill r;
600
601// Memory leak
602
603ring r=integer,t,ds;
604ideal I = 2,3-t;
605int i;
606bigint m=memory(0);
607for(i=1; i<100; i++)
608{
609  ideal II = std(I);
610  II = reduce(I,II);
611  kill II;
612}
613m-memory(0); //  should be 0
614
615// Singular Trac: #747
616
617ring rng = (integer),(x,y),(Ds(2),C);
618intvec op = 67108864,16787537;
619ideal I = 6-x^2*y, -2-3*x+2*y;
620ideal gI =  std(I);
621lead(gI);
622lead(std(gI));
623reduce(I,gI);
624kill rng;
625
626//  Github Adi: #29
627
628ring rng = (integer),(x,y,z),lp;
629option(redSB);
630option(redTail);
631ideal I = -4*z+10;
632ideal J = -7*x+4*z,-13;
633ideal gI = std(I);
634ideal gJ = std(J);
635ideal IJ   = intersect(I,J) ;
636ideal gIJ = intersect( gI,gJ ); // falsch!
637ASSUME(0, idealsEqual(IJ , gIJ )   ) ;
638lead(gI);
639lead(gJ);
640lead(IJ);
641lead(gIJ);
642attrib(gI,"isSB",0);
643attrib(gJ,"isSB",0);
644ideal ngIJ = intersect(gI,gJ ); // wrong
645kill rng;
646
647//  Janko's Beispiel
648
649ring r =integer,(z,x,y),lp;
650poly f = z7x5 + z2xy9 + zy11 + x3y9;
651ideal ii = diff(f,x),diff(f,y),diff(f,z);
652ideal I = std(ii);
653lead(I);
654kill r;
655
656tst_status(1);$;
657
Note: See TracBrowser for help on using the repository browser.