#this file isn't used too much anymore
#it was for the old method of generating the psi functions, were we
#just inverted the matrix of phi functions.
#in this method we had to type a lot of the functions in by hand
#we did gthat in this file.


#this is [seq(-1,1..num)], we didn't knwo seq yet...
makeEpsilons:=proc(num)
	local i,l;
	l:=[];
	for i from 1 to num do
		l:=add_one(l,-1);
	end do;
	RETURN(l);
end;

#this copies inf to outf at level lvl

copyfunc:=proc(T,lvl,inf,outf)
	local vm,ver;
	vm=vertices(lvl);
	for ver in vm do
		T[ver,outf]:=T[ver,inf];
	end do;
end;

setfuncsNewMethod:=proc(T,lambdas,startlevels)
	local i,j;
	for i from 0 to 2 do
		for j from 12 to 13 do
			T[[i],j]:=0;
		end do;
	end do;

	lambdas[12]:=evalf((5-sqrt(5))/2,20);
	lambdas[13]:=evalf((5-sqrt(5))/2,20);

	startlevels[12]:=1;
	startlevels[13]:=1;
	
	i:=12;
	T[[0,1],i]:=0;
	T[[0,2],i]:=evalf(1/sqrt(2),20);
	T[[1,2],i]:=evalf(-1/sqrt(2),20);
	
	i:=13;
	T[[0,1],i]:=evalf(-2/sqrt(6),20);
	T[[0,2],i]:=evalf(1/sqrt(6),20);
	T[[1,2],i]:=evalf(1/sqrt(6),20);

end;

setfuncs:=proc(T,lambdas,startlevels)
	local i,j,l;

	lambdas[12]:=(5-sqrt(17))/2;
	lambdas[13]:=(5+sqrt(17))/2;
	lambdas[14]:=(5-sqrt(5))/2;
	lambdas[15]:=(5+sqrt(5))/2;
	lambdas[16]:=(5-sqrt(5))/2;
	lambdas[17]:=(5+sqrt(5))/2;#used to be sqrt(17)
	lambdas[18]:=3;
	lambdas[19]:=3;
	lambdas[20]:=3;
	lambdas[21]:=(5-sqrt(5))/2;
	lambdas[22]:=(5-sqrt(5))/2;
	lambdas[23]:=(5-sqrt(5))/2;
	
	for i from 12 to 23 do
		lambdas[i]:=evalf(lambdas[i],20);
	end do;

	for i from 12 to 17 do
		startlevels[i]:=1;
	end do;

	for i from 18 to 23 do
		startlevels[i]:=2;
	end do;
	
	for i from 12 to 23 do
		for j from 0 to 2 do
			T[[j],i]:=0;
		end do;
	end do;

	for i from 12 to 13 do
		for l in [[0,1],[0,2],[1,2]] do
			T[l,i]:=1;
		end do;
	end do;

	for i from 14 to 15 do
		T[[0,1],i]:=0;
		T[[0,2],i]:=-1;
		T[[1,2],i]:=1;
	end do;

	for i from 16 to 17 do
		T[[0,1],i]:=-1;
		T[[0,2],i]:=0;
		T[[1,2],i]:=1;
	end do;

	T[[0,1],18]:=0;
	T[[0,2],18]:=0;
	T[[1,2],18]:=2;
	T[[0,0,1],18]:=0;
	T[[0,0,2],18]:=0;
	T[[0,1,2],18]:=0;
	T[[1,0,1],18]:=1;
	T[[1,0,2],18]:=-1;
	T[[1,1,2],18]:=-1;
	T[[2,0,1],18]:=-1;
	T[[2,0,2],18]:=1;
	T[[2,1,2],18]:=-1;

	
	T[[0,1],19]:=2;
	T[[0,2],19]:=0;
	T[[1,2],19]:=0;
	T[[0,0,1],19]:=-1;
	T[[0,0,2],19]:=1;
	T[[0,1,2],19]:=-1;
	T[[1,0,1],19]:=-1;
	T[[1,0,2],19]:=-1;
	T[[1,1,2],19]:=1;
	T[[2,0,1],19]:=0;
	T[[2,0,2],19]:=0;
	T[[2,1,2],19]:=0;


	T[[0,1],20]:=0;
	T[[0,2],20]:=2;
	T[[1,2],20]:=0;
	T[[0,0,1],20]:=1;
	T[[0,0,2],20]:=-1;
	T[[0,1,2],20]:=-1;
	T[[1,0,1],20]:=0;
	T[[1,0,2],20]:=0;
	T[[1,1,2],20]:=0;
	T[[2,0,1],20]:=-1;
	T[[2,0,2],20]:=-1;
	T[[2,1,2],20]:=1;

	i:=21;
	T[[0,1],i]:=0;
	T[[0,2],i]:=0;
	T[[1,2],i]:=0;
	T[[0,0,1],i]:=0;
	T[[0,0,2],i]:=0;
	T[[0,1,2],i]:=0;
	T[[1,0,1],i]:=-1;
	T[[1,0,2],i]:=1;
	T[[1,1,2],i]:=0;
	T[[2,0,1],i]:=-1;
	T[[2,0,2],i]:=1;
	T[[2,1,2],i]:=0;

	i:=22;
	T[[0,1],i]:=0;
	T[[0,2],i]:=0;
	T[[1,2],i]:=0;
	T[[0,0,1],i]:=0;
	T[[0,0,2],i]:=-1;
	T[[0,1,2],i]:=1;
	T[[1,0,1],i]:=0;
	T[[1,0,2],i]:=-1;
	T[[1,1,2],i]:=1;
	T[[2,0,1],i]:=0;
	T[[2,0,2],i]:=0;
	T[[2,1,2],i]:=0;


	i:=23;
	T[[0,1],i]:=0;
	T[[0,2],i]:=0;
	T[[1,2],i]:=0;
	T[[0,0,1],i]:=1;
	T[[0,0,2],i]:=0;
	T[[0,1,2],i]:=-1;
	T[[1,0,1],i]:=0;
	T[[1,0,2],i]:=0;
	T[[1,1,2],i]:=0;
	T[[2,0,1],i]:=1;
	T[[2,0,2],i]:=0;
	T[[2,1,2],i]:=-1;
