All gpart now find all of their interactions following correct boundary condition check in runner_do_grav_long_range()