Output energy_density causes a segfault

The output energy_density causes a segfault when computing the xc contribution.

This output is also not tested in the testsuite at the moment. A test should be added when this bug is fixed.