Add ofdft_fe
Use it from ofmd.f90
Converge SCF to 1e-9, correct final energies
Use the same free_energy_min()

And expose Nq, quad_type and energy_eps.