read "itall.txt":


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


ver:=ourVertices(writelvl);

for i from 1 to numpsis do
	readFunc(T,12,cat(path,phatlvl,"phat",lvl,"psi.",i));
	fd:=fopen(cat(path,phatlvl,"phat",lvl,"psivals.",i),WRITE);
	fprintf(fd,"Values for psi %d generated from lvl %d on lvl %d\n",i,lvl,writelvl);
	for vertex in ver do
		fprintf(fd,"[");
		fprintf(fd,"%d",vertex[1]);
		for j from 2 to nops(vertex) do
			fprintf(fd,",%d",vertex[j]);
		end do;
		fprintf(fd,"]\t\t%20.20g\n",T[vertex,12]);
	end do;
	fclose(fd);
end do;
print("done printing psis");