function [h,hh] = genstrict(P,m,n)
%P = P + 1e-12;
%P = P/sum(P);
hh=[];
for i=1:n
  hh(end+1) = stricti(P,m,n,i);
end
h = max(hh);
return

function h = stricti(P,m,n,i)
ii = setdiff(1:n,i);
PP = [];

dims = repmat([m],[1 n-1]);
for xi = 1:m^(n-1)
  X = ind2coord(xi,dims);
  [pX,Xinds] = getPX(P,m,n,X,ii);
  Pi = P(Xinds)/(pX+realmin);
  PP = [PP Pi];
end
h = strict(PP);
return


function h = strict(A)
[n,n] = size(A);
h = 0;
for i=1:n
  for j=i+1:n
    x = A(:,i);
    y = A(:,j);
    d = .5*sum(abs(x-y));
    if d>h
      h=d;
    end
  end
end
return