jengelh-datetimespielwiese
Last change on this file since ad3dfe was ad3dfe, checked in by Hans Schönemann <hannes@…>, 14 years ago
*hannes: new zeroset.lib git-svn-id: file:///usr/local/Singular/svn/trunk@11689 2c84dea3-7e68-4137-9b89-c4e89433aadc
• Property mode set to `100644`
File size: 3.2 KB
Line
1// Test File for 'zeroset.lib"
2//
3// Implementation by : Thomas Bayer
5// Institut fuer Informatik, Technische Universitaet Muenchen
6// www:    http://wwwmayr.informatik.tu-muenchen.de/personen/bayert/
7// email : bayert@in.tum.de
8//
9// Last change 10.12.2000
10//
11// written in the frame of the diploma thesis (advisor: Prof. Gert-Martin Greuel)
12// "Computations of moduli spaces of semiquasihomogenous singularities and an
13//  implementation in Singular"
14// Arbeitsgruppe Algebraische Geometrie, Fachbereich Mathematik,
15// Universitaet Kaiserslautern
16///////////////////////////////////////////////////////////////////////////////
17
18LIB "tst.lib";
19LIB "zeroset.lib";
20
21tst_init();
22tst_ignore("CVS ID \$Id: zeroset1.tst,v 1.2 2009-04-14 10:29:45 Singular Exp \$");
23
24
25
26proc TestzerosetLIB(S)
27// plug in the elements from the zero-set in the generators of the ideal
28// result should be 0,0,...,0.
29{
30
31        setring S;
32        map F;
33
34        print(" #Points = " + string(vdim(std(id))) + "; #Solutions = " + string(size(theZeroset)));
35        print(" minpoly = " + string(minpoly));
36        print(" ideal = " + string(id));
37
38
39        for(int i = 1; i <= size(theZeroset); i++) {
40                F = theZeroset[i];
41                print(string(theZeroset[i]) + " --> " + string(F(id)));
42        }
43}
44
45// example 1
46
47ring R = 0,(x,y,z), lp;
48ideal I = x2-1,y3-1,z-2x;
49
50def S = zeroSet(I);TestzerosetLIB(S);
51setring S;
52minpoly;
53id;
54theZeroset;
55
56kill R;
57kill S;
58//tst_status();
59
60// example 2
61
62ring R = 0,(x,y,z), lp;
63ideal I = x3-1,3y-x,z4-1;
64
65def S = zeroSet(I);TestzerosetLIB(S);
66setring S;
67minpoly;
68id;
69theZeroset;
70
71kill R;
72kill S;
73//tst_status();
74
75// example 3
76
77ring R = 0, (x(1..5)), lp;
78ideal I = x(1)^2+1,x(2)^4-1,x(3)-x(4),x(4)^2+1,x(5)-x(4);
79
80def S = zeroSet(I);TestzerosetLIB(S);
81setring S;
82minpoly;
83id;
84theZeroset;
85
86kill R;
87kill S;
88//tst_status();
89
90// example 4
91
92ring R = (0,a),(x,y), lp;
93minpoly = a2+1;
94ideal I = x2-a,y-2x;
95
96def S = zeroSet(I);TestzerosetLIB(S);
97setring S;
98minpoly;
99id;
100theZeroset;
101
102kill R;
103kill S;
104//tst_status();
105
106// example 5
107
108ring R = (0,a),(x,y,z), lp;
109minpoly = a2+a+1;
110ideal I = x3-1,a*y^2-x,z^2-a*y;
111
112def S = zeroSet(I);TestzerosetLIB(S);
113setring S;
114minpoly;
115id;
116theZeroset;
117
118kill R;
119kill S;
120//tst_status();
121
122// example 6
123
124ring R = (0,a),(x,y,z), lp;
125minpoly = a2+1;
126ideal I = x2-2,a*y^2-2*x,z^2 - a; // 2
127
128def S = zeroSet(I);TestzerosetLIB(S);
129setring S;
130minpoly;
131id;
132theZeroset;
133
134kill R;
135kill S;
136//tst_status();
137
138// example 7
139
140ring R = (0,a),(x(1..5)), lp;
141minpoly = a2+a+1;
142ideal I = x(1)^2 - 2,x(2) -a*x(1), x(3)^2 + 1, x(4)^2 + a, a*x(5) - x(1);
143
144def S = zeroSet(I);TestzerosetLIB(S);
145setring S;
146minpoly;
147id;
148theZeroset;
149
150kill S;
151
152setring R;
153def S = zeroSet(I, 1);TestzerosetLIB(S);
154setring S;
155minpoly;
156id;
157theZeroset;
158
159kill R;
160kill S;
161
162//tst_status();
163
164// example 8
165
166ring R = (0,a),(s(1..4)), lp;
167minpoly = a2+1;
168ideal I = s(3)*s(4),s(2)*s(4),s(1)*s(3),s(1)*s(2),s(3)^8+s(4)^8-1,s(2)^8+s(4)^8-1,s(1)^8-s(4)^8,s(4)^9-s(4),s(1)*s(4)^8-s(1);
169
170def S = zeroSet(I);TestzerosetLIB(S);
171setring S;
172minpoly;
173id;
174theZeroset;
175
176kill S;
177
178setring R;
179def S = zeroSet(I, 1);TestzerosetLIB(S);
180setring S;
181minpoly;
182id;
183theZeroset;
184
185kill S;
186
187tst_status(1); \$
Note: See TracBrowser for help on using the repository browser.