read "itall.txt":

lvl:=1:
phatlvl:=8:
numpsis:=(3^(lvl+1)-3)/2:
path:=cat("../data/","level",lvl,"/"):
T:=gasket(phatlvl):


clearFunc(T,phatlvl, 13):
for i from 1 to numpsis do
	readFunc(T,12,cat(path,phatlvl,"phat",lvl,"psi.",i));
	funcAbs(T,phatlvl,12,12);
	funcSumWithConstant(T,phatlvl,1,12,13,13);
end do;
saveFunc(T,13,phatlvl,cat(path,phatlvl,"phatsumabspsi.txt"));