source: git/doc/orbitsG25.gp @ 39d022e

spielwiese
Last change on this file since 39d022e was 6cfc68, checked in by Janko Boehm <boehm@…>, 7 years ago
Added GAP script to compute orbit decomposition
  • Property mode set to 100644
File size: 2.5 KB
Line 
1G := Group(
2 (2, 3)(5, 6)(9, 10), (1, 5, 9, 10, 3)(2, 7, 8, 4, 6)
3);;
4Size(G);
5dimQ:=5;
6
7minidx:=SmallestMovedPoint(G);
8maxidx:=LargestMovedPoint(G);
9XZorbitsRepresentatives:=[];;
10for k in [dimQ..(maxidx-minidx+1)] do
11  Print("Considering faces of cardinality ",k,"\n");
12  XZ := Combinations([minidx..maxidx],k);;
13  Bahnen := OrbitsDomain(G,XZ,OnSets);;
14  for i in [1..Size(Bahnen)] do
15    Append(XZorbitsRepresentatives,[Bahnen[i][1]]);
16  od;
17od;
18
19PrintTo("simplexOrbitRepresentativesG25.sing","list simplexOrbitRepresentatives = ");
20for k in [1..Size(XZorbitsRepresentatives)-1] do
21  s:=String(XZorbitsRepresentatives[k]);
22  s:=s{[2..Size(s)-1]};
23  AppendTo ("simplexOrbitRepresentativesG25.sing","intvec(",s,"),\n");
24od;
25s:=String(XZorbitsRepresentatives[Size(XZorbitsRepresentatives)]);
26s:=s{[2..Size(s)-1]};
27AppendTo ("simplexOrbitRepresentativesG25.sing","intvec(",s,");\n");
28
29
30PrintTo("simplexSymmetryGroupG25.sing","list simplexSymmetryGroup = ");
31elementsG:=Elements(G);
32for i in [1..Size(elementsG)-1] do
33   sigma:=elementsG[i];
34   l:=ListPerm(sigma,maxidx);
35   l:=l{[minidx..maxidx]};
36   s:=String(l);
37   s:=s{[2..Size(s)-1]};
38   AppendTo ("simplexSymmetryGroupG25.sing","permutationFromIntvec(intvec(",s,")),\n");
39od;
40sigma:=elementsG[Size(elementsG)];
41l:=ListPerm(sigma,maxidx);
42l:=l{[minidx..maxidx]};
43s:=String(l);
44s:=s{[2..Size(s)-1]};
45AppendTo ("simplexSymmetryGroupG25.sing","permutationFromIntvec(intvec(",s,"));\n");
46
47
48PrintTo("elementsInTermsOfGeneratorsG25.sing","list generatorsG = ");
49L:=GeneratorsOfGroup(G);
50for i in [1..Size(L)-1] do
51   sigma:=L[i];
52   l:=ListPerm(sigma,maxidx);
53   l:=l{[minidx..maxidx]};
54   s:=String(l);
55   s:=s{[2..Size(s)-1]};
56   AppendTo ("elementsInTermsOfGeneratorsG25.sing","permutationFromIntvec(intvec(",s,")),\n");
57od;
58   sigma:=L[Size(L)];
59   l:=ListPerm(sigma,maxidx);
60   l:=l{[minidx..maxidx]};
61   s:=String(l);
62   s:=s{[2..Size(s)-1]};
63   AppendTo ("elementsInTermsOfGeneratorsG25.sing","permutationFromIntvec(intvec(",s,"));\n");
64
65AppendTo("elementsInTermsOfGeneratorsG25.sing","list elementsInTermsOfGenerators = ");
66hom:=EpimorphismFromFreeGroup(G);
67   AppendTo ("elementsInTermsOfGeneratorsG25.sing","\"\",\n");
68for i in [2..Size(elementsG)-1] do
69   sigma:=elementsG[i];
70   l:=PreImagesRepresentative(hom,sigma);
71   s:=String(l);
72   AppendTo ("elementsInTermsOfGeneratorsG25.sing","\"",s,"\",\n");
73od;
74sigma:=elementsG[Size(elementsG)];
75l:=PreImagesRepresentative(hom,sigma);
76s:=String(l);
77AppendTo ("elementsInTermsOfGeneratorsG25.sing","\"",s,"\";\n");
78
79
Note: See TracBrowser for help on using the repository browser.