function [Ainf,Binf] = lipabcw(m,n,CW)
mn = m^n;
L = mn;
Ainf = zeros(L,mn);
Binf = zeros(L,1);
MN = repmat([m],[1 n]);
rind = 0;
for ind1 = 1:mn
  for ind2 = ind1+1:mn
    x = ind2coord(ind1,MN);
    y = ind2coord(ind2,MN);
    ii = find(x~=y);
    if length(ii)==1  % Hamming dist == 1
      Cii = CW(:,:,ii);
      cxy = Cii(x(ii),y(ii));
      for s = [-1,1]
        rind = rind+1;
        Ainf(rind,ind1) =  s;
        Ainf(rind,ind2) = -s;
        Binf(rind,1) = cxy;
      end
    end
  end
end
return