source: git/Tst/Short/equising.tst @ 379305

spielwiese
Last change on this file since 379305 was 379305, checked in by Hans Schönemann <hannes@…>, 23 years ago
*hannes: added equising examples git-svn-id: file:///usr/local/Singular/svn/trunk@4783 2c84dea3-7e68-4137-9b89-c4e89433aadc
  • Property mode set to 100644
File size: 1.7 KB
Line 
1LIB "tst.lib";
2tst_init();
3
4LIB "sing.lib";
5LIB "deform.lib";
6LIB "equising.lib";
7
8///////////////////////////////////////////////////////////////////////////////////////////////
9
10proc 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
42list 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
53tst(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/////////////////////////////////////////////////
67ring rr = 11,(A,B,C,a,b,x,y),(dp(3),dp(2),dp);
68
69poly F = (y-ax-b*x^2)^5-x^12;
70ideal id = (1-A)*a-1,(1+B+A^2)*b-1;
71F = F*(1-A)^5*(1+B+A^2)^5;
72F = reduce(F,std(id));
73
74esStratum(F);
75
76isEquising(F);
77
78example esStratum;
79example isEquising;
80
81tst_status(1);$
Note: See TracBrowser for help on using the repository browser.