source: git/Tst/Short/symodstd.tst @ 2f6c055

spielwiese
Last change on this file since 2f6c055 was 2f6c055, checked in by Christian Eder, 12 years ago
added Tst/short/Short/symodstd.* it is rest of commit 8979f91: modstd.lib: Example in modpStd has been changed. symodstd.lib: Library for computing Groebner basis of symmetric ideals. Conflicts: Singular/LIB/modstd.lib Singular/LIB/symodstd.lib
  • Property mode set to 100644
File size: 1.9 KB
Line 
1LIB "tst.lib";
2tst_init();
3
4LIB "symodstd.lib";
5
6example genSymId;
7
8example isSymmetric;
9
10example primRoot;
11
12example eigenvalues;
13
14intvec pi, tau, sigma;
15
16ring R1 = 0, x(1..4), dp;
17ideal I = cyclic(4);
18pi = 2,3,4,1;
19symmStd(I,pi);
20
21ring R2 = 31, (x,y,z), dp;
22ideal J;
23J[1] = xy-y2+xz;
24J[2] = xy+yz-z2;
25J[3] = -x2+xz+yz;
26tau = 3,1,2;
27symmStd(J,tau);
28
29ring R3 = 30817, x(1..8), dp;
30ideal I = cyclic(8);
31sigma = 2,3,4,5,6,7,8,1;
32symmStd(I,sigma);
33
34ring R4 = 0, (x,y,z), dp;
35ideal I;
36I[1] = -2xyz4+xz5+xz;
37I[2] = -2xyz4+yz5+yz;
38pi = 2,1,3;
39syModStd(I,pi);
40
41ring R5 = 0, x(1..4), dp;
42ideal I = cyclic(4);
43tau = 2,3,4,1;
44ideal sJ1 = syModStd(I,tau,3);
45ideal sJ2 = syModStd(I,tau,1,1);
46size(reduce(sJ1,sJ2));
47size(reduce(sJ2,sJ1));
48
49ring R6 = 0, (a(1..3),b(1..3),c(1..2)), dp;
50ideal I;
51I[1]=4*a(1)^5*a(2)+2*b(1)*b(2)*b(3)*c(2);
52I[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);
53I[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);
54I[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);
55I[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);
56I[6]=20*a(1)+4*a(2)+2*b(1)-2*b(2)-3*b(3)-5;
57I[7]=4*b(1)^5*b(2)+2*a(1)*a(2)*a(3)*c(2);
58I[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);
59I[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);
60I[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);
61I[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);
62I[12]=2*a(1)-2*a(2)-3*a(3)+20*b(1)+4*b(2)-5;
63sigma = 4,5,6,1,2,3,7,8;
64ideal sI1 = syModStd(I,sigma);
65ideal sI2 = modStd(I);
66size(reduce(sI1,sI2));
67size(reduce(sI2,sI1));
68
69tst_status(1);$
Note: See TracBrowser for help on using the repository browser.