In runner_doself1_density_vec() do not check for r2 > 0. This is never done in the scalar code either.