1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
|
######################################################
######################################################
InstallGlobalFunction(PolymakeFaceLattice,
function(arg)
local polygon, toggle, F, FF, FA, D, FL, i, j, cnt, lst,bool;
#WARNING: This fixed is used in HAPcryst. Don't fiddle with toggle
polygon:=arg[1];
if Length(arg)=2 then toggle:=false;
else toggle:=true;
fi;
bool:=false;
FF:=Polymake(polygon,"FACES");
FA:=Polymake(polygon,"ADJACENCY");
D:=Polymake(polygon,"F_VECTOR");
if not Size(Filtered(FF,x->Size(x)=1))=D[1] then D:=Reversed(D); fi;
if not FF[1]=[] then FF:=Reversed(FF); FA:=Reversed(FA); bool:=true; fi;
if toggle then F:=FF; else F:=FA; fi;
FL:=[];;
FL[1]:=[F[1]];
cnt:=1;
for i in [1..Length(D)] do
FL[i+1]:=[];
lst:=[cnt+1..cnt+D[i]];
if bool then lst:=Reversed(lst); fi;
for j in lst do
Add(FL[i+1],F[j]);
od;
cnt:=cnt+D[i];
od;
if FL[1]=[[]] then FL:=Reversed(FL); else Add(FL,[[]]); fi;
if not toggle then Remove(FL,1); fi;
return FL;
end);
############################################################
############################################################
|
