function linmask(phi,m,n)
phi = round(phi);
phi = phi-min(phi);

dims = repmat([m],[1 n]);


A = zeros(m^n,m*n);
B = phi;

for i=1:m^n
  x = ind2coord(i,dims);
  xstr = [];
  for j=1:n
    xj = zeros(1,m);
    xj(x(j)) = 1;
    xstr = [xstr xj];
  end
  A(i,:) = xstr;
end

keyboard
return