Improved check of if only forces are computed in do_pairs()