read "itall.txt":

PRINTEIGENS:=0:
PRINTPSIS:=0:
PRINTDIRS:=0:

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



if (PRINTEIGENS=1) then 
	for i from 1 to numpsis do
	        readFunc(T,12,cat(path,phatlvl,"phat",lvl,"eigen.",i)):  
		displayToPs(cat("/tmp/richard/",phatlvl,"phat",lvl,"eigen",i,".ps"));
	        #ourPlot(T,phatlvl,12,cat(phatlvl,"phat",lvl,"eigen.",i));
	        ourPlot(T,phatlvl,12,"");
	end do;
	print("done printing eigens");
end if;



#if(PRINTPSIS=1) then 
	for i from 1 to numpsis do
		readFunc(T,12,cat(path,phatlvl,"phat",lvl,"psi.",i));
		#ourPlot(T,phatlvl,12,cat(phatlvl,"phat",lvl,"psi.",i));
		ourPlotBlack(T,phatlvl,12,cat(phatlvl,"phat",lvl,"psi.",i));
	end do;
	print("done printing psis");
#end if;

if(PRINTDIRS=1) then
        for i from 1 to numpsis do
                readFunc(T,12,cat(path,phatlvl,"phat",lvl,"dirichlet.",i));
                ourPlot(T,phatlvl,12,cat(phatlvl,"phat",lvl,"dirichlet.",i));
        end do;
        print("done printing dirichlet kernels");
end if;