22#include "system/System.hpp"
33 boost::mpi::communicator
const &comm)
const {
37 std::vector<vel_type> local_velocities{};
43 if (
auto const vel = lb.get_interpolated_velocity(
pos)) {
51 local_positions.emplace_back(pos_cyl);
52 local_velocities.emplace_back(vel_cyl);
56 auto const [global_positions, global_velocities] =
57 detail::gather(comm, local_positions, local_velocities);
59 if (comm.rank() != 0) {
64 detail::accumulate(histogram, global_positions, global_velocities);
65 return detail::normalize_by_bin_size(histogram);
__shared__ int pos[MAXDEPTH *THREADS5/WARPSIZE]
std::vector< Utils::Vector3d > sampling_positions
std::vector< double > operator()(boost::mpi::communicator const &comm) const override
std::shared_ptr< Utils::CylindricalTransformationParameters > transform_params
Histogram in cylindrical coordinates.
Vector3d transform_vector_cartesian_to_cylinder(Vector3d const &vec, Vector3d const &axis, Vector3d const &pos)
Vector transformation from Cartesian to cylindrical coordinates.
Vector3d transform_coordinate_cartesian_to_cylinder(Vector3d const &pos)
Coordinate transformation from Cartesian to cylindrical coordinates.
auto get_lattice_speed() const
Get the lattice speed (agrid/tau).