1 | LIB "tst.lib"; |
---|
2 | tst_init(); |
---|
3 | |
---|
4 | LIB "symodstd.lib"; |
---|
5 | |
---|
6 | example genSymId; |
---|
7 | |
---|
8 | example isSymmetric; |
---|
9 | |
---|
10 | example primRoot; |
---|
11 | |
---|
12 | example eigenvalues; |
---|
13 | |
---|
14 | intvec pi, tau, sigma; |
---|
15 | |
---|
16 | ring R1 = 0, x(1..4), dp; |
---|
17 | ideal I = cyclic(4); |
---|
18 | pi = 2,3,4,1; |
---|
19 | symmStd(I,pi); |
---|
20 | |
---|
21 | ring R2 = 31, (x,y,z), dp; |
---|
22 | ideal J; |
---|
23 | J[1] = xy-y2+xz; |
---|
24 | J[2] = xy+yz-z2; |
---|
25 | J[3] = -x2+xz+yz; |
---|
26 | tau = 3,1,2; |
---|
27 | symmStd(J,tau); |
---|
28 | |
---|
29 | ring R3 = 30817, x(1..8), dp; |
---|
30 | ideal I = cyclic(8); |
---|
31 | sigma = 2,3,4,5,6,7,8,1; |
---|
32 | symmStd(I,sigma); |
---|
33 | |
---|
34 | ring R4 = 0, (x,y,z), dp; |
---|
35 | ideal I; |
---|
36 | I[1] = -2xyz4+xz5+xz; |
---|
37 | I[2] = -2xyz4+yz5+yz; |
---|
38 | pi = 2,1,3; |
---|
39 | syModStd(I,pi); |
---|
40 | |
---|
41 | ring R5 = 0, x(1..4), dp; |
---|
42 | ideal I = cyclic(4); |
---|
43 | tau = 2,3,4,1; |
---|
44 | ideal sJ1 = syModStd(I,tau,3); |
---|
45 | ideal sJ2 = syModStd(I,tau,1,1); |
---|
46 | size(reduce(sJ1,sJ2)); |
---|
47 | size(reduce(sJ2,sJ1)); |
---|
48 | |
---|
49 | ring R6 = 0, (a(1..3),b(1..3),c(1..2)), dp; |
---|
50 | ideal I; |
---|
51 | I[1]=4*a(1)^5*a(2)+2*b(1)*b(2)*b(3)*c(2); |
---|
52 | I[2]=4*a(1)^5+20*a(1)^4*a(2)+3*b(1)*b(2)*b(3)*c(1)+b(1)*b(2)*c(2)-4*b(2)*b(3)*c(2); |
---|
53 | I[3]=20*a(1)^4+40*a(1)^3*a(2)+4*b(1)*b(2)*b(3)+2*b(1)*b(2)*c(1)+b(1)*b(3)*c(1)-3*b(2)*b(3)*c(1)-b(1)*c(2)-5*b(2)*c(2)-6*b(3)*c(2); |
---|
54 | I[4]=40*a(1)^3+40*a(1)^2*a(2)+5*b(1)*b(2)*b(3)+3*b(1)*b(2)+2*b(1)*b(3)-2*b(2)*b(3)-4*b(2)*c(1)-5*b(3)*c(1)-7*c(2); |
---|
55 | I[5]=40*a(1)^2+20*a(1)*a(2)+4*b(1)*b(2)+3*b(1)*b(3)-b(2)*b(3)+b(1)-3*b(2)-4*b(3)-6*c(1); |
---|
56 | I[6]=20*a(1)+4*a(2)+2*b(1)-2*b(2)-3*b(3)-5; |
---|
57 | I[7]=4*b(1)^5*b(2)+2*a(1)*a(2)*a(3)*c(2); |
---|
58 | I[8]=4*b(1)^5+20*b(1)^4*b(2)+3*a(1)*a(2)*a(3)*c(1)+a(1)*a(2)*c(2)-4*a(2)*a(3)*c(2); |
---|
59 | I[9]=20*b(1)^4+40*b(1)^3*b(2)+4*a(1)*a(2)*a(3)+2*a(1)*a(2)*c(1)+a(1)*a(3)*c(1)-3*a(2)*a(3)*c(1)-a(1)*c(2)-5*a(2)*c(2)-6*a(3)*c(2); |
---|
60 | I[10]=5*a(1)*a(2)*a(3)+40*b(1)^3+40*b(1)^2*b(2)+3*a(1)*a(2)+2*a(1)*a(3)-2*a(2)*a(3)-4*a(2)*c(1)-5*a(3)*c(1)-7*c(2); |
---|
61 | I[11]=4*a(1)*a(2)+3*a(1)*a(3)-a(2)*a(3)+40*b(1)^2+20*b(1)*b(2)+a(1)-3*a(2)-4*a(3)-6*c(1); |
---|
62 | I[12]=2*a(1)-2*a(2)-3*a(3)+20*b(1)+4*b(2)-5; |
---|
63 | sigma = 4,5,6,1,2,3,7,8; |
---|
64 | ideal sI1 = syModStd(I,sigma); |
---|
65 | ideal sI2 = modStd(I); |
---|
66 | size(reduce(sI1,sI2)); |
---|
67 | size(reduce(sI2,sI1)); |
---|
68 | |
---|
69 | tst_status(1);$ |
---|