function [h,hh] = neostrict2(P,m,n)
h = 0;
hh = [];
for t=1:n-1
  hh(end+1) = strictt(P,m,n,t);
end
h = max(hh);
return

function h = strictt(P,m,n,t)
PP = [];
dimt = repmat([m],[1 t]);
for xi = 1:m^t
  X = ind2coord(xi,dimt);
  Pt = getcondP(P,m,n,t+1,1:t,X);
  PP = [PP Pt];
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