function h = partstrict(A,t,i1,i2)
x = A(:,i1,t);
y = A(:,i2,t);
h = .5*sum(abs(x-y));
return