read "itall.txt";

lvl:=4;
phatlvl:=8;

T:=gasket(phatlvl);

inpath:=cat("../data/level",lvl,"/");
outpath:="../data/";

verts:=ourVertices(lvl):
verts:=subsop(1=NULL, 2=NULL, 3=NULL, verts):

for k from 1 to nops(verts) do
	readFunc(T,12,cat(inpath,phatlvl,"phat",lvl,"psi.",k));
	clearFunc(T,phatlvl,13);
	copyRotateScaleChunk(T,lvl,phatlvl,12,13,verts[k]);
	saveFunc(T,13,(phatlvl+2-lvl),cat(outpath,phatlvl,"phat",lvl,"spike.",k));
end do;

readFunc(T,12,cat(outpath,phatlvl,"phat",lvl,"spike.",1));

for k from 1 to nops(verts) do
	readFunc(T,13,cat(outpath,phatlvl,"phat",lvl,"spike.",k));
	clearFunc(T,(phatlvl+2-lvl),14);
	funcDifference(T,(phatlvl+2-lvl),12,13,14);
        saveFunc(T,14,(phatlvl+2-lvl),cat(outpath,phatlvl,"phat",lvl,"spikediffwith1.",k));
end do;