r11=f(iw(1,1),lw(1,1))+f(iw(2,1),lw(1,2))+f(iw(3,1),lw(1,3))+f(iw(4,1),lw(1,4));
 R11=ff(r11);
 r21=f(iw(1,2),lw(1,1))+f(iw(2,2),lw(1,2))+f(iw(3,2),lw(1,3))+f(iw(4,2),lw(1,4));
 R21=ff(r21);
 r31=f(iw(1,3),lw(1,1))+f(iw(2,3),lw(1,2))+f(iw(3,3),lw(1,3))+f(iw(4,3),lw(1,4));
 R31=ff(r31);
 r41=f(iw(1,4),lw(1,1))+f(iw(2,4),lw(1,2))+f(iw(3,4),lw(1,3))+f(iw(4,4),lw(1,4));
 R41=ff(r41);
 r51=f(iw(1,5),lw(1,1))+f(iw(2,5),lw(1,2))+f(iw(3,5),lw(1,3))+f(iw(4,5),lw(1,4));
 R51=ff(r51);
 r61=f(iw(1,6),lw(1,1))+f(iw(2,6),lw(1,2))+f(iw(3,6),lw(1,3))+f(iw(4,6),lw(1,4));
 R61=ff(r61);
 r71=f(iw(1,7),lw(1,1))+f(iw(2,7),lw(1,2))+f(iw(3,7),lw(1,3))+f(iw(4,7),lw(1,4));
 R71=ff(r71);
 r81=f(iw(1,8),lw(1,1))+f(iw(2,8),lw(1,2))+f(iw(3,8),lw(1,3))+f(iw(4,8),lw(1,4));
 R81=ff(r81);
 S11=fff(R11,R11,R21,R31,R41,R51,R61,R71,R81);
 S21=fff(R21,R11,R21,R31,R41,R51,R61,R71,R81);
 S31=fff(R31,R11,R21,R31,R41,R51,R61,R71,R81);
 S41=fff(R41,R11,R21,R31,R41,R51,R61,R71,R81);
 S51=fff(R51,R11,R21,R31,R41,R51,R61,R71,R81);
 S61=fff(R61,R11,R21,R31,R41,R51,R61,R71,R81);
 S71=fff(R71,R11,R21,R31,R41,R51,R61,R71,R81);
 S81=fff(R81,R11,R21,R31,R41,R51,R61,R71,R81);
 S=[S11,S21,S31,S41,S51,S61,S71,S81]