working on PF, haven't broken gibbs yet
PF now working for single particle, checked against reference DM