read "itall.txt":

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

startlvl:=1:
endlvl:=5:
vertex:=1:
phatlvl:=8:
numpsis:=(3^(lvl+1)-3)/2:
path:=cat("../data2/"):
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 startlvl to endlvl do
		readFunc(T,12,cat(path,"/level",i,"/",phatlvl,"phat",i,"psi.",vertex));
		#ourPlot(T,phatlvl,12,cat(phatlvl,"phat",lvl,"psi.",i));
		ourPlotBlack(T,phatlvl,12,cat(phatlvl,"phat",i,"psi.",vertex));
	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;