Numerical minimisation of Gutzwiller energy functionals