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