read "itall.txt":


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


ver:=ourVertices(lvl);
ver:=subsop(1=NULL, 2=NULL, 3=NULL, ver);

fd:=fopen(cat("../data/",lvl,"vertexlist.txt"),WRITE);

k:=1;

for vertex in ver do
	fprintf(fd,"%d\t",k);
	fprintf(fd,"[");
	fprintf(fd,"%d",vertex[1]);
	for j from 2 to nops(vertex) do
		fprintf(fd,",%d",vertex[j]);
	end do;
	fprintf(fd,"]\n");
	k:=k+1;
end do;

fclose(fd);