N=128;%length
M=log2(N)-1;ii=2^(M-1)-1;
LS=0;LS1=0;
fid=fopen('LS_64.txt','r');%original data file
LS=fscanf(fid,'%d',[N/2,N/2]);
fclose(fid);
for cnti=0:ii
   for cntn=1:N/4
      LS1(4*cnti+1,cntn)=LS(2*cnti+1,cntn);
      LS1(4*cnti+1,cntn+N/4)=LS(2*cnti+2,cntn);
      LS1(4*cnti+1,cntn+N/2)=LS(2*cnti+1,cntn+N/4);
      LS1(4*cnti+1,cntn+3*N/4)=LS(2*cnti+2,cntn+N/4);
   end   
   for cntn=1:N/4
      LS1(4*cnti+2,cntn)=LS(2*cnti+1,cntn);
      LS1(4*cnti+2,cntn+N/4)=-LS(2*cnti+2,cntn);
      LS1(4*cnti+2,cntn+N/2)=LS(2*cnti+1,cntn+N/4);
      LS1(4*cnti+2,cntn+3*N/4)=-LS(2*cnti+2,cntn+N/4);
   end   
   for cntn=1:N/4
      LS1(4*cnti+3,cntn)=LS(2*cnti+2,cntn);
      LS1(4*cnti+3,cntn+N/4)=LS(2*cnti+1,cntn);
      LS1(4*cnti+3,cntn+N/2)=LS(2*cnti+2,cntn+N/4);
      LS1(4*cnti+3,cntn+3*N/4)=LS(2*cnti+1,cntn+N/4);
   end   
   for cntn=1:N/4
      LS1(4*cnti+4,cntn)=LS(2*cnti+2,cntn);
      LS1(4*cnti+4,cntn+N/4)=-LS(2*cnti+1,cntn);
      LS1(4*cnti+4,cntn+N/2)=LS(2*cnti+2,cntn+N/4);
      LS1(4*cnti+4,cntn+3*N/4)=-LS(2*cnti+1,cntn+N/4);
   end   
end

fid=fopen('LS_128.txt','w');%destination data file
for i=1:N
   fprintf(fid,'%1d\n',LS1(:,i));
end
fclose(fid);
fid=fopen('LS_128_1.txt','w');%destination data file
for i=1:N
   fprintf(fid,'%1d ',LS1(i,:));
   fprintf(fid,'\n ');
end
fclose(fid);