1 | //deform_l.tst |
---|
2 | //long tests for deform.lib |
---|
3 | //--------------------------- |
---|
4 | LIB "tst.lib"; |
---|
5 | tst_init(); |
---|
6 | LIB "deform.lib"; |
---|
7 | printlevel = 2; |
---|
8 | ring r0 = 0,(x,y,z),ls; |
---|
9 | ideal i2 = y4-xy,yz-xz,z4-xz; |
---|
10 | versal(i2);kill_rings(); |
---|
11 | //------------------------------------------------------------------------------ |
---|
12 | ring r1 = 0,(x,y,z,u),dp; |
---|
13 | ideal i5 = x3-y2,xz,xu,yz,yu,z3-u2; |
---|
14 | versal(i5);kill_rings(); |
---|
15 | ideal i6 = z2-xz,zu-xu,u2-yu,yz-xu; |
---|
16 | versal(i6);kill_rings(); |
---|
17 | //------- cone over rational normal curve ofe degree d ------------------------- |
---|
18 | int d = 4; |
---|
19 | ring r2 = 0,(x(1..d+1)),ls; |
---|
20 | ideal i7 = maxideal(1); |
---|
21 | matrix m[2][d] = i7[1..d],i7[2..d+1]; |
---|
22 | i7 = minor(m,2); |
---|
23 | i7 = minbase(i7); |
---|
24 | versal(i7);kill_rings(); |
---|
25 | //------- L_d_d: d generic lines in A^d ________________________________________ |
---|
26 | d = 4; |
---|
27 | ring r3=0,(x(1..d)),ls; |
---|
28 | int n,m; ideal i8; |
---|
29 | for (n=1;n<d;n=n+1) |
---|
30 | { for (m=n+1;m<=d;m=m+1) {i8=i8,x(n)*x(m);} |
---|
31 | } |
---|
32 | i8 = simplify(i8,10); |
---|
33 | versal(i8,3,"myring");kill_rings("myring"); |
---|
34 | //------- two elliptic singularity (t1=7 and 13,t2=3) --------------------------- |
---|
35 | ring r6 = 0,(x,y,z,u,v,w),ds; |
---|
36 | ideal i9 = uv-zw,u2-yw,zu-yv,yu-xw,xu+v2-uw,yz-xv,y2+v2-uw,xy+zv-yw,x2+z2-yu; |
---|
37 | versal(i9);kill_rings(); |
---|
38 | ring r7 = 0,(x,y,z,u,v),dp; |
---|
39 | ideal I1 = xy,xz,xu,xv,yz,yu,yv,zu,zv,uv; |
---|
40 | versal(I1);kill_rings(); |
---|
41 | ring r8 = 0,(w(1..9)),dp; |
---|
42 | ideal I3 = w(8)^2-w(6)*w(9), |
---|
43 | w(7)*w(8)-w(5)*w(9), |
---|
44 | w(5)*w(8)-w(3)*w(9), |
---|
45 | w(4)*w(8)-w(2)*w(9), |
---|
46 | w(7)^2-w(4)*w(9), |
---|
47 | w(6)*w(7)-w(3)*w(9), |
---|
48 | w(5)*w(7)-w(2)*w(9), |
---|
49 | w(4)*w(7)-w(1)*w(9), |
---|
50 | w(3)*w(7)-w(2)*w(8), |
---|
51 | w(2)*w(7)-w(1)*w(8), |
---|
52 | w(5)*w(6)-w(3)*w(8), |
---|
53 | w(4)*w(6)-w(2)*w(8), |
---|
54 | w(5)^2-w(3)*w(7), |
---|
55 | w(4)*w(5)-w(2)*w(7), |
---|
56 | w(3)*w(5)-w(2)*w(6), |
---|
57 | w(2)*w(5)-w(1)*w(6), |
---|
58 | w(4)^2-w(1)*w(7), |
---|
59 | w(3)*w(4)-w(2)*w(5), |
---|
60 | w(2)*w(4)-w(1)*w(5), |
---|
61 | w(2)^2-w(1)*w(3), |
---|
62 | w(2)*w(6)*w(8)-w(3)^2*w(9), |
---|
63 | w(1)*w(6)*w(8)-w(2)*w(3)*w(9), |
---|
64 | w(2)*w(6)^2-w(3)^2*w(8), |
---|
65 | w(1)*w(6)^2-w(2)*w(3)*w(8); |
---|
66 | versal(I3);kill_rings(); |
---|
67 | kill r0,r1,r2,r3,r6,r7,r8; |
---|
68 | ////////////////////////////////////////////////////////////////////////////////// |
---|
69 | ring Po = 0,(x,y),dp; |
---|
70 | ideal Io = std(x^4+y^3); |
---|
71 | matrix Mo; |
---|
72 | //============= rk 1 ====================================== |
---|
73 | module k =[x],[y]; |
---|
74 | module m(0)=k; |
---|
75 | module m(1)=[x,y],[-y2,x3]; |
---|
76 | module m(2)=[x3,y],[-y2,x]; |
---|
77 | module m(3)=[x2,y],[-y2,x2]; |
---|
78 | module m(4)=[y,x,0],[-x2,0,-y],[0,-y,-x]; |
---|
79 | //============= rk 2 ====================================== |
---|
80 | module m(5)=[xy,-y2,x2],[-y2,-x3,-xy],[-x3,x2y,y2]; |
---|
81 | module m(6)=[xy,-yy,xx,0],[-yy,-xxx,-xy,0],[-xx,0,0,-y],[0,x2,y,x]; |
---|
82 | int i'; |
---|
83 | int ii=6; |
---|
84 | for (i'=4;i'<ii;i'=i'+1) |
---|
85 | { |
---|
86 | Mo=m(i'); |
---|
87 | mod_versal(Mo,Io); |
---|
88 | kill_rings(); |
---|
89 | } |
---|
90 | mod_versal(m(2),Io,4,"qq","W","dp","oo"); |
---|
91 | listvar(ring); |
---|
92 | qqQx; |
---|
93 | tst_status(1);$ |
---|