read "itall.txt";

newNormFactor:=proc(oldlambda,newlambda)
	local factor, num,den;
	
	den:=(2-newlambda)^2*(5-newlambda)^2;
	num:=4*(4-newlambda)^2 + 8 + (6-newlambda)*(4-newlambda)*(4-oldlambda);

	factor:=(num/den)+1;

	return(factor);
end;

newestNormFactor:=proc(newlambda)
	local factor, num, den;

	num:=(newlambda^2 - 10*newlambda + 20);
	den:=(5 - newlambda)*(2-newlambda);
	factor:=1+num/den;
	
	return(factor);
end;


T:=gasket(7);


readFunc(T,12,"../data/ortheigen2level1.1");
newAlgoWrapper(T,1,12,2,[1,-1,1,-1,1]);

newlambda:=makeLambda(2,[1]);

p:=discInnerProduct(T,1,12,12);
evalf(p*newestNormFactor(newlambda));
discInnerProduct(T,2,12,12);

newlambda:=makeLambda(2,[1,-1]);

p:=discInnerProduct(T,2,12,12);
evalf(p*newestNormFactor(newlambda));
discInnerProduct(T,3,12,12);

newlambda:=makeLambda(2,[1,-1,1]);

p:=discInnerProduct(T,3,12,12);
evalf(p*newestNormFactor(newlambda));
discInnerProduct(T,4,12,12);

newlambda:=makeLambda(2,[1,-1,1,-1]);

p:=discInnerProduct(T,4,12,12);
evalf(p*newestNormFactor(newlambda));
discInnerProduct(T,5,12,12);

newlambda:=makeLambda(2,[1,-1,1,-1,1]);

p:=discInnerProduct(T,5,12,12);
evalf(p*newestNormFactor(newlambda));
discInnerProduct(T,6,12,12);