fragment"><pre>0001 <a name="_sub0" href="#_subfunctions" class="code">function [ Xhat ] = Axz( y,xhat,zhat,T )</a>
0002 <span class="comment">%</span>
0003 <span class="comment">% Axz operator A_{xz}(x) \equiv A(x \otimes y \otimes z) fixing variables x</span>
0004 <span class="comment">% and z</span>
0005 
0006 n=sqrt(size(xhat,1));m=n;
0007 yhat=1/n*fft2(reshape(y,n,m));yhat=yhat(:);
0008 
0009 Xhat=zeros(size(T,1),1);
0010 
0011 
0012 <span class="keyword">for</span> i=1:size(T,1)
0013     
0014     Xha