read "itall.txt":


TESTFN:=12:

MAKEBIHARMONIC:=false:
	BLAP1:=0:	
	BLAP2:=50:
	BLAP3:=-100:
MAKESQUAREBIHARMONIC:=true:
MAKESQUAREEIGEN:=false:
	SELVL:=4:
	SENO:=95:
MAKEEIGEN:=false:
	TELVL:=5:
	TENO:=349:
APPFUNCTION:=true:
	APPFN:=13:
	TMPFN:=14:
lvl:=3:
phatlvl:=8:
path:="../data2/":

T:=gasket(phatlvl);

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

if (MAKEBIHARMONIC) then
	T[[0],8]:=0:		
	T[[1],8]:=0:
	T[[2],8]:=0:
	T[[0],9]:=BLAP1;
	T[[1],9]:=BLAP2;
	T[[2],9]:=BLAP3;
	
	n_harm(T,R,2,phatlvl,[0],[1],[2],8):
	copyFunc(T,phatlvl,8,TESTFN):

	print("Done Making Bi-Harmonic Function");
end if;

if (MAKESQUAREBIHARMONIC) then
	T[[0],8]:=0:		
	T[[1],8]:=0:
	T[[2],8]:=0:
	T[[0],9]:=BLAP1;
	T[[1],9]:=BLAP2;
	T[[2],9]:=BLAP3;
	
	n_harm(T,R,2,phatlvl,[0],[1],[2],8):
	copyFunc(T,phatlvl,8,TESTFN):
	funcProduct(T,phatlvl,TESTFN,TESTFN,TESTFN):	

	print("Done Making Square Bi-Harmonic Function");
end if;

if (MAKESQUAREEIGEN) then
	readFunc(T,TESTFN,cat(path,"level",SELVL,"/",phatlvl,"phat",SELVL,"eigen.",SENO));
	funcProduct(T,phatlvl,TESTFN,TESTFN,TESTFN);
end if:

if (MAKEEIGEN) then
	readFunc(T,TESTFN,cat(path,"level",TELVL,"/",phatlvl,"phat",TELVL,"eigen.",TENO));
end if:

if (APPFUNCTION) then
	clearFunc(T,phatlvl,APPFN);
	for i from 1 to nops(verts) do
		co:=T[verts[i],TESTFN];
		readFunc(T,TMPFN,cat(path,"level",lvl,"/",phatlvl,"phat",lvl,"psi.",i));
		funcSumWithConstant(T,phatlvl,co,APPFN,TMPFN,APPFN);
	end do;
end if:

print("test function");
ourPlotBlack(T,phatlvl,TESTFN,"Test Function");
print(cat("Max=",findFuncMax(T,phatlvl,TESTFN),", Min=",findFuncMin(T,phatlvl,TESTFN)));
print("approximation of test function");
ourPlotBlack(T,phatlvl,APPFN,"Approximation of Test Function");
print(cat("Max=",findFuncMax(T,phatlvl,APPFN),", Min=",findFuncMin(T,phatlvl,APPFN)));
funcDifference(T,phatlvl,APPFN,TESTFN,TMPFN);
print("approximation - test function");
ourPlotBlack(T,phatlvl,TMPFN,"Approximation - Test Function");
print(cat("Max=",findFuncMax(T,phatlvl,TMPFN),", Min=",findFuncMin(T,phatlvl,TMPFN)));

print(cat("level= ",lvl));