Double checking the Elser Processor

proxtoolbox/experiments/phase/Elser_processor.m

@@ -30,4 +30,14 @@ for(i=M/2+1:M)
         Fmag(j,i) = fmag(i,j);
     end
 end
-Fmag(1,M/2:M)=fmag(M/2:M,1);
\ No newline at end of file
+Fmag(1,M/2:M)=fmag(M/2:M,1);
+
+% Now try this a different way:
+Fmag2=zeros(M,M);
+fmag2=sqrt(x)/M;
+fmag2(:,M/2+1)=0;
+Fmag2(:,1:M/2+1)=fmag2;
+Fmag2(2:M/2,M/2+2:M)=fmag2(M:-1:M/2+2,M/2-1:-1:1);
+Fmag2(M/2+2:M,M/2+2:M)=fmag2(M/2:-1:2,M/2:-1:2);
+Fmag2(1,M/2:M)=fmag2(M/2:M,1);
+datapow2=sum(sum(Fmag2.^2))*M;