read "itall.txt";

phatlvl:=8; 
lvl:=4; 
path:="../data/"; 
T:=gasket(phatlvl);
tol:=.000000000000000001;

vl:=ourVertices(lvl):
print(cat("nops vl at first= ",nops(vl)));
vl:=subsop(1=NULL,2=NULL,3=NULL,vl):

vp:=ourVertices(phatlvl):

print(cat("nops vl= ",nops(vl)));

for i from 1 to nops(vl) do
	clearFunc(T,phatlvl,13);
	for j from 1 to nops(vl) do
		readFunc(T,12,cat(path,phatlvl,"phat",lvl,"eigen.",j));
		print(cat("reading function: ",path,phatlvl,"phat",lvl,"eigen.",j));
		co:=T[vl[i],12];
		if (abs(co) > tol) then 	
			funcSumWithConstant(T,phatlvl,co,13,12,13)
		end if;
	end do;
	saveFunc(T,13,phatlvl,cat(path,phatlvl,"phat",lvl,"psi.",i));
	print(cat("saving function: ",path,phatlvl,"phat",lvl,"psi.",i));
end do;