read "itall.txt":

lvl:=3;
numeigens:=(3^(lvl+1)-3)/2;
path:="../data/";
T:=gasket(lvl);

fd:=fopen(cat(path,lvl,"squarednorms.txt"),WRITE);

for i from 1 to numeigens do
        readFunc(T,12,cat(path,lvl,"eigen.",i));
        fprintf(fd,"%20.20g\n",discInnerProduct(T,lvl,12,12));
        saveFunc(T,12,lvl,cat(path,lvl,"eigen.",i));
end do;

fclose(fd);