The function `SystemState::conservation_energy` needs to be rewritten to use the cabana particle data structure and cabana parallel for abstraction. It may be the case to look if Cabana and/or kokkos has a "parallel reduce" for this.