function s = sumkm(K,m,n,mask,t)
global locm
locm = mask;

s = 0;
dims = repmat([m],[1 n]);
for xind = 1:m^n
  x = ind2coord(xind,dims);
  c = K(xind);
  
  sx = 0;
  %for t=1:n
    y = mu(t,x(t));
    sx = sx + y;
  %end
  
  %xstr = num2str(x-1);
  %xstr = strrep(xstr,' ','');
  %fprintf('%s  :  %2.5f \n',xstr,sx);
  
  s = s + c*sx;
end
return

function y=mu(n,k)
global locm
mm = locm{n};
y = ismember(k,mm);
return