  function [output]=xops(W,n)
% function output=xops(x,w,n)
%Xu is the symmetric data while w is the general input to xops
%    Xu is usMat(x)
global Xu;
output=sum((Xu.*W),2) + W(1:n+1:end)';    % column vector
