source: git/Tst/Long/deform_l.tst @ 6263f2

spielwiese
Last change on this file since 6263f2 was 6fe3a0, checked in by Hans Schönemann <hannes@…>, 19 years ago
*hannes/lossen: new libs, new tests git-svn-id: file:///usr/local/Singular/svn/trunk@7891 2c84dea3-7e68-4137-9b89-c4e89433aadc
  • Property mode set to 100644
File size: 3.2 KB
Line 
1//deform_l.tst
2//long  tests for deform.lib
3//---------------------------
4LIB "tst.lib";
5tst_init();
6LIB "deform.lib";
7 printlevel =  2;
8ring  r0  = 0,(x,y,z),ls;
9 ideal i2 = y4-xy,yz-xz,z4-xz;   
10 list L=versal(i2);
11kill L;
12//-----------------------------------------------------------------------------
13ring  r1  = 0,(x,y,z,u),dp;
14 ideal i5 = x3-y2,xz,xu,yz,yu,z3-u2;
15 list L=versal(i5);
16 ideal i6 = z2-xz,zu-xu,u2-yu,yz-xu;
17 L=versal(i6); kill L;
18//------- cone over rational normal curve ofe degree d -------------------------
19 int    d = 4;
20ring r2   = 0,(x(1..d+1)),ls;
21 ideal i7 = maxideal(1);
22 matrix m[2][d] = i7[1..d],i7[2..d+1];
23       i7 = minor(m,2);
24       i7 = minbase(i7);
25 list L=versal(i7);
26 def Px=L[1]; setring Px; Fs; print(Js); print(Rs);
27 kill L,Px;
28//------- L_d_d: d generic lines in A^d ---------------------------------------
29        d = 4;
30ring r3=0,(x(1..d)),ls;
31 int n,m; ideal i8;
32 for (n=1;n<d;n=n+1)
33 { for (m=n+1;m<=d;m=m+1) {i8=i8,x(n)*x(m);}
34 }
35        i8 = simplify(i8,10);
36 list L=versal(i8,3,"myring");
37 def Px=L[1]; setring Px; Fs;
38 kill L,Px;
39//------- two elliptic singularity  (t1=7 and 13,t2=3) ------------------------
40ring r6   = 0,(x,y,z,u,v,w),ds;
41 ideal i9 = uv-zw,u2-yw,zu-yv,yu-xw,xu+v2-uw,yz-xv,y2+v2-uw,xy+zv-yw,x2+z2-yu;
42 list L=versal(i9);
43 def Px=L[1]; setring Px; Fs;
44 kill L,Px;
45ring r7   = 0,(x,y,z,u,v),dp;
46 ideal I1 = xy,xz,xu,xv,yz,yu,yv,zu,zv,uv;   
47 list L=versal(I1);
48 def Px=L[1]; setring Px; Fs;
49 kill L,Px;
50
51ring r8   = 0,(w(1..9)),dp;
52 ideal I3 = w(8)^2-w(6)*w(9),
53            w(7)*w(8)-w(5)*w(9),
54            w(5)*w(8)-w(3)*w(9),
55            w(4)*w(8)-w(2)*w(9),
56            w(7)^2-w(4)*w(9),
57            w(6)*w(7)-w(3)*w(9),
58            w(5)*w(7)-w(2)*w(9),
59            w(4)*w(7)-w(1)*w(9),
60            w(3)*w(7)-w(2)*w(8),
61            w(2)*w(7)-w(1)*w(8),
62            w(5)*w(6)-w(3)*w(8),
63            w(4)*w(6)-w(2)*w(8),
64            w(5)^2-w(3)*w(7),
65            w(4)*w(5)-w(2)*w(7),
66            w(3)*w(5)-w(2)*w(6),
67            w(2)*w(5)-w(1)*w(6),
68            w(4)^2-w(1)*w(7),
69            w(3)*w(4)-w(2)*w(5),
70            w(2)*w(4)-w(1)*w(5),
71            w(2)^2-w(1)*w(3),
72            w(2)*w(6)*w(8)-w(3)^2*w(9),
73            w(1)*w(6)*w(8)-w(2)*w(3)*w(9),
74            w(2)*w(6)^2-w(3)^2*w(8),
75            w(1)*w(6)^2-w(2)*w(3)*w(8);
76list L= versal(I3);
77def Px=L[1]; setring Px; Fs;
78kill L,Px;
79kill r0,r1,r2,r3,r6,r7,r8,m,n,d;
80//////////////////////////////////////////////////////////////////////////////////
81ring   Po = 0,(x,y),dp;
82ideal  Io = std(x^4+y^3);
83matrix Mo;
84//============= rk 1 ======================================
85module k =[x],[y];
86module m(0)=k;
87module m(1)=[x,y],[-y2,x3];
88module m(2)=[x3,y],[-y2,x];
89module m(3)=[x2,y],[-y2,x2];
90module m(4)=[y,x,0],[-x2,0,-y],[0,-y,-x];
91//============= rk 2 ======================================
92module m(5)=[xy,-y2,x2],[-y2,-x3,-xy],[-x3,x2y,y2];
93module m(6)=[xy,-yy,xx,0],[-yy,-xxx,-xy,0],[-xx,0,0,-y],[0,x2,y,x];
94int i';
95int ii=6;
96for (i'=4;i'<ii;i'=i'+1)
97{
98  Mo=m(i');
99  list L=mod_versal(Mo,Io);
100  def Qx=L[2];
101  print(Ms);
102  print(Ls);
103  print(Js);
104  kill L,Qx;
105}
106list L=mod_versal(m(2),Io,4,"qq","W","dp","oo");
107def Px=L[1]; def Qx=L[2]; def So=L[3]; def Ox=L[4];
108listvar(ring);
109show(Qx);
110tst_status(1);$
Note: See TracBrowser for help on using the repository browser.