Fixed virial when biasing ENERGY

Virial was not scaled together with forces. The contribution is not
negligible, and I verified it is correct by testing with gromacs
a box of water in NPT and biasing the energy, then comparing the
-kt*log(histogram) biased and unbiased. It seems after this fix agreement
is much better