1 | LIB "tst.lib"; |
---|
2 | tst_init(); |
---|
3 | |
---|
4 | LIB("primdec.lib"); |
---|
5 | |
---|
6 | |
---|
7 | //Singular Trac #636 |
---|
8 | |
---|
9 | ring 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); |
---|
10 | ideal inI= |
---|
11 | 3, |
---|
12 | gh+2fi+ej, |
---|
13 | dh+2ci+bj, |
---|
14 | df+2cg+aj, |
---|
15 | de+2bg+ai, |
---|
16 | ce+2bf+ah; |
---|
17 | std(inI); |
---|
18 | |
---|
19 | //Github Adi #39 |
---|
20 | |
---|
21 | |
---|
22 | ring rng = (integer),(x,y,z),(ls(3),C); |
---|
23 | ideal J = -2xz3+x2yz,3z-2xy-3xy2; |
---|
24 | ideal gJ = std(J); |
---|
25 | "Zweite:"; |
---|
26 | ideal ggJ = std(gJ); |
---|
27 | gJ;ggJ; |
---|
28 | |
---|
29 | |
---|
30 | |
---|
31 | //Github Adi #32 |
---|
32 | |
---|
33 | ring rng = (integer),(x,y,z),dp; |
---|
34 | ideal I = -25*x^3+104*x*y^3-23*y^2*z-128; |
---|
35 | ideal 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; |
---|
36 | ideal gI = std(I); |
---|
37 | ideal gJ = std(J); |
---|
38 | reduce (I,gI ); |
---|
39 | reduce (J,gJ ); |
---|
40 | |
---|
41 | option(redSB); |
---|
42 | |
---|
43 | |
---|
44 | ring rng = (integer),(x,y,z),(dp(3),C); |
---|
45 | ideal I = -27*x+25*z,-46*x*z+35*x,21*x-42; |
---|
46 | ideal gI = std(I); |
---|
47 | ideal ggI = std(gI); |
---|
48 | reduce (I,gI ); |
---|
49 | |
---|
50 | |
---|
51 | //Github Adi #33 |
---|
52 | |
---|
53 | |
---|
54 | ring rng = (integer),(x,y,z),dp; |
---|
55 | option(redSB); |
---|
56 | option(redTail); |
---|
57 | ideal I = -5*z^2-y*z^2, x^3*y*z-z^2+1, y^2*z^2+5*x^2*y*z; |
---|
58 | ideal gI = std(I); |
---|
59 | ideal ggI = std(gI); |
---|
60 | gI; |
---|
61 | ggI; |
---|
62 | |
---|
63 | |
---|
64 | |
---|
65 | ring rng = (integer),(x,y,z),(dp(3),C); |
---|
66 | option(redSB); |
---|
67 | option(redTail); |
---|
68 | ideal 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; |
---|
69 | ideal gI = std(I); |
---|
70 | ideal ggI = std(gI); |
---|
71 | size(gI)==size(ggI); |
---|
72 | |
---|
73 | |
---|
74 | ring rng = (integer),(x,y,z),(dp(3),C); |
---|
75 | option(redSB); |
---|
76 | option(redTail); |
---|
77 | ideal 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; |
---|
78 | ideal gI = std(I); |
---|
79 | ideal ggI = std(gI); |
---|
80 | size(gI)==size(ggI); |
---|
81 | |
---|
82 | |
---|
83 | option(redSB); |
---|
84 | |
---|
85 | |
---|
86 | ring rng = (integer),(x,y,z),(dp(3),C); |
---|
87 | ideal I = -49*y^2+40*z^2+17,-34*z^2-3,-25*y*z-28*z; |
---|
88 | ideal gI = std(I); |
---|
89 | ideal ggI = std(gI); |
---|
90 | gI; |
---|
91 | ggI; |
---|
92 | |
---|
93 | //Github Adi #36 |
---|
94 | |
---|
95 | ring rng = (0),(x,y,z),(ls(3),C); |
---|
96 | ideal I = |
---|
97 | (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; |
---|
98 | ideal J = (3/53)*x*z+(47/44)*x*y^2*z,(-40/107)*x*z^2+22*x*z^3; |
---|
99 | ideal gI = std(I); |
---|
100 | ideal ggI = std(gI); |
---|
101 | ideal gJ = std(J); |
---|
102 | ideal ggJ = std(gJ); |
---|
103 | 0== size( std(reduce (I,gI ) )); |
---|
104 | 0== size( std(reduce (J,gJ ) )); |
---|
105 | |
---|
106 | |
---|
107 | //Github Adi #37 |
---|
108 | |
---|
109 | ring rng = (integer),(x,y,z),dp; |
---|
110 | ideal i = -5*z^2-y*z^2, x^3*y*z-z^2+1, y^2*z^2+5*x^2*y*z; |
---|
111 | ideal gI = std(i); |
---|
112 | ideal ggI = std(gI); |
---|
113 | gI;ggI; |
---|
114 | |
---|
115 | |
---|
116 | |
---|
117 | //Github Adi #39 |
---|
118 | //Problem mit ls aber kein problem mit ds,Ds,lp,dp |
---|
119 | //ls solved: in posinL11Ring one should take the Fdeg instead of pLmCmp |
---|
120 | |
---|
121 | option(redSB); |
---|
122 | |
---|
123 | |
---|
124 | ring rng = (integer),(x,y,z),(ls(3),C); |
---|
125 | ideal I = 1+7*y*z^3,-4*x+2*x*z+14*x^2,-9-10*z^2+8*x*z; |
---|
126 | ideal J = -14*x*z^3+5*x^2*y*z,15*z-7*x*y+12*x*y*z^2-11*x*y^2; |
---|
127 | ideal gI = std(I); |
---|
128 | ideal ggI = std(gI); |
---|
129 | ideal gJ = std(J); |
---|
130 | ideal ggJ = std(gJ); |
---|
131 | 0== size( std(reduce (I,gI ) ) ); |
---|
132 | 0== size( std(reduce (J,gJ ) ) ); |
---|
133 | idealsEqual(gI,ggI ); |
---|
134 | idealsEqual(gJ,ggJ ); |
---|
135 | |
---|
136 | |
---|
137 | |
---|
138 | //Github Adi #39 |
---|
139 | //redtail hat Problemen wenn man im posinL11Ring pLmCmp nimmt statt FDeg |
---|
140 | |
---|
141 | |
---|
142 | ring rng = (integer),(x,y,z),(ls(3),C); |
---|
143 | option("redTail"); |
---|
144 | option(redSB); |
---|
145 | ideal 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; |
---|
146 | ideal gI = std(I); |
---|
147 | ideal ggI = std(gI); |
---|
148 | idealsEqual(gI,ggI ); |
---|
149 | gI; |
---|
150 | |
---|
151 | |
---|
152 | //Github Adi #40 |
---|
153 | |
---|
154 | ring rng = (integer),(x,y,z),(ls(3),C); |
---|
155 | ideal 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; |
---|
156 | ideal J = -1+7*y^3*z,-15+14*x+14*x*y*z,15*z+7*x-10*x*z+5*x*y^2; |
---|
157 | ideal gI = std(I); |
---|
158 | |
---|
159 | ring rng = (integer),(x,y,z),(ls(3),C); |
---|
160 | ideal I = 7*x*z+3*x*z^2-x^2,4*z+6*z^2+2*y; |
---|
161 | ideal J = 5-5*z^2+3*y^2*z; |
---|
162 | ideal gI = std(I); |
---|
163 | |
---|
164 | |
---|
165 | |
---|
166 | |
---|
167 | //Github Adi #27 |
---|
168 | |
---|
169 | ring rng = (integer),(y(1..4),x(1..3)),dp; |
---|
170 | ideal 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); |
---|
171 | ideal J = std(I); |
---|
172 | |
---|
173 | |
---|
174 | //Github Adi #20 |
---|
175 | |
---|
176 | ring rng = (integer),(xw,xq,xk),(lp(2),lp(1),C); |
---|
177 | ideal I = -4*xk+10; |
---|
178 | ideal J = -7*xw+4*xk,-13; |
---|
179 | idealsEqual( intersect(I,J) , intersect(std(I),std(J) ) ); |
---|
180 | |
---|
181 | |
---|
182 | //Github Adi #42 |
---|
183 | |
---|
184 | ring rng = (integer),(x,y,z),(dp(1),dp(2),C); |
---|
185 | ideal I = 20*x*z+7*y^2; |
---|
186 | ideal J = 6*x*y+20*z+6,-5*x-16*y-11,8*x*y+2*y^2; |
---|
187 | ideal gI = std(I); |
---|
188 | ideal ggI = std(gI); |
---|
189 | ideal gJ = std(J); |
---|
190 | ideal ggJ = std(gJ); |
---|
191 | std(reduce (I,gI ) ); |
---|
192 | std(reduce (J,gJ ) ); |
---|
193 | size(gJ)== size(ggJ ); |
---|
194 | |
---|
195 | |
---|
196 | |
---|
197 | |
---|
198 | |
---|
199 | //Github Adi #38 |
---|
200 | |
---|
201 | ring rng = (integer),(x(1),x(3),x(4),x(5)),(dp(4),C); |
---|
202 | ideal Q = x(4)^2-x(3)*x(5), |
---|
203 | x(3)*x(4)-x(5), |
---|
204 | x(3)^2-x(4), |
---|
205 | x(1)^2-121*x(5)^2-11*x(5)-1; |
---|
206 | Q = std(Q); |
---|
207 | qring qr = Q; |
---|
208 | ideal a = x(5)-2; |
---|
209 | ideal j = 6,x(5)-2,x(4)+2,x(3)-2,x(1)+3; |
---|
210 | ideal aJJ = quotient(a*j,j); // ? |
---|
211 | aJJ; |
---|
212 | std(aJJ); |
---|
213 | aJJ; |
---|
214 | |
---|
215 | //segfaults |
---|
216 | //ok mit neue posinL11Ring |
---|
217 | |
---|
218 | ring rng = (integer),(x(1),x(2),x(3)),(ls(3),C); |
---|
219 | ideal I = -3*x(1)-10*x(1)*x(3); |
---|
220 | ideal 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); |
---|
221 | ideal gI = std(I); |
---|
222 | ideal ggI = std(gI); |
---|
223 | ideal gJ = std(J); |
---|
224 | ideal ggJ = std(gJ); |
---|
225 | |
---|
226 | |
---|
227 | //ok now |
---|
228 | |
---|
229 | ring rng = (integer),(x(1),x(2)),(ls(2),C); |
---|
230 | ideal I = -41*x(1)-10*x(1)^3*x(2),-22*x(2)^2-12*x(1)^3; |
---|
231 | ideal J = 14,35*x(1)*x(2)+12*x(1)*x(2)^3; |
---|
232 | ideal gI = std(I); |
---|
233 | ideal ggI = std(gI); |
---|
234 | ideal gJ = std(J); |
---|
235 | |
---|
236 | //ok now |
---|
237 | |
---|
238 | ring rng = (integer),(x(1),x(2)),(ls(2),C); |
---|
239 | ideal I = -22-24*x(1)^3,4*x(1)*x(2)^2-19*x(1)^2*x(2); |
---|
240 | ideal J = 16*x(1)*x(2),12-31*x(1)-43*x(1)*x(2); |
---|
241 | ideal gI = std(I); |
---|
242 | ideal ggI = std(gI); |
---|
243 | ideal gJ = std(J); |
---|
244 | |
---|
245 | //in the new posinL11Ring works smooth, in the other one takes forever |
---|
246 | |
---|
247 | ring rng = (integer),(x(1),x(2)),(ls(2),C); |
---|
248 | ideal 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; |
---|
249 | ideal J = -2*x(1)-11*x(1)^2+10*x(1)^2*x(2),-13*x(1)*x(2)-10*x(1)^3; |
---|
250 | ideal gI = std(I); |
---|
251 | |
---|
252 | //ok now |
---|
253 | |
---|
254 | ring rng = (integer),(x(1),x(2)),(ls(2),C); |
---|
255 | ideal I = -10,15*x(1)*x(2)+8*x(1)*x(2)^2; |
---|
256 | ideal J = 2*x(1)^2*x(2),-3+2*x(1)*x(2); |
---|
257 | ideal gI = std(I); |
---|
258 | |
---|
259 | //ok now |
---|
260 | |
---|
261 | ring rng = (integer),(x(1)),(ls(1),C); |
---|
262 | ideal I = -12+2*x(1)^2; |
---|
263 | ideal J = -6*x(1)^2+11*x(1)^4,-15; |
---|
264 | ideal gI = std(I); |
---|
265 | ideal ggI = std(gI); |
---|
266 | ideal gJ = std(J); |
---|
267 | |
---|
268 | //Github Adi #47 |
---|
269 | |
---|
270 | ring R = integer,x,dp; |
---|
271 | reduce(x,std(2*x)); |
---|
272 | |
---|
273 | ring S = integer,x,ls; |
---|
274 | reduce(x,std(2*x)); |
---|
275 | |
---|
276 | //Github Adi #46 |
---|
277 | |
---|
278 | ring rng = integer,x,ds; |
---|
279 | ideal J = -7, -3-15*x^2, -6+18*x-x^2; |
---|
280 | ideal gJ = std(J); |
---|
281 | |
---|
282 | //Github Adi #49 |
---|
283 | |
---|
284 | ring rng = (integer),(x,y,z),(dp(3),C); |
---|
285 | ideal Q = x,yz2+5; |
---|
286 | qring QR = std(Q); |
---|
287 | poly a = z; |
---|
288 | ideal J = 5,z,x; |
---|
289 | ideal aJJ = quotient( a*J, J ); |
---|
290 | aJJ; |
---|
291 | |
---|
292 | //Github Adi #50 |
---|
293 | |
---|
294 | ring rng = (integer),(x(1)),(lp(1),C); |
---|
295 | ideal I = 4*x(1)^3+35*x(1)^2+24*x(1),45*x(1),35; |
---|
296 | ideal J = 18*x(1)^3+43*x(1)^2,-7*x(1)-10,-12*x(1)^2-20*x(1)+13; |
---|
297 | ideal I_J = intersect(I,J); |
---|
298 | size(I_J); |
---|
299 | |
---|
300 | //Github Adi #21 |
---|
301 | |
---|
302 | ring r = integer,(a,d,g,h,i,j),ws(-1,-1,-1,-1,-1,-1); |
---|
303 | ideal inI= gh+2i, dh+j, d+2g+aj, d+2g+ai; |
---|
304 | ideal g= std(inI); |
---|
305 | /// |
---|
306 | ring rng = (integer),(x,y,z),(ws(-1,-1,-1),C); |
---|
307 | ideal I = 6yz-6,6xy-6yz-8y; |
---|
308 | ideal gI =std(I); |
---|
309 | /// |
---|
310 | ring rng = (integer),(x,y,z),(ws(-1,-1,-1),C); |
---|
311 | ideal I = -10x+15y,4y-7; |
---|
312 | ideal gI =std(I); |
---|
313 | /// |
---|
314 | ring rng = (integer),(x,y,z),(ws(-1,-1,-1),C); |
---|
315 | ideal I = 9x-14y,15z-2; |
---|
316 | ideal gI =std(I); |
---|
317 | /// |
---|
318 | ring rng = (integer),(x,y),(ws(-1,-1),C); |
---|
319 | ideal I = 8x+7,-6x-13y; |
---|
320 | ideal gI =std(I); |
---|
321 | |
---|
322 | // Singular bug #648 |
---|
323 | |
---|
324 | ring r=integer,(t,x(1..4)),ws(1,-1,-11,-3,-19); |
---|
325 | option(prot); |
---|
326 | ideal inJ = |
---|
327 | 3, |
---|
328 | t*x(3)*x(4)+2*t*x(1)*x(2), |
---|
329 | t*x(1)*x(2)^2+2*t^3*x(1)*x(2)*x(3), |
---|
330 | x(1)^2*x(2)*x(3)*x(4)+2*t^2*x(1)^3*x(2)*x(3); |
---|
331 | ideal J = |
---|
332 | 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), |
---|
333 | 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), |
---|
334 | t*x(3)*x(4)-3*t^2*x(3)*x(4)-t*x(1)*x(2)-2*x(1)^2, |
---|
335 | 3-t; |
---|
336 | reduce(inJ,J); |
---|
337 | |
---|
338 | //Github Adi #48 |
---|
339 | |
---|
340 | ring rng = integer,x,ds; |
---|
341 | ideal I = -8*x+2*x^2, -16+9*x^2-x^3; |
---|
342 | ideal gI = std(I); |
---|
343 | gI; |
---|
344 | std(gI); |
---|
345 | |
---|
346 | //Github Adi #43 |
---|
347 | |
---|
348 | LIB "primdec.lib"; |
---|
349 | ring rng = integer,(x),lp; |
---|
350 | short = 0 ; |
---|
351 | ideal J = -10*x^2-9*x,-8*x^2-4; |
---|
352 | ideal gJ = std(J); |
---|
353 | ideal ggJ = std(gJ); |
---|
354 | ASSUME(0, 0== size( std(reduce (J,gJ ) )) ); |
---|
355 | ASSUME(0, idealsEqual(gJ,ggJ )); |
---|
356 | gJ; |
---|
357 | ggJ; |
---|
358 | |
---|
359 | ring rng = integer,x,lp; |
---|
360 | option(noredSB) ; |
---|
361 | option(noredTail) ; |
---|
362 | ideal I = -8*x^3+7*x^2-7*x, 8*x^2+10*x-2; |
---|
363 | ideal gI = std(I); |
---|
364 | ideal ggI = std(gI); |
---|
365 | reduce (I,gI ); // <>0 !! |
---|
366 | gI; |
---|
367 | ggI; |
---|
368 | |
---|
369 | LIB("primdec.lib"); |
---|
370 | ring rng = integer,x,lp; |
---|
371 | short = 0 ; |
---|
372 | ideal I = 8*x^3-4,-10*x^3+9*x; |
---|
373 | ideal gI = std(I); |
---|
374 | ideal ggI = std(gI); |
---|
375 | ASSUME(0, 0== size( std(reduce (I,gI ) )) ); |
---|
376 | ASSUME(0, idealsEqual(gI,ggI )); |
---|
377 | size(gI)== size(ggI) ; |
---|
378 | gI; |
---|
379 | ggI; |
---|
380 | |
---|
381 | //Github Adi #25 |
---|
382 | |
---|
383 | ring r = integer,(a,d,g,h,i,j),ws(-1,-1,-1,-1,-1,-1); |
---|
384 | ideal inI= gh+2i, dh+j, d+2g+aj, d+2g+ai; |
---|
385 | ideal g= std(inI); |
---|
386 | inI; |
---|
387 | g; |
---|
388 | |
---|
389 | ring rng = (integer),(x,y,z),(ws(-1,-1,-1),C); |
---|
390 | ideal I = 6yz-6,6xy-6yz-8y; |
---|
391 | ideal gI =std(I); |
---|
392 | I;gI; |
---|
393 | |
---|
394 | |
---|
395 | |
---|
396 | |
---|
397 | tst_status(1);$; |
---|
398 | |
---|