function phi = fillgmask(mask,m)
n = length(mask);
phi = zeros(1,m^n);
dims = repmat([m],[1 n]);
for xind = 1:m^n
  x = ind2coord(xind,dims);
  phi(xind) = applymask(mask,x);
end
return


function r = applymask(mask,x)
n = length(x);
r = 0;
for i=1:n
  r = r + ismember(x(i),mask{i});
end
return