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;

T:=gasket(7);


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

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

newNormFactor(oldlambda,newlambda);

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

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

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

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

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

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

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

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

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