end;


setfuncs3:=proc(T,lambdas,startlevels)
	local i,j;
	setfuncs(T,lambdas,startlevels);
	for i from 12 to 17 do
		copyfunc(T,1,i,i+12);
		lambdas[i+12]:=lambdas[i];
	end do;
	for i from 21 to 23 do
		copyfunc(T,2,i,i+9);
		lambdas[i+9]:=lambdas[i];
	end do;

	i:=33;
	T[[0],i]:=0;
	T[[1],i]:=0;
	T[[2],i]:=0;
	T[[0,1],i]:=0;
	T[[0,2],i]:=0;
	T[[1,2],i]:=0;
	T[[0,0,1],i]:=0;
	T[[0,0,2],i]:=0;
	T[[0,1,2],i]:=0;
	T[[0,0,0,1],i]:=1;
	T[[0,0,0,2],i]:=2;
	T[[0,0,1,2],i]:=0;
	T[[0,1,0,1],i]:=-1;
	T[[0,1,0,2],i]:=0;
	T[[0,1,1,2],i]:=1;
	T[[0,2,0,1],i]:=0;
	T[[0,2,0,2],i]:=1;
	T[[0,2,1,2],i]:=-1;
	T[[1,0,1],i]:=0;
	T[[1,0,2],i]:=0;
	T[[1,1,2],i]:=0;
	T[[1,0,0,1],i]:=0;
	T[[1,0,0,2],i]:=0;
	T[[1,0,1,2],i]:=0;
	T[[1,1,0,1],i]:=0;
	T[[1,1,0,2],i]:=0;
	T[[1,1,1,2],i]:=0;
	T[[1,2,0,1],i]:=0;
	T[[1,2,0,2],i]:=0;
	T[[1,2,1,2],i]:=0;
	T[[2,0,1],i]:=0;
	T[[2,0,2],i]:=0;
	T[[2,1,2],i]:=0;
	T[[2,0,0,1],i]:=0;
	T[[2,0,0,2],i]:=0;
	T[[2,0,1,2],i]:=0;
	T[[2,1,0,1],i]:=0;
	T[[2,1,0,2],i]:=0;
	T[[2,1,1,2],i]:=0;
	T[[2,2,0,1],i]:=0;
	T[[2,2,0,2],i]:=0;
	T[[2,2,1,2],i]:=0;

	funcRotate(T,3,[1,2,0],i,i+1);
	funcRotate(T,3,[2,0,1],i,i+2);
	


	i:=36;
	T[[0],i]:=0;
	T[[1],i]:=0;
	T[[2],i]:=0;
	T[[0,1],i]:=0;
	T[[0,2],i]:=0;
	T[[1,2],i]:=0;
	T[[0,0,1],i]:=0;
	T[[0,0,2],i]:=0;
	T[[0,1,2],i]:=0;
	T[[0,0,0,1],i]:=0;
	T[[0,0,0,2],i]:=0;
	T[[0,0,1,2],i]:=0;
	T[[0,1,0,1],i]:=0;
	T[[0,1,0,2],i]:=0;
	T[[0,1,1,2],i]:=0;
	T[[0,2,0,1],i]:=0;
	T[[0,2,0,2],i]:=0;
	T[[0,2,1,2],i]:=0;
	T[[1,0,1],i]:=0;
	T[[1,0,2],i]:=0;
	T[[1,1,2],i]:=0;
	T[[1,0,0,1],i]:=0;
	T[[1,0,0,2],i]:=0;
	T[[1,0,1,2],i]:=0;
	T[[1,1,0,1],i]:=-1;
	T[[1,1,0,2],i]:=1;
	T[[1,1,1,2],i]:=0;
	T[[1,2,0,1],i]:=-1;
	T[[1,2,0,2],i]:=1;
	T[[1,2,1,2],i]:=0;
	T[[2,0,1],i]:=0;
	T[[2,0,2],i]:=0;
	T[[2,1,2],i]:=0;
	T[[2,0,0,1],i]:=0;
	T[[2,0,0,2],i]:=0;
	T[[2,0,1,2],i]:=0;
	T[[2,1,0,1],i]:=-1;
	T[[2,1,0,2],i]:=1;
	T[[2,1,1,2],i]:=0;
	T[[2,2,0,1],i]:=-1;
	T[[2,2,0,2],i]:=1;
	T[[2,2,1,2],i]:=0;

	funcRotate(T,3,[1,2,0],i,i+1);
	funcRotate(T,3,[2,0,1],i,i+2);

	
	i:=39;
	T[[0],i]:=0;
	T[[1],i]:=0;
	T[[2],i]:=0;
	T[[0,1],i]:=0;
	T[[0,2],i]:=0;
	T[[1,2],i]:=0;
	T[[0,0,1],i]:=0;
	T[[0,0,2],i]:=0;
	T[[0,1,2],i]:=2;
	T[[0,0,0,1],i]:=0;
	T[[0,0,0,2],i]:=0;
	T[[0,0,1,2],i]:=0;
	T[[0,1,0,1],i]:=1;
	T[[0,1,0,2],i]:=-1;
	T[[0,1,1,2],i]:=-1;
	T[[0,2,0,1],i]:=-1;
	T[[0,2,0,2],i]:=1;
	T[[0,2,1,2],i]:=-1;
	T[[1,0,1],i]:=0;
	T[[1,0,2],i]:=0;
	T[[1,1,2],i]:=0;
	T[[1,0,0,1],i]:=0;
	T[[1,0,0,2],i]:=0;
	T[[1,0,1,2],i]:=0;
	T[[1,1,0,1],i]:=0;
	T[[1,1,0,2],i]:=0;
	T[[1,1,1,2],i]:=0;
	T[[1,2,0,1],i]:=0;
	T[[1,2,0,2],i]:=0;
	T[[1,2,1,2],i]:=0;
	T[[2,0,1],i]:=0;
	T[[2,0,2],i]:=0;
	T[[2,1,2],i]:=0;
	T[[2,0,0,1],i]:=0;
	T[[2,0,0,2],i]:=0;
	T[[2,0,1,2],i]:=0;
	T[[2,1,0,1],i]:=0;
	T[[2,1,0,2],i]:=0;
	T[[2,1,1,2],i]:=0;
	T[[2,2,0,1],i]:=0;
	T[[2,2,0,2],i]:=0;
	T[[2,2,1,2],i]:=0;

	funcRotate(T,3,[1,2,0],i,i+1);
	funcRotate(T,3,[2,0,1],i,i+2);

	i:=42;
	T[[0],i]:=0;
	T[[1],i]:=0;
	T[[2],i]:=0;
	T[[0,1],i]:=0;
	T[[0,2],i]:=0;
	T[[1,2],i]:=0;
	T[[0,0,1],i]:=0;
	T[[0,0,2],i]:=2;
	T[[0,1,2],i]:=0;
	T[[0,0,0,1],i]:=1;
	T[[0,0,0,2],i]:=-1;
	T[[0,0,1,2],i]:=-1;
	T[[0,1,0,1],i]:=0;
	T[[0,1,0,2],i]:=0;
	T[[0,1,1,2],i]:=0;
	T[[0,2,0,1],i]:=-1;
	T[[0,2,0,2],i]:=-1;
	T[[0,2,1,2],i]:=1;
	T[[1,0,1],i]:=0;
	T[[1,0,2],i]:=0;
	T[[1,1,2],i]:=0;
	T[[1,0,0,1],i]:=0;
	T[[1,0,0,2],i]:=0;
	T[[1,0,1,2],i]:=0;
	T[[1,1,0,1],i]:=0;
	T[[1,1,0,2],i]:=0;
	T[[1,1,1,2],i]:=0;
	T[[1,2,0,1],i]:=0;
	T[[1,2,0,2],i]:=0;
	T[[1,2,1,2],i]:=0;
	T[[2,0,1],i]:=0;
	T[[2,0,2],i]:=0;
	T[[2,1,2],i]:=0;
	T[[2,0,0,1],i]:=0;
	T[[2,0,0,2],i]:=0;
	T[[2,0,1,2],i]:=0;
	T[[2,1,0,1],i]:=0;
	T[[2,1,0,2],i]:=0;
	T[[2,1,1,2],i]:=0;
	T[[2,2,0,1],i]:=0;
	T[[2,2,0,2],i]:=0;
	T[[2,2,1,2],i]:=0;

	funcRotate(T,3,[1,2,0],i,i+1);
	funcRotate(T,3,[2,0,1],i,i+2);

	i:=45;
	T[[0],i]:=0;
	T[[1],i]:=0;
	T[[2],i]:=0;
	T[[0,1],i]:=0;
	T[[0,2],i]:=0;
	T[[1,2],i]:=0;
	T[[0,0,1],i]:=2;
	T[[0,0,2],i]:=0;
	T[[0,1,2],i]:=0;
	T[[0,0,0,1],i]:=-1;
	T[[0,0,0,2],i]:=1;
	T[[0,0,1,2],i]:=-1;
	T[[0,1,0,1],i]:=-1;
	T[[0,1,0,2],i]:=-1;
	T[[0,1,1,2],i]:=1;
	T[[0,2,0,1],i]:=0;
	T[[0,2,0,2],i]:=0;
	T[[0,2,1,2],i]:=0;
	T[[1,0,1],i]:=0;
	T[[1,0,2],i]:=0;
	T[[1,1,2],i]:=0;
	T[[1,0,0,1],i]:=0;
	T[[1,0,0,2],i]:=0;
	T[[1,0,1,2],i]:=0;
	T[[1,1,0,1],i]:=0;
	T[[1,1,0,2],i]:=0;
	T[[1,1,1,2],i]:=0;
	T[[1,2,0,1],i]:=0;
	T[[1,2,0,2],i]:=0;
	T[[1,2,1,2],i]:=0;
	T[[2,0,1],i]:=0;
	T[[2,0,2],i]:=0;
	T[[2,1,2],i]:=0;
	T[[2,0,0,1],i]:=0;
	T[[2,0,0,2],i]:=0;
	T[[2,0,1,2],i]:=0;
	T[[2,1,0,1],i]:=0;
	T[[2,1,0,2],i]:=0;
	T[[2,1,1,2],i]:=0;
	T[[2,2,0,1],i]:=0;
	T[[2,2,0,2],i]:=0;
	T[[2,2,1,2],i]:=0;

	funcRotate(T,3,[1,2,0],i,i+1);
	funcRotate(T,3,[2,0,1],i,i+2);

	i:=48;
	T[[0],i]:=0;
	T[[1],i]:=0;
	T[[2],i]:=0;
	T[[0,1],i]:=0;
	T[[0,2],i]:=0;
	T[[1,2],i]:=2;
	T[[0,0,1],i]:=0;
	T[[0,0,2],i]:=0;
	T[[0,1,2],i]:=0;
	T[[0,0,0,1],i]:=0;
	T[[0,0,0,2],i]:=0;
	T[[0,0,1,2],i]:=0;
	T[[0,1,0,1],i]:=0;
	T[[0,1,0,2],i]:=0;
	T[[0,1,1,2],i]:=0;
	T[[0,2,0,1],i]:=0;
	T[[0,2,0,2],i]:=0;
	T[[0,2,1,2],i]:=0;
	T[[1,0,1],i]:=0;
	T[[1,0,2],i]:=0;
	T[[1,1,2],i]:=0;
	T[[1,0,0,1],i]:=0;
	T[[1,0,0,2],i]:=0;
	T[[1,0,1,2],i]:=0;
	T[[1,1,0,1],i]:=0;
	T[[1,1,0,2],i]:=0;
	T[[1,1,1,2],i]:=0;
	T[[1,2,0,1],i]:=1;
	T[[1,2,0,2],i]:=-1;
	T[[1,2,1,2],i]:=-1;
	T[[2,0,1],i]:=0;
	T[[2,0,2],i]:=0;
	T[[2,1,2],i]:=0;
	T[[2,0,0,1],i]:=0;
	T[[2,0,0,2],i]:=0;
	T[[2,0,1,2],i]:=0;
	T[[2,1,0,1],i]:=-1;
	T[[2,1,0,2],i]:=1;
	T[[2,1,1,2],i]:=-1;
	T[[2,2,0,1],i]:=0;
	T[[2,2,0,2],i]:=0;
	T[[2,2,1,2],i]:=0;

	funcRotate(T,3,[1,2,0],i,i+1);
	funcRotate(T,3,[2,0,1],i,i+2);

end;