Perhaps a lumped-mass projection (as demonstrated in the function lumped_project here) of the cell-wise energy to FunctionSpace(mesh,"CG",1) would be a good candidate. That would produce a sort of weighted average of the values in cells surrounding each mesh vertex.
1 Like