1 | LIB "tst.lib"; |
---|
2 | tst_init(); |
---|
3 | |
---|
4 | LIB "mondromy.lib"; |
---|
5 | |
---|
6 | list unimodal= |
---|
7 | "P[8]","x,y,z","x3+y3+z3+xyz", |
---|
8 | "X[9]","x,y","x4+y4+x2y2", |
---|
9 | "J[10]","x,y","x3+y6+x2y2", |
---|
10 | "T[3,4,5]","x,y,z","x3+y4+z5+xyz", |
---|
11 | "T[3,4,6]","x,y,z","x3+y4+z6+xyz", |
---|
12 | "T[4,5,6]","x,y,z","x4+y5+z6+xyz", |
---|
13 | "E[12]","x,y","x3+y7+xy5", |
---|
14 | "E[13]","x,y","x3+xy5+y8", |
---|
15 | "E[14]","x,y","x3+y8+xy6", |
---|
16 | "Z[11]","x,y","x3y+y5+xy4", |
---|
17 | "Z[12]","x,y","x3y+xy4+x2y3", |
---|
18 | "Z[13]","x,y","x3y+y6+xy5", |
---|
19 | "W[12]","x,y","x4+y5+x2y3", |
---|
20 | "W[13]","x,y","x4+xy4+y6", |
---|
21 | "Q[10]","x,y,z","x3+y4+yz2+xy3", |
---|
22 | "Q[11]","x,y,z","x3+y2z+xz3+z5", |
---|
23 | "Q[12]","x,y,z","x3+y5+yz2+xy4", |
---|
24 | "S[11]","x,y,z","x4+y2z+xz2+x3z", |
---|
25 | "S[12]","x,y,z","x2y+y2z+xz3+z5", |
---|
26 | "U[12]","x,y,z","x3+y3+z4+xyz2"; |
---|
27 | |
---|
28 | list bimodal= |
---|
29 | //"J[3,0]","x,y","x3+x2y3+y9+xy7", |
---|
30 | "J[3,1]","x,y","x3+x2y3+y10", |
---|
31 | "J[3,2]","x,y","x3+x2y3+y11", |
---|
32 | //"Z[1,0]","x,y","x3y+x2y3+xy6+y7", |
---|
33 | "Z[1,1]","x,y","x3y+x2y3+y8", |
---|
34 | "Z[1,2]","x,y","x3y+x2y3+y9", |
---|
35 | "W[1,0]","x,y","x4+x2y3+y6", |
---|
36 | "W[1,1]","x,y","x4+x2y3+y7", |
---|
37 | "W[1,2]","x,y","x4+x2y3+y8", |
---|
38 | "W#[1,1]","x,y","(x2+y3)^2+xy5", |
---|
39 | "W#[1,2]","x,y","(x2+y3)^2+x2y4", |
---|
40 | "W#[1,3]","x,y","(x2+y3)^2+xy6", |
---|
41 | //"W#[1,4]","x,y","(x2+y3)^2+x2y5", |
---|
42 | "Q[2,0]","x,y,z","x3+yz2+x2y2+xy4", |
---|
43 | "Q[2,1]","x,y,z","x3+yz2+x2y2+y7", |
---|
44 | "Q[2,2]","x,y,z","x3+yz2+x2y2+y8", |
---|
45 | "S[1,0]","x,y,z","x2z+yz2+y5+zy3", |
---|
46 | "S[1,1]","x,y,z","x2z+yz2+x2y2+y6", |
---|
47 | "S[1,2]","x,y,z","x2z+yz2+x2y2+y7", |
---|
48 | "S#[1,1]","x,y,z","x2z+yz2+zy3+xy4", |
---|
49 | "S#[1,2]","x,y,z","x2z+yz2+zy3+x2y3", |
---|
50 | "S#[1,3]","x,y,z","x2z+yz2+zy3+xy5", |
---|
51 | //"S#[1,4]","x,y,z","x2z+yz2+zy3+x2y4", |
---|
52 | "U[1,0]","x,y,z","x3+xz2+xy3+y3z", |
---|
53 | "U[1,1]","x,y,z","x3+xz2+xy3+y2z2", |
---|
54 | "U[1,2]","x,y,z","x3+xz2+xy3+y4z", |
---|
55 | "U[1,3]","x,y,z","x3+xz2+xy3+y3z2", |
---|
56 | "U[1,4]","x,y,z","x3+xz2+xy3+y5z", |
---|
57 | "E[18]","x,y","x3+y10+xy7", |
---|
58 | "E[19]","x,y","x3+xy7+y11", |
---|
59 | "E[20]","x,y","x3+y11+xy8", |
---|
60 | "Z[17]","x,y","x3y+y8+xy6", |
---|
61 | "Z[18]","x,y","x3y+xy6+y9", |
---|
62 | "Z[19]","x,y","x3y+y9+xy7", |
---|
63 | "W[17]","x,y","x4+xy5+y7", |
---|
64 | "W[18]","x,y","x4+y7+x2y4", |
---|
65 | "Q[16]","x,y,z","x3+yz2+y7+xy5", |
---|
66 | //"Q[17]","x,y,z","x2z+yz2+xy5+y8", |
---|
67 | "Q[18]","x,y,z","x3+yz2+y8+xy6", |
---|
68 | "S[16]","x,y,z","x2z+yz2+xy4+y6", |
---|
69 | "S[17]","x,y,z","x2z+yz2+y6+zy4", |
---|
70 | "U[16]","x,y,z","x3+xz2+y5+x2y2"; |
---|
71 | |
---|
72 | proc tst_monodromy(list data) |
---|
73 | { |
---|
74 | int i; |
---|
75 | string s,typ; |
---|
76 | for(i=1;i<=size(data);i=i+3) |
---|
77 | { |
---|
78 | s="typ=\""+data[i]+"\";"; |
---|
79 | execute(s); |
---|
80 | s="ring R=0,("+data[i+1]+"),ds;"; |
---|
81 | execute(s); |
---|
82 | export(R); |
---|
83 | s="poly f="+data[i+2]+";"; |
---|
84 | execute(s); |
---|
85 | typ; |
---|
86 | jordan(monodromyB(f)); |
---|
87 | tst_status(); |
---|
88 | kill R; |
---|
89 | } |
---|
90 | } |
---|
91 | |
---|
92 | tst_monodromy(unimodal); |
---|
93 | |
---|
94 | tst_status(1); $ |
---|