% junk

function Lpd = getpairwisedif(m,n)
T = n-1;

adim = m^2*T;

Lpd = zeros(0,adim);
rind=0;
for t=1:T
  for i1=1:m
    for i2=i1+1:m
      rind = rind+1;
      for j=1:m
        ind1 = coord2ind([j i1 t],[m m T]);
        ind2 = coord2ind([j i2 t],[m m T]);
        Lpd(rind,ind1) =  1;
        Lpd(rind,ind2) = -1;
      end
    end  
  end
end
return