function a = correlate(labels1, labels2)
  a = op(labels1,'==',min(labels1):max(labels1))' ...
      * op(labels2,'==',min(labels2):max(labels2));
  