If total energy (not energy density) is returned, by convention Exc should be reduced over MPI ranks before returning.
If total energy (not energy density) is returned, by convention Exc should be reduced over MPI ranks before returning.