function H = totstrict(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));
    H = H + d;
  end
end
return