85 auto pdfs =
block->getData<field::GhostLayerField<float, 19>>(
pdfsID);
86 auto force =
block->getData<field::GhostLayerField<float, 3>>(
forceID);
88 auto &rho_0 = this->
rho_0_;
89 WALBERLA_ASSERT_GREATER_EQUAL(0, -int_c(force->nrOfGhostLayers()))
90 float *
RESTRICT const _data_force = force->dataAt(0, 0, 0, 0);
91 WALBERLA_ASSERT_EQUAL(force->layout(), field::fzyx)
92 WALBERLA_ASSERT_GREATER_EQUAL(0, -int_c(pdfs->nrOfGhostLayers()))
93 float *
RESTRICT _data_pdfs = pdfs->dataAt(0, 0, 0, 0);
94 WALBERLA_ASSERT_EQUAL(pdfs->layout(), field::fzyx)
95 WALBERLA_ASSERT_GREATER_EQUAL(0, -int_c(
velocity->nrOfGhostLayers()))
97 WALBERLA_ASSERT_EQUAL(
velocity->layout(), field::fzyx)
98 WALBERLA_ASSERT_GREATER_EQUAL(force->xSizeWithGhostLayer(), int64_t(int64_c(force->xSize()) + 0))
99 const int64_t _size_force_0 = int64_t(int64_c(force->xSize()) + 0);
100 WALBERLA_ASSERT_EQUAL(force->layout(), field::fzyx)
101 WALBERLA_ASSERT_GREATER_EQUAL(force->ySizeWithGhostLayer(), int64_t(int64_c(force->ySize()) + 0))
102 const int64_t _size_force_1 = int64_t(int64_c(force->ySize()) + 0);
103 WALBERLA_ASSERT_EQUAL(force->layout(), field::fzyx)
104 WALBERLA_ASSERT_GREATER_EQUAL(force->zSizeWithGhostLayer(), int64_t(int64_c(force->zSize()) + 0))
105 const int64_t _size_force_2 = int64_t(int64_c(force->zSize()) + 0);
106 WALBERLA_ASSERT_EQUAL(force->layout(), field::fzyx)
107 const int64_t _stride_force_0 = int64_t(force->xStride());
108 const int64_t _stride_force_1 = int64_t(force->yStride());
109 const int64_t _stride_force_2 = int64_t(force->zStride());
110 const int64_t _stride_force_3 = int64_t(1 * int64_t(force->fStride()));
111 const int64_t _stride_pdfs_0 = int64_t(pdfs->xStride());
112 const int64_t _stride_pdfs_1 = int64_t(pdfs->yStride());
113 const int64_t _stride_pdfs_2 = int64_t(pdfs->zStride());
114 const int64_t _stride_pdfs_3 = int64_t(1 * int64_t(pdfs->fStride()));
115 const int64_t _stride_velocity_0 = int64_t(
velocity->xStride());
116 const int64_t _stride_velocity_1 = int64_t(
velocity->yStride());
117 const int64_t _stride_velocity_2 = int64_t(
velocity->zStride());
118 const int64_t _stride_velocity_3 = int64_t(1 * int64_t(
velocity->fStride()));
119 internal_b8085d63d6b7e842485134abbac511e8::initialpdfssettersingleprecision_initialpdfssettersingleprecision(_data_force, _data_pdfs, _data_velocity, _size_force_0, _size_force_1, _size_force_2, _stride_force_0, _stride_force_1, _stride_force_2, _stride_force_3, _stride_pdfs_0, _stride_pdfs_1, _stride_pdfs_2, _stride_pdfs_3, _stride_velocity_0, _stride_velocity_1, _stride_velocity_2, _stride_velocity_3, rho_0);
124 CellInterval ci = globalCellInterval;
125 CellInterval blockBB = blocks->getBlockCellBB(*
block);
126 blockBB.expand(ghostLayers);
127 ci.intersect(blockBB);
128 blocks->transformGlobalToBlockLocalCellInterval(ci, *
block);
133 auto pdfs =
block->getData<field::GhostLayerField<float, 19>>(
pdfsID);
134 auto force =
block->getData<field::GhostLayerField<float, 3>>(
forceID);
136 auto &rho_0 = this->
rho_0_;
137 WALBERLA_ASSERT_GREATER_EQUAL(ci.xMin(), -int_c(force->nrOfGhostLayers()))
138 WALBERLA_ASSERT_GREATER_EQUAL(ci.yMin(), -int_c(force->nrOfGhostLayers()))
139 WALBERLA_ASSERT_GREATER_EQUAL(ci.zMin(), -int_c(force->nrOfGhostLayers()))
140 float *
RESTRICT const _data_force = force->dataAt(ci.xMin(), ci.yMin(), ci.zMin(), 0);
141 WALBERLA_ASSERT_EQUAL(force->layout(), field::fzyx)
142 WALBERLA_ASSERT_GREATER_EQUAL(ci.xMin(), -int_c(pdfs->nrOfGhostLayers()))
143 WALBERLA_ASSERT_GREATER_EQUAL(ci.yMin(), -int_c(pdfs->nrOfGhostLayers()))
144 WALBERLA_ASSERT_GREATER_EQUAL(ci.zMin(), -int_c(pdfs->nrOfGhostLayers()))
145 float *
RESTRICT _data_pdfs = pdfs->dataAt(ci.xMin(), ci.yMin(), ci.zMin(), 0);
146 WALBERLA_ASSERT_EQUAL(pdfs->layout(), field::fzyx)
147 WALBERLA_ASSERT_GREATER_EQUAL(ci.xMin(), -int_c(
velocity->nrOfGhostLayers()))
148 WALBERLA_ASSERT_GREATER_EQUAL(ci.yMin(), -int_c(
velocity->nrOfGhostLayers()))
149 WALBERLA_ASSERT_GREATER_EQUAL(ci.zMin(), -int_c(
velocity->nrOfGhostLayers()))
150 float *
RESTRICT const _data_velocity =
velocity->dataAt(ci.xMin(), ci.yMin(), ci.zMin(), 0);
151 WALBERLA_ASSERT_EQUAL(
velocity->layout(), field::fzyx)
152 WALBERLA_ASSERT_GREATER_EQUAL(force->xSizeWithGhostLayer(), int64_t(int64_c(ci.xSize()) + 0))
153 const int64_t _size_force_0 = int64_t(int64_c(ci.xSize()) + 0);
154 WALBERLA_ASSERT_EQUAL(force->layout(), field::fzyx)
155 WALBERLA_ASSERT_GREATER_EQUAL(force->ySizeWithGhostLayer(), int64_t(int64_c(ci.ySize()) + 0))
156 const int64_t _size_force_1 = int64_t(int64_c(ci.ySize()) + 0);
157 WALBERLA_ASSERT_EQUAL(force->layout(), field::fzyx)
158 WALBERLA_ASSERT_GREATER_EQUAL(force->zSizeWithGhostLayer(), int64_t(int64_c(ci.zSize()) + 0))
159 const int64_t _size_force_2 = int64_t(int64_c(ci.zSize()) + 0);
160 WALBERLA_ASSERT_EQUAL(force->layout(), field::fzyx)
161 const int64_t _stride_force_0 = int64_t(force->xStride());
162 const int64_t _stride_force_1 = int64_t(force->yStride());
163 const int64_t _stride_force_2 = int64_t(force->zStride());
164 const int64_t _stride_force_3 = int64_t(1 * int64_t(force->fStride()));
165 const int64_t _stride_pdfs_0 = int64_t(pdfs->xStride());
166 const int64_t _stride_pdfs_1 = int64_t(pdfs->yStride());
167 const int64_t _stride_pdfs_2 = int64_t(pdfs->zStride());
168 const int64_t _stride_pdfs_3 = int64_t(1 * int64_t(pdfs->fStride()));
169 const int64_t _stride_velocity_0 = int64_t(
velocity->xStride());
170 const int64_t _stride_velocity_1 = int64_t(
velocity->yStride());
171 const int64_t _stride_velocity_2 = int64_t(
velocity->zStride());
172 const int64_t _stride_velocity_3 = int64_t(1 * int64_t(
velocity->fStride()));
173 internal_b8085d63d6b7e842485134abbac511e8::initialpdfssettersingleprecision_initialpdfssettersingleprecision(_data_force, _data_pdfs, _data_velocity, _size_force_0, _size_force_1, _size_force_2, _stride_force_0, _stride_force_1, _stride_force_2, _stride_force_3, _stride_pdfs_0, _stride_pdfs_1, _stride_pdfs_2, _stride_pdfs_3, _stride_velocity_0, _stride_velocity_1, _stride_velocity_2, _stride_velocity_3, rho_0);
static FUNC_PREFIX void initialpdfssettersingleprecision_initialpdfssettersingleprecision(float *RESTRICT const _data_force, float *RESTRICT _data_pdfs, float *RESTRICT const _data_velocity, int64_t const _size_force_0, int64_t const _size_force_1, int64_t const _size_force_2, int64_t const _stride_force_0, int64_t const _stride_force_1, int64_t const _stride_force_2, int64_t const _stride_force_3, int64_t const _stride_pdfs_0, int64_t const _stride_pdfs_1, int64_t const _stride_pdfs_2, int64_t const _stride_pdfs_3, int64_t const _stride_velocity_0, int64_t const _stride_velocity_1, int64_t const _stride_velocity_2, int64_t const _stride_velocity_3, float rho_0)