33 boost::mpi::communicator
const &comm)
const {
37 std::vector<vel_type> local_velocities{};
41 lb.ghost_communication_vel();
44 if (
auto const vel = lb.get_interpolated_velocity(pos)) {
52 local_positions.emplace_back(pos_cyl);
53 local_velocities.emplace_back(vel_cyl);
57 auto const [global_positions, global_velocities] =
58 detail::gather(comm, local_positions, local_velocities);
60 if (comm.rank() != 0) {
65 detail::accumulate(histogram, global_positions, global_velocities);
66 return detail::normalize_by_bin_size(histogram);