read "itall.txt":


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


fd:=fopen(cat(path,phatlvl,"phat",lvl,"psiminmax.txt"),WRITE);

for i from 1 to numpsis do
	readFunc(T,12,cat(path,phatlvl,"phat",lvl,"psi.",i));
	ma:=findFuncMax2(T,phatlvl,12);
	mi:=findFuncMin2(T,phatlvl,12);
	fprintf(fd,"function %d:\nmax - [%d",i,ma[2][1]);
	for j from 2 to nops(ma[2]) do
		fprintf(fd,",%d",ma[2][j]);
	end do;
	fprintf(fd,"] - %20.20g\n",ma[1]);

        fprintf(fd,"min - [%d",mi[2][1]);
        for j from 2 to nops(mi[2]) do
                fprintf(fd,",%d",mi[2][j]);
        end do;
        fprintf(fd,"] - %20.20g\n",mi[1]);


	#fprintf(fd,"%d has max\t%20.20g and min\t%20.20g\n",i,ma,mi);
end do;
fclose(fd);
print("done printing psis");