37 std::vector<vel_type> local_velocities{};
41 lb.ghost_communication_vel();
44 if (
auto const vel = lb.get_interpolated_velocity(pos)) {
45 local_positions.emplace_back(pos);
46 local_velocities.emplace_back((*vel) * vel_conv);
50 auto const [global_positions, global_velocities] =
51 detail::gather(comm, local_positions, local_velocities);
53 if (comm.rank() != 0) {
58 detail::accumulate(histogram, global_positions, global_velocities);
61 }
catch (detail::empty_bin_exception
const &) {
62 throw std::runtime_error(
63 "Decrease sampling delta(s), some bins have no hit");