function output=zops(Qy,W,n)
%          output=zops(Qy,W,n);
%   Qy=diag(y)-Q  (sparse)
global indsu ;

ZW=(Qy/sqrt(2))*W;
ZWp=ZW';
output=ZW(indsu)+ZWp(indsu);