read "everything.txt":
read "../setfuncs.txt":
read "../algo24.mws":
with(linalg):
read "../domatrix.txt":
read "../func_file.txt":
read "../funcarith.txt":
T:=gasket(10):
setfuncs(T,startlambdas,startlevels):
for i from 12 to 17 do
algo24Wrapper(T,i,startlevels[i],startlambdas[i],makeEpsilons(8)):
saveFunc(T,i,10,"../output/psi10");
end do;
for i from 18 to 23 do
algo24Wrapper(T,i,startlevels[i],startlambdas[i],makeEpsilons(7)):
saveFunc(T,i,10,"../output/psi10");
end do;
makeAllPsi(T,2,10);
for i from 24 to 35 do
saveFunc(T,i,10,"../output/psi10");
end do;
for i from 24 to 35 do
for j from i to 35 do
funcProduct(T,10,i,j,36);
err:=simp_rule(T,10,36);
fd:=fopen("../output/psiinnerprod.txt",APPEND);
fprintf(fd,"%d %d %20.20g\n",i,j,err);
fclose(fd);
end do;
end do;
quit