read "itall.txt":

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

fd:=fopen(cat(path,phatlvl,"phat",lvl,"dirinnner.txt"),WRITE);
fprintf(fd,"l1, l2, and integral data for the %d phat %d dirichelt\n",phatlvl,lvl);

for i from 1 to numpsis do
	readFunc(T,12,cat(path,phatlvl,"phat",lvl,"dirichlet.",i));
	ans:=discInnerProduct(T,phatlvl,12,12)/nops(ourVertices(phatlvl));	
	fprintf(fd,"%d: %20.20g\n",i,ans);
end do;
fclose(fd);