function [yop] = y_operator(data)

sz = size(data);
ii = (sz(2)/2);

yop = ([1:sz]'*0+1) * [-ii:(ii-1)];
