Compute the moments and .
psi(phi_1) - psi(sum_d phi_{1,d})
Standard gradient given the gradient with respect to the moments, that is, given the Riemannian gradient :