Update how P is initialized
Put mu.setb() back--mistakenly deleted in last commmit
Make d0Hist work with arbitrary nf