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
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
|
#C 2008 Graham Ellis
#####################################################################
InstallMethod(HAPPrintTo,
"Method for writing HAP resolutions to a file",
[IsString,IsHapResolution],
function(file,R)
local n,k,g,AppendTo,PrintTo;
AppendTo:=HAP_AppendTo;
PrintTo:=HAP_PrintTo;
#Create a gap function
PrintTo(file,"HAPTEMPORARYFUNCTION:=function() local TYPE,RANKS, HIGHEST_DEGREE,BOUNDARIES,CONTRACTING_HOMOTOPY,ELEMENTS,SPECIAL_PROPERTIES;\n\n");
#Export type of HAP object
AppendTo(file,"TYPE:=\nHapResolution;\n\n");
#Export length of resolution
AppendTo(file,"HIGHEST_DEGREE:= \n");
AppendTo(file,Length(R),";\n\n");
#Export ranks of the resolution in each degree.
AppendTo(file,"RANKS:=\n");
AppendTo(file,List([0..Length(R)],i->Dimension(R)(i)),";\n\n");
#Export boundary words
AppendTo(file,"#BOUNDARIES[n][k] is the boundary of the kth free generator in degree n\n");
AppendTo(file,"BOUNDARIES:=","\n[\n");
for n in [1..Length(R)] do
AppendTo(file,"#Degree ", n,"\n","[\n\n" );
for k in [1..R!.dimension(n)] do
AppendTo(file,R!.boundary(n,k));
if k<R!.dimension(n) then AppendTo(file,",\n\n");
else AppendTo(file,"\n\n");
fi;
od;
if n<Length(R) then
AppendTo(file,"],\n\n");
else
AppendTo(file,"]\n];\n\n");
fi;
od;
#Export the contracting homotopy
AppendTo(file,"#CONTRATCTING_HOMOTOPY[n][k][g] is the image of the free\n");
AppendTo(file,"#abelian group generator g(e_k)^n in degree n.\n");
if R!.homotopy=fail then
AppendTo(file, "CONTRACTING_HOMOTOPY:=\nfail;\n\n");
else
AppendTo(file, "CONTRACTING_HOMOTOPY:=\n[\n\n");
for n in [0..Length(R)-2] do
AppendTo(file,"#Degree ", n,"\n","[\n\n" );
for k in [1..R!.dimension(n)] do
AppendTo(file,"#Free Generator ",k," in degree ",n," \n[\n\n");
for g in [1..Length(R!.elts)] do
AppendTo(file,R!.homotopy(n,[k,g]));
AppendTo(file,",\n\n");
od;
AppendTo(file,"\n],\n\n");
od;
if n<Length(R)-2 then
AppendTo(file,"],\n\n");
else
AppendTo(file,"]\n];\n\n");
fi;
od;
fi;
#Export group elements
AppendTo(file, "ELEMENTS:=\n");
AppendTo(file, R!.elts,";\n\n");
#Export special properties of the resolution
AppendTo(file,"SPECIAL_PROPERTIES:=\n", R!.properties,";\n\n");
AppendTo(file, "return rec( type:=TYPE, highest_degree:=HIGHEST_DEGREE, ranks:=RANKS, boundaries:=BOUNDARIES,contracting_homotopy:=CONTRACTING_HOMOTOPY, elements:=ELEMENTS,special_properties:=SPECIAL_PROPERTIES);\n\n");
AppendTo(file,"end;\n");
end);
#####################################################################
#####################################################################
InstallOtherMethod(HAPPrintTo,
"Method for writing a HAP regular CW-complex to a file",
[IsString, IsHapRegularCWComplex],
function(file,X)
local Dims;
#Create a gap function
PrintTo(file,"HAPTEMPORARYFUNCTION:=Objectify(HapRegularCWComplex,rec()); \n\n");
AppendTo(file,"HAPTEMPORARYFUNCTION!.boundaries:=\n");
AppendTo(file,X!.boundaries);
AppendTo(file,";\n\n");
AppendTo(file,"HAPTEMPORARYFUNCTION!.coboundaries:=\n");
AppendTo(file,X!.coboundaries);
AppendTo(file,";\n\n");
AppendTo(file,"HAPTEMPORARYFUNCTION!.vectorField:=\n");
AppendTo(file,X!.vectorField);
AppendTo(file,";\n\n");
AppendTo(file,"HAPTEMPORARYFUNCTION!.inverseVectorField:=\n");
AppendTo(file,X!.inverseVectorField);
AppendTo(file,";\n\n");
AppendTo(file,"HAPTEMPORARYFUNCTION!.criticalCells:=\n");
AppendTo(file,X!.criticalCells);
AppendTo(file,";\n\n");
AppendTo(file,"HAPTEMPORARYFUNCTION!.orientation:=\n");
AppendTo(file,X!.orientation);
AppendTo(file,";\n\n");
AppendTo(file,"HAPTEMPORARYFUNCTION!.properties:=\n");
AppendTo(file,X!.properties);
AppendTo(file,";\n\n");
Dims:=List([0..Dimension(X)],n->X!.nrCells(n));
AppendTo(file,"DIMS:=\n");
AppendTo(file,Dims);
AppendTo(file,";\n\n");
AppendTo(file,"HAPTEMPORARYFUNCTION!.nrCells:=");
AppendTo(file,"function(n) if IsBound(DIMS[n+1]) then return DIMS[n+1]; fi; return 0; end;");
AppendTo(file,";\n\n");
end);
#####################################################################
|
