read "itall.txt":

lvl:=2:
phatlvl:=4:
sl:=12:
rphatlvl:=8;
path:=cat("../data/level",lvl,"/");


verts:=subsop(1=NULL,2=NULL,3=NULL,ourVertices(lvl));
phatverts:=subsop(1=NULL,2=NULL,3=NULL,ourVertices(phatlvl));

badlist:=[];

for j from 1 to nops(phatverts) do
	if(nops(phatverts[j])=(lvl+2)) then
		badlist:=[op(badlist),j=NULL];
	end if;
end do;

phatverts:=subsop(op(badlist),phatverts);

locverts:=[];

for j from 1 to nops(verts) do
	locverts:=[op(locverts),[]];
	badlist:=[];
	for k from 1 to nops(phatverts) do
		if(sameCell(lvl+1,phatverts[k],verts[j])) then
			locverts[j]:=[op(locverts[j]),phatverts[k]];		
			badlist:=[op(badlist),k=NULL];
		end if;
	end do;
	phatverts:=subsop(op(badlist),phatverts);
end do;

for i from 1 to nops(verts) do
	readFunc(T,sl+i-1,cat(path,rphatlvl,"phat",lvl,"psi.",i));
end do;


lister:=rp([seq(i,i=1..nops(locverts[1]))],nops(locverts));
M:=matrix(nops(lister[1]),nops(lister[1]));


vals:=[];

for l in lister do
	for i from 1 to nops(l) do
		for j from 1 to nops(l) do
			M[i,j]:=T[ locverts[i][l[i]] , sl+j-1];
		end do;
	end do;
	vals:=[op(vals),[det(M),l]];
end do;


sorterfunc:=proc(a,b)
	if(a[1]<b[1]) then return(true);
	else return(false);
	end if;
end;

vals:=sort(vals,sorterfunc);

fd:=fopen(cat(path,phatlvl,"phat",lvl,"dets.txt"),WRITE);
fprintf(fd,"%d phat %d determinants\n",phatlvl,lvl);
for val in vals do
	fprintf(fd,"%20.20g\n",val[1]);
	l:=val[2];
	fprintf(fd,"\t");
	for i from 1 to nops(l) do
		fprintf(fd,"|%a-%a|=%d",locverts[i][l[i]],verts[i],verDist(phatlvl,locverts[i][l[i]],verts[i]));
		if i=nops(l) then
			fprintf(fd,"\n");
		else fprintf(fd,"\t");
		end if;
	end do;
	fprintf(fd,"\t");
	for i from 1 to (nops(l)-1) do
		for j from (i+1) to nops(l) do
			fprintf(fd,"|%a-%a|=%d",locverts[i][l[i]],locverts[j][l[j]],verDist(phatlvl,locverts[i][l[i]],locverts[j][l[j]]));
			if not ( (i=nops(l)) and (j=nops(l)) ) then
				fprintf(fd,"\t");
			end if;
		end do;
	end do;
	fprintf(fd,"\n");
end do;
fclose(fd);