1 | // $Id: makedbm.lib,v 1.4 1997-11-10 19:18:46 krueger Exp $ |
---|
2 | //========================================================================= |
---|
3 | // |
---|
4 | // Please send bugs and comments to krueger@mathematik.uni-kl.de |
---|
5 | // |
---|
6 | //============================================================================= |
---|
7 | |
---|
8 | LIBRARY: makedbm.lib some usefull tools needed by the Arnold-Classifier. |
---|
9 | |
---|
10 | dbm_read(l); read all entries from a DBM-databaes pointed by l |
---|
11 | dbm_getnext(l); read next entry from a DBM-databaes pointed by l |
---|
12 | |
---|
13 | //============================================================================= |
---|
14 | |
---|
15 | proc dbm_read (link l) |
---|
16 | { |
---|
17 | string s=""; |
---|
18 | s=read(l); |
---|
19 | while( s != "" ) |
---|
20 | { |
---|
21 | s,"=",read(l,s); |
---|
22 | s=read(l); |
---|
23 | } |
---|
24 | } |
---|
25 | |
---|
26 | //============================================================================= |
---|
27 | proc dbm_getnext (link l) |
---|
28 | { |
---|
29 | string s=""; |
---|
30 | s=read(l); |
---|
31 | if( s != "" ) { s,"=",read(l,s); } |
---|
32 | } |
---|
33 | |
---|
34 | //============================================================================= |
---|
35 | proc create_sing_dbm |
---|
36 | { |
---|
37 | link l="DBM:rw NFlist"; |
---|
38 | write(l, "A[k]", "x^(k+1)"); |
---|
39 | write(l, "D[k]", "x2y+y^(k-1)"); |
---|
40 | write(l, "E[6k]", "x3+y^(3*k+1)+a*x*(y^(2*k+1))"); |
---|
41 | write(l, "E[6k+1]", "x3+x*(y^(2*k+1))+a*(y^(3*k+2))"); |
---|
42 | write(l, "E[6k+2]", "x3+y^(3*k+2)+a*x*(y^(2*k+2))"); |
---|
43 | write(l, "J[k,0]", "x3+b*x^2*y^k+y^(3*k)+c*x*y^(2*k+1)"); |
---|
44 | write(l, "J[k,r]", "x3+x2*y^k+a*y^(3*k+r)"); |
---|
45 | write(l, "X[1,0]", "x4+a*x2y2+y4"); |
---|
46 | write(l, "X[1,r]", "x4+x2y2+a*y^(4+r)"); |
---|
47 | write(l, "X[k,0]", "x4+b*x3y^k+a*x2y^(2*k) + xy^(3*k)"); |
---|
48 | write(l, "X[k,r]", "x4+a*x3*y^k+x2y^(2*k)+b*(y^(4*k+r))"); |
---|
49 | write(l, "W[12k]", "x4+y^(4*k+1)+a*x*(y^(3*k+1))+c*x2*(y^(2*k+1))"); |
---|
50 | write(l, "W[12k+1]", "x4+x*(y^(3*k+1))+a*x2*(y^(2*k+1))+c*y^(4*k+2)"); |
---|
51 | write(l, "W[12k+5]", "x4+x*(y^(3*k+2))+a*x2*(y^(2*k+2))+b*y^(4*k+3)"); |
---|
52 | write(l, "W[12k+6]", "x4+y^(4*k+3)+a*x*(y^(3*k+3))+b*x2*(y^(2*k+2))"); |
---|
53 | write(l, "W[k,0]", "x4+b*x2*(y^(2*k+1))+a*x*(y^(3*k+2))+y^(4*k+2)"); |
---|
54 | write(l, "W[k,r]", "x4+a*x3*(y^(k+1))+x2*(y^(2*k+1))+b*y^(4*k+2+r)"); |
---|
55 | write(l, "W#[k,2r-1]", "(x2+y^(2*k+1))^2+b*x*(y^(3*k+1+r))+a*y^(4*k+2+r)"); |
---|
56 | write(l, "W#[k,2r]", "(x2+y^(2*k+1))^2+b*x2*(y^(2*k+1+r))+a*x*(y^(3*k+2+r))"); |
---|
57 | write(l, "Y[1,r,s]", "x^(4+r)+a*x2*(y2)+y^(4+s)"); |
---|
58 | write(l, "Y[k,r,s]", "((x+a*y^k)^2 + b*y^(2*k+s))*( x2 + y^(2*k+r))"); |
---|
59 | write(l, "Z[1,0]", "x3y + x2y3 + xy6 +y7"); |
---|
60 | write(l, "Z[1,r]", "x3y + x2*(y^3) + a*y^(7+r)"); |
---|
61 | write(l, "Z[k,r,s]", "(x2+a*x*(y^k)+b*y^(2*k+r))*(x2+y^(2*k+2*r+s))"); |
---|
62 | write(l, "Z[k,r,0]", "(x+a*(y^k))*(x3+d*x2*(y^(k+r))+c*x*(y^(2*k+2*r+1))+y^(3*k+3*r))"); |
---|
63 | write(l, "Z[k,12k+6r-1]", "(x+a*(y^k))*(x3+b*x*(y^(2*k+2*r+1))+y^(3*k+3*r+1))"); |
---|
64 | write(l, "Z[k,12k+6r]", "(x+a*(y^k))*(x3+x*(y^(2*k+2*r+1))+b*y^(3*k+3*r+2))"); |
---|
65 | write(l, "Z[k,12k+6r+1]", "(x+a*(y^k))*(x3+b*x*(y^(2*k+2*r+2))+y^(3*k+3*r+2))"); |
---|
66 | write(l, "Z[k,0]", "y*(x3+d*x2*(y^k)+c*x*(y^(2*k+1))+y^(3*k))"); |
---|
67 | write(l, "Z[k,r]", "y*(x3+x2y^(k+1)+b*(y^(3*k+r+3)))"); |
---|
68 | write(l, "Z[6k+5]", "y*(x3+b*x*(y^(2*k+1))+y^(3*k+1))"); |
---|
69 | write(l, "Z[6k+6]", "y*(x3+x*(y^(2*k+1))+b*y^(3*k+2))"); |
---|
70 | write(l, "Z[6k+7]", "y*(x3+b*x*(y^(2*k+2))+y^(3*k+2))"); |
---|
71 | write(l, "Q[k,0]", "x3+z2y+b*x2*(y^k)+x*(y^(2*k))"); |
---|
72 | write(l, "Q[k,r]", "x3+z2y+x2*(y^k)+b*y^(3*k+r)"); |
---|
73 | write(l, "Q[6k+4]", "x3+z2y+y^(3*k+1)+b*x*(y^(2*k+1))"); |
---|
74 | write(l, "Q[6k+5]", "x3+z2y+x*(y^(2*k+1))+b*y^(3*k+2)"); |
---|
75 | write(l, "Q[6k+6]", "x3+z2y+y^(3*k+2)+b*x*(y^(2*k+2))"); |
---|
76 | write(l, "S[12k-1]", "x2z+z2y+y^(4*k)+a*x*(y^(3*k))+c*z*(y^(2*k+1))"); |
---|
77 | write(l, "S[12k]", "x2z+z2y+x*(y^(3*k))+c*y^(4*k+1)+a*z*(y^(2*k+1))"); |
---|
78 | write(l, "S[k,0]", "x2z+z2y+y^(4*k+1)+a*x*(y^(3*k+1))+b*z*(y^(2*k+1))"); |
---|
79 | write(l, "S[k,r]", "x2z+z2y+x2*(y^(2*k))+a*x3*(y^k)+b*y^(4*k+r+1)"); |
---|
80 | write(l, "S#[k,2r-1]", "x2z+z2y+z*(y^(2*k+1))+b*x*(y^(3*k+r))+a*(y^(4*k+r+1))"); |
---|
81 | write(l, "S#[k,2r]", "x2z+z2y+z*(y^(2*k+1))+b*x2*(y^(2*k+r))"); |
---|
82 | write(l, "S[12k+4]", "x2z+z2y+x*(y^(3*k+1))+a*z*(y^(2*k+2))+b*y^(4*k+2)"); |
---|
83 | write(l, "S[12k+5]", "x2z+z2y+y^(4*k+2)+a*x*(y^(3*k+2))+b*z*(y^(2*k+2))"); |
---|
84 | write(l, "U[12k]", "x3+z2x+y^(3*k+1)+a*x*(y^(2*k+1))+b*z*(y^(2*k+1))+d*x2*(y^(k+1))"); |
---|
85 | write(l, "U[k,2r-1]", "x3+z2x+x*(y^(2*k+1))+a*x2*(y^(k+1))+b*(y^(3*k+r+2))+c*z*(y^(2*k+r+1))"); |
---|
86 | write(l, "U[k,2r]", "x3+z2x+x*(y^(2*k+1))+a*x2*(y^(k+1))+b*z*(y^(2*k+r+1))+c*z2*(y^(k+r))"); |
---|
87 | write(l, "U[12k+4]", "x3+z2x+y^(3*k+2)+a*x*(y^(2*k+2))+b*z*(y^(2*k+2))+c*x2*(y^(k+1))"); |
---|
88 | write(l, "V[1,0]", "x2y+z4+a*z3y+b*z2y2+y3z"); |
---|
89 | write(l, "V[1,r]", "x2y+z4+b*z3y+z2y2+a*(y^(r+4))"); |
---|
90 | write(l, "V#[1,2r-1]", "x2y+z3y+a*z2y2+y4+b*x*(z^(r+2))"); |
---|
91 | write(l, "V#[1,2r]", "x2y+z3y+a*z2y2+y4+b*(z^(r+4))"); |
---|
92 | write(l, "T[k,r,s]", "x^k+y^r+z^s+xyz"); |
---|
93 | |
---|
94 | ////////////////////////////////////////////////////////////////// |
---|
95 | // DatenFormat: crk=#; Mu=#; MlrCd=#; |
---|
96 | string s; |
---|
97 | s ="crk=1; Mu=k; MlnCd=k;"; |
---|
98 | write(l, "I_A[k]", s); |
---|
99 | s = "crk=2; Mu=k; MlnCd=1,1,k-3"; |
---|
100 | write(l, "I_D[k]", s); |
---|
101 | s = "crk=2; Mu=6*k; MlnCd=1,2*k,2*k-1"; |
---|
102 | write(l, "I_E[6k]", s); |
---|
103 | s = "crk=2; Mu=6*k+1; MlnCd=1,2*k,2*k"; |
---|
104 | write(l, "I_E[6k+1]", s); |
---|
105 | //"I_E[6k+1]=", read(l, "I_E[6k+1]"); |
---|
106 | s = "crk=2; Mu=6*k+2; MlnCd=1,2*k+1,2*k-1"; |
---|
107 | write(l, "I_E[6k+2]", s); |
---|
108 | s = "crk=2; Mu=6*k-2; MlnCd=1,2*k-1,2*k-1"; |
---|
109 | write(l, "I_J[k,0]", s); |
---|
110 | //"I_J[k,0]=", read(l, "I_J[k,0]"); |
---|
111 | s = "crk=2; Mu=6*k-2+r; MlnCd=1,2*k-1,2*k-1+r"; |
---|
112 | write(l, "I_J[k,r]", s); |
---|
113 | //"I_J[k,r]=", read(l, "I_J[k,r]"); |
---|
114 | s = "crk=2; Mu=9; MlnCd=1,1,1,1,1"; |
---|
115 | write(l, "I_X[1,0]", s); |
---|
116 | s = "crk=2; Mu=9+r; MlnCd=1,1,1,1,1+r"; |
---|
117 | write(l, "I_X[1,r]", s); |
---|
118 | s = "crk=2; Mu=12*k-3; MlnCd=1,1,2*k-1,2*k-1,2*k-1"; |
---|
119 | write(l, "I_X[k,0]", s); |
---|
120 | s = "crk=2; Mu=12*k-3+r; MlnCd=1,1,2*k-1,2*k-1,2*k-1+r"; |
---|
121 | write(l, "I_X[k,r]", s); |
---|
122 | s = "crk=2; Mu=12*k;"; |
---|
123 | write(l, "I_W[12k]", s); |
---|
124 | s = "crk=2; Mu=12*k+1;"; |
---|
125 | write(l, "I_W[12k+1]", s); |
---|
126 | s = "crk=2; Mu=12*k+5;"; |
---|
127 | write(l, "I_W[12k+5]", s); |
---|
128 | s = "crk=2; Mu=12*k+6;"; |
---|
129 | write(l, "I_W[12k+6]", s); |
---|
130 | s = "crk=2; Mu=12*k+3;"; |
---|
131 | write(l, "I_W[k,0]", s); |
---|
132 | s = "crk=2; Mu=12*k+3+r;"; |
---|
133 | write(l, "I_W[k,r]", s); |
---|
134 | s = "crk=2; Mu=12*k+2+2*r;"; |
---|
135 | write(l, "I_W#[k,2r-1]", s); |
---|
136 | s = "crk=2; Mu=12*k+3+2*r;"; |
---|
137 | write(l, "I_W#[k,2r]", s); |
---|
138 | s = "crk=2; Mu=9+r+s;"; |
---|
139 | write(l, "I_Y[1,r,s]", s); |
---|
140 | s = "crk=2; Mu=12*k-3+r+s;"; |
---|
141 | write(l, "I_Y[k,r,s]", s); |
---|
142 | s = "crk=2; Mu=15;"; |
---|
143 | write(l, "I_Z[1,0]", s); |
---|
144 | s = "crk=2; Mu=15+r;"; |
---|
145 | write(l, "I_Z[1,r]", s); |
---|
146 | s = "crk=2; Mu=9+6*k+r;"; |
---|
147 | write(l, "I_Z[k,r]", s); |
---|
148 | s = "crk=2; Mu=12*k+6*r-3;"; |
---|
149 | write(l, "I_Z[k,r,0]", s); |
---|
150 | s = "crk=2; Mu=12*k+6*r+s-3;"; |
---|
151 | write(l, "I_Z[k,r,s]", s); |
---|
152 | s = "crk=2; Mu=12*k+6*r-1;"; |
---|
153 | write(l, "I_Z[k,12k+6r-1]", s); |
---|
154 | s = "crk=2; Mu=12*k+6*r;"; |
---|
155 | write(l, "I_Z[k,12k+6r]", s); |
---|
156 | s = "crk=2; Mu=12*k+6*r+1;"; |
---|
157 | write(l, "I_Z[k,12k+6r+1]", s); |
---|
158 | s = "crk=2; Mu=9+6*k;"; |
---|
159 | write(l, "I_Z[k,0]", s); |
---|
160 | s = "crk=2; Mu=6*(r+1)-1;"; |
---|
161 | write(l, "I_Z[6k+5]", s); |
---|
162 | s = "crk=2; Mu=6*(r+1);"; |
---|
163 | write(l, "I_Z[6k+6]", s); |
---|
164 | s = "crk=2; Mu=6*(r+1)+1;"; |
---|
165 | write(l, "I_Z[6k+7]", s); |
---|
166 | s = "crk=3; Mu=6*k+2;"; |
---|
167 | write(l, "I_Q[k,0]", s); |
---|
168 | s = "crk=3; Mu=6*k+2+r;"; |
---|
169 | write(l, "I_Q[k,r]", s); |
---|
170 | s = "crk=3; Mu=6*k+4;"; |
---|
171 | write(l, "I_Q[6k+4]", s); |
---|
172 | s = "crk=3; Mu=6*k+5;"; |
---|
173 | write(l, "I_Q[6k+5]", s); |
---|
174 | s = "crk=3; Mu=6*k+6;"; |
---|
175 | write(l, "I_Q[6k+6]", s); |
---|
176 | s = "crk=3; Mu=12*k-1;"; |
---|
177 | write(l, "I_S[12k-1]", s); |
---|
178 | s = "crk=3; Mu=12*k;"; |
---|
179 | write(l, "I_S[12k]", s); |
---|
180 | s = "crk=3; Mu=12*k+2;"; |
---|
181 | write(l, "I_S[k,0]", s); |
---|
182 | s = "crk=3; Mu=12*k+2+r;"; |
---|
183 | write(l, "I_S[k,r]", s); |
---|
184 | s = "crk=3; Mu=12*k+2*r+1;"; |
---|
185 | write(l, "I_S#[k,2r-1]", s); |
---|
186 | s = "crk=3; Mu=12*k+2*r+2;"; |
---|
187 | write(l, "I_S#[k,2r]", s); |
---|
188 | s = "crk=3; Mu=12*k+4;"; |
---|
189 | write(l, "I_S[12k+4]", s); |
---|
190 | s = "crk=3; Mu=12*k+5;"; |
---|
191 | write(l, "I_S[12k+5]", s); |
---|
192 | s = "crk=3; Mu=12*k;"; |
---|
193 | write(l, "I_U[12k]", s); |
---|
194 | s = "crk=3; Mu=12*k+4;"; |
---|
195 | write(l, "I_U[12k+4]", s); |
---|
196 | s = "crk=3; Mu=12*k+1+2*r;"; |
---|
197 | write(l, "I_U[k,2r-1]", s); |
---|
198 | s = "crk=3; Mu=12*k+2+2*r;"; |
---|
199 | write(l, "I_U[k,2r]", s); |
---|
200 | s = "crk=3; Mu=15;"; |
---|
201 | write(l, "I_V[1,0]", s); |
---|
202 | s = "crk=3; Mu=15+r;"; |
---|
203 | write(l, "I_V[1,r]", s); |
---|
204 | s = "crk=3; Mu=14+2*r;"; |
---|
205 | write(l, "I_V#[1,2r-1]", s); |
---|
206 | s = "crk=3; Mu=15+2*r;"; |
---|
207 | write(l, "I_V#[1,2r]", s); |
---|
208 | s = "crk=3; Mu=0;"; |
---|
209 | write(l, "I_T[k,r,s]", s); |
---|
210 | close(l); |
---|
211 | } |
---|
212 | |
---|
213 | proc read_sing_dbm |
---|
214 | { |
---|
215 | link l="DBM: NFlist"; |
---|
216 | "A[k] = "+read(l, "A[k]"); |
---|
217 | "D[k] = "+read(l, "D[k]"); |
---|
218 | "E[6k] = "+read(l, "E[6k]"); |
---|
219 | "E[6k+1] = "+read(l, "E[6k+1]"); |
---|
220 | "E[6k+2] = "+read(l, "E[6k+2]"); |
---|
221 | "J[k,0] = "+read(l, "J[k,0]"); |
---|
222 | "J[k,r] = "+read(l, "J[k,r]"); |
---|
223 | "X[1,0] = "+read(l, "X[1,0]"); |
---|
224 | "X[1,r] = "+read(l, "X[1,r]"); |
---|
225 | "X[k,0] = "+read(l, "X[k,0]"); |
---|
226 | "X[k,r] = "+read(l, "X[k,r]"); |
---|
227 | "W[12k] = "+read(l, "W[12k]"); |
---|
228 | "W[12k+1] = "+read(l, "W[12k+1]"); |
---|
229 | "W[12k+5] = "+read(l, "W[12k+5]"); |
---|
230 | "W[12k+6] = "+read(l, "W[12k+6]"); |
---|
231 | "W[k,0] = "+read(l, "W[k,0]"); |
---|
232 | "W[k,r] = "+read(l, "W[k,r]"); |
---|
233 | "W#[k,2r-1] = "+read(l, "W#[k,2r-1]"); |
---|
234 | "W#[k,2r] = "+read(l, "W#[k,2r]"); |
---|
235 | "Y[1,r,s] = "+read(l, "Y[1,r,s]"); |
---|
236 | "Y[k,r,s] = "+read(l, "Y[k,r,s]"); |
---|
237 | "Z[1,0] = "+read(l, "Z[1,0]"); |
---|
238 | "Z[1,r] = "+read(l, "Z[1,r]"); |
---|
239 | "Z[k,r,s] = "+read(l, "Z[k,r,s]"); |
---|
240 | "Z[k,r,0] = "+read(l, "Z[k,r,0]"); |
---|
241 | "Z[k,12k+6r-1]= "+read(l, "Z[k,12k+6r-1]"); |
---|
242 | "Z[k,12k+6r] = "+read(l, "Z[k,12k+6r]"); |
---|
243 | "Z[k,12k+6r+1]= "+read(l, "Z[k,12k+6r+1]"); |
---|
244 | "Z[k,0] = "+read(l, "Z[k,0]"); |
---|
245 | "Z[k,r] = "+read(l, "Z[k,r]"); |
---|
246 | "Z[6k+5] = "+read(l, "Z[6k+5]"); |
---|
247 | "Z[6k+6] = "+read(l, "Z[6k+6]"); |
---|
248 | "Z[6k+7] = "+read(l, "Z[6k+7]"); |
---|
249 | "Q[k,0] = "+read(l, "Q[k,0]"); |
---|
250 | "Q[k,r] = "+read(l, "Q[k,r]"); |
---|
251 | "Q[6k+4] = "+read(l, "Q[6k+4]"); |
---|
252 | "Q[6k+5] = "+read(l, "Q[6k+5]"); |
---|
253 | "Q[6k+6] = "+read(l, "Q[6k+6]"); |
---|
254 | "S[12k-1] = "+read(l, "S[12k-1]"); |
---|
255 | "S[12k] = "+read(l, "S[12k]"); |
---|
256 | "S[k,0] = "+read(l, "S[k,0]"); |
---|
257 | "S[k,r] = "+read(l, "S[k,r]"); |
---|
258 | "S#[k,2r-1] = "+read(l, "S#[k,2r-1]"); |
---|
259 | "S#[k,2r] = "+read(l, "S#[k,2r]"); |
---|
260 | "S[12k+4] = "+read(l, "S[12k+4]"); |
---|
261 | "S[12k+5] = "+read(l, "S[12k+5]"); |
---|
262 | "U[12k] = "+read(l, "U[12k]"); |
---|
263 | "U[k,2r-1]= "+read(l, "U[k,2r-1]"); |
---|
264 | "U[k,2r] = "+read(l, "U[k,2r]"); |
---|
265 | "U[12k+4] = "+read(l, "U[12k+4]"); |
---|
266 | "V[1,0] = "+read(l, "V[1,0]"); |
---|
267 | "V[1,r] = "+read(l, "V[1,r]"); |
---|
268 | "V#[1,2r-1] = "+read(l, "V#[1,2r-1]"); |
---|
269 | "V#[1,2r] = "+read(l, "V#[1,2r]"); |
---|
270 | "T[k,r,s] = "+read(l, "T[k,r,s]"); |
---|
271 | close(l); |
---|
272 | } |
---|