  function output=xop(dy,n)
% function output=xop(dy,n);
%    Xu is usMat(x)
%global Xu Q q2kron indsu ;
global Xu;
%D=spdiags(dy,0,n,n);
D=sparse(diag(dy));
output= D*Xu+D;
output=output(:);
