debug: ignoring regularization terms in gradients (mf.pmf)
mf.pmf: eliminate np.float when integer division
mf.pmf: update test
+ 3.1.0