Digits:=20:
read "everything.txt":
read "func_file.txt":
read "funcarith.txt":
read "autoset.txt":

lvl:=4;
path:="../data/";
numEigens:=(3^(lvl+1) -3)/2;
startEigen:=12;
tol1:=.0000000001;
tol2:=.0000000001;
tol0:=.0000000001;

fd:=fopen(cat(path,lvl,"lambdalist.txt"), READ);
#fscanf(fd,"%g\n")[1]


for i from 1 to numEigens do
	readFunc(T,i+startEigen-1,cat(path,lvl,"eigen.",i));
end do;

if 1=1
then
for i from 1 to numEigens do
	for j from i to numEigens do
		m:=discInnerProduct(T,lvl,i+startEigen-1,j+startEigen-1);
                #fd:=fopen(cat(path,"neigensvalues.txt"),APPEND);
			
		if (i=j) then
			if(abs(m-1) > tol0) then
			printf("(%d,%d)=%20.20g\n",i,j,m);
			end if;
		else 	if(abs(m) > tol1) then                
			printf("(%d,%d)=%20.20g\n",i,j,m);
			end if;
		end if;

                #fclose(fd);
	end do;
end do;
end if;

for i from 1 to numEigens do
	discLaplac(T,lvl,i+startEigen-1,startEigen+numEigens);
	funcQuotient(T,lvl,startEigen+numEigens,i+startEigen-1,startEigen+numEigens);
	lambda:=fscanf(fd,"%g\n")[1];
	lambda:=-1*lambda;
	m:=compareFuncList(T,lvl,startEigen+numEigens,[0,lambda]);
	#m:=findFuncMaxDifferenceWithoutZero(T,lvl,startEigen+numEigens,.0000000001);
        #fd:=fopen(cat(path,"neigensvalues.txt"),APPEND);
        
	if (abs(m) > tol2) then 
		printf("max variance of (laplacian(%d)/%d) =%20.20g\n",i,i,m);
        end if;

	#fclose(fd);
end do;

fclose(fd);

print("All done");