Digits:=20;
read "everything.txt":
read "func_file.txt":
read "funcarith.txt":
read "autoset.txt":

T:=gasket(9):
eigenlvl:=4;
phatlvl:=8;
numEigens:=(3^(eigenlvl+1)-3)/2;
startEigen:=12;
path:="../data/";

for i from 1 to numEigens do
	readFunc(T,i+startEigen-1,cat(path,phatlvl,"phat",eigenlvl,"eigen.",i));
end do;
#for i from 1 to numEigens do
#	normalizeFunc(T,phatlvl,i+startEigen-1);
#end do;
ver:=ourVertices(phatlvl);
smallver:=ourVertices(eigenlvl);
#ver:=subsop(1=NULL,2=NULL,3=NULL,ver);
smallver:=subsop(1=NULL,2=NULL,3=NULL,smallver);
funcNum:=startEigen+numEigens+1;
for y from 1 to numEigens do
	clearFunc(T,phatlvl,funcNum);
	for x in ver do
		for j from 1 to numEigens do

T[x,funcNum]:=T[x,funcNum]+T[smallver[y],j+startEigen-1]*T[x,j+startEigen-1];
		end do;
	end do;
	print(cat(path,phatlvl,"phat",eigenlvl,"psi.",y));
	saveFunc(T,funcNum,phatlvl,cat(path,phatlvl,"phat",eigenlvl,"psi.",y));
end do;