read "itall.txt":

setupdisp:=proc(T,phatlvl,path,lvl,i)
	local start_points,plus_points,minus_points,pl,pnt;
	readFunc(T,12,cat(path,phatlvl,"phat",lvl,"psi.",i));
	make_points(T,phatlvl):
	start_points:=SG_plot(T,phatlvl,12):
	plus_points:=[]:
	minus_points:=[]:
	pl:=[]:
	for pnt in start_points do
		if pnt[3]>=0 then
			pnt[3]:=0:
			plus_points:=[op(plus_points),pnt]:
		else
			pnt[3]:=0;
			minus_points:=[op(minus_points),pnt]:
		end if:
	end do:
	pl:=[pointplot3d(plus_points,color=black,symbol=POINT,orientation=[0,0]),pointplot3d(minus_points,color=grey,symbol=POINT,orientation=[0,0])]:
	RETURN(pl);
end:


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


make_points(T,phatlvl):
ver:=subsop(1=NULL,2=NULL,3=NULL,vertices(lvl)):
for i from 1 to numrand do
	k:=rand() mod numpsis:
	print(k);
	print(ver[k]);
	display3d(setupdisp(T,phatlvl,path,lvl,k));
end do;
print("done printing psis");