Y=X(pv,:); J1=pdist2(Y,Y);

J1=J0(pv,pv);