1 | LIB "tst.lib"; |
---|
2 | tst_init(); |
---|
3 | |
---|
4 | LIB "sing.lib"; |
---|
5 | LIB "deform.lib"; |
---|
6 | LIB "equising.lib"; |
---|
7 | |
---|
8 | /////////////////////////////////////////////////////////////////////////////////////////////// |
---|
9 | |
---|
10 | proc tst(list d) |
---|
11 | { |
---|
12 | int i; |
---|
13 | for(i=1;i<=size(d);i++) |
---|
14 | { |
---|
15 | execute(d[i]); |
---|
16 | |
---|
17 | poly f = y8+4x3y6+6x6y4+2x5y5+4x9y2+4x8y3+x12+2x11y+2x10y2+x13; |
---|
18 | |
---|
19 | versal(f); |
---|
20 | setring Px; |
---|
21 | poly F = Fs[1,1]; |
---|
22 | |
---|
23 | list EsList = esStratum(F); |
---|
24 | EsList; |
---|
25 | ideal ES = EsList[1]; |
---|
26 | |
---|
27 | qring q = std(ES); |
---|
28 | poly F = imap(Px,F); |
---|
29 | |
---|
30 | isEquising(F); |
---|
31 | |
---|
32 | kill(Px); |
---|
33 | kill(Qx); |
---|
34 | kill(Ox); |
---|
35 | kill(So); |
---|
36 | kill(q); |
---|
37 | kill(r); |
---|
38 | } |
---|
39 | } |
---|
40 | |
---|
41 | |
---|
42 | list l = |
---|
43 | "ring r = 0,(x,y),ds", |
---|
44 | "ring r = (0,a),(x,y),ds; minpoly = a^2-1;", |
---|
45 | "ring r = (0, t(1..4)),(x,y),ds", |
---|
46 | "ring r = 7,(x,y),ds", |
---|
47 | "ring r = (11,a(1)),(x,y),ds; minpoly = a(1)^2-1;", |
---|
48 | "ring r = (3,t(1..8)),(x,y),ds", |
---|
49 | "ring r = (27,a),(x,y),ds" |
---|
50 | ; |
---|
51 | |
---|
52 | |
---|
53 | tst(l); |
---|
54 | |
---|
55 | ////////////////////////////////////////////////////// |
---|
56 | // Die folgenden beiden Beispiele sind "aquisingul"ar. |
---|
57 | |
---|
58 | ring r = 0,(A,B,C,D,x,y),dp; |
---|
59 | poly F =D3x12+3CD2x11+3C2Dx10+3BD2x10+x13+C3x9+6BCDx9+3AD2x9+3BC2x8+3B2Dx8+6ACDx8+3D2x9-3D2x8y+3B2Cx7+3AC2x7+6ABDx7+6CDx8-6CDx7y+B3x6+6ABCx6+3A2Dx6+3C2x7+6BDx7-3C2x6y-6BDx6y+3AB2x5+3A2Cx5+6BCx6+6ADx6-6BCx5y-6ADx5y+3A2Bx4+3B2x5+6ACx5+3Dx6-3B2x4y-6ACx4y-6Dx5y+3Dx4y2+A3x3+6ABx4+3Cx5-6ABx3y-6Cx4y+3Cx3y2+3A2x3+3Bx4-3A2x2y-6Bx3y+3Bx2y2+3Ax3-6Ax2y+3Axy2+x3-3x2y+3xy2-y3; |
---|
60 | |
---|
61 | esStratum(F); |
---|
62 | |
---|
63 | isEquising(F); |
---|
64 | |
---|
65 | |
---|
66 | ///////////////////////////////////////////////// |
---|
67 | ring rr = 11,(A,B,C,a,b,x,y),(dp(3),dp(2),dp); |
---|
68 | |
---|
69 | poly F = (y-ax-b*x^2)^5-x^12; |
---|
70 | ideal id = (1-A)*a-1,(1+B+A^2)*b-1; |
---|
71 | F = F*(1-A)^5*(1+B+A^2)^5; |
---|
72 | F = reduce(F,std(id)); |
---|
73 | |
---|
74 | esStratum(F); |
---|
75 | |
---|
76 | isEquising(F); |
---|
77 | |
---|
78 | example esStratum; |
---|
79 | example isEquising; |
---|
80 | |
---|
81 | tst_status(1);$ |
---|