read "itall.txt":

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


for i from 1 to numpsis do
	readFunc(T,startfunc-1+i,cat(path,phatlvl,"phat",lvl,"psi.",i));
end do;

phatver:=subsop(1=NULL,2=NULL,3=NULL,ourVertices(phatlvl)):
ver:=ourVertices(lvl):
ver:=subsop(1=NULL,2=NULL,3=NULL,ver):
M:=matrix(nops(ver),nops(ver)):
for i from 1 to nops(ver) do
	for j from 1 to nops(ver) do
		M[i,j]:=T[ver[i],startfunc-1+j];
	end do;
end do;

for vertex in phatver do
	for j from 1 to nops(ver) do
		M[1,j]:=T[vertex,startfunc-1+j];
	end do;
	fd:=fopen(cat(path,phatlvl,"phat",lvl,"dets.txt"),APPEND);
	fprintf(fd,"1->%a = %20.20g\n",vertex,det(M));
	fclose(fd);
end do:
#print("det M=");
#print(det(M));