66 auto rho_0 =
block->getData<field::GhostLayerField<double, 1>>(
rho_0ID);
67 auto rho_1 =
block->getData<field::GhostLayerField<double, 1>>(
rho_1ID);
74 WALBERLA_ASSERT_GREATER_EQUAL(0, -int_c(rho_0->nrOfGhostLayers()))
75 double *
RESTRICT _data_rho_0 = rho_0->dataAt(0, 0, 0, 0);
76 WALBERLA_ASSERT_GREATER_EQUAL(0, -int_c(rho_1->nrOfGhostLayers()))
77 double *
RESTRICT _data_rho_1 = rho_1->dataAt(0, 0, 0, 0);
78 WALBERLA_ASSERT_GREATER_EQUAL(rho_0->xSizeWithGhostLayer(), int64_t(int64_c(rho_0->xSize()) + 0))
79 const int64_t _size_rho_0_0 = int64_t(int64_c(rho_0->xSize()) + 0);
80 WALBERLA_ASSERT_GREATER_EQUAL(rho_0->ySizeWithGhostLayer(), int64_t(int64_c(rho_0->ySize()) + 0))
81 const int64_t _size_rho_0_1 = int64_t(int64_c(rho_0->ySize()) + 0);
82 WALBERLA_ASSERT_GREATER_EQUAL(rho_0->zSizeWithGhostLayer(), int64_t(int64_c(rho_0->zSize()) + 0))
83 const int64_t _size_rho_0_2 = int64_t(int64_c(rho_0->zSize()) + 0);
84 const int64_t _stride_rho_0_0 = int64_t(rho_0->xStride());
85 const int64_t _stride_rho_0_1 = int64_t(rho_0->yStride());
86 const int64_t _stride_rho_0_2 = int64_t(rho_0->zStride());
87 const int64_t _stride_rho_1_0 = int64_t(rho_1->xStride());
88 const int64_t _stride_rho_1_1 = int64_t(rho_1->yStride());
89 const int64_t _stride_rho_1_2 = int64_t(rho_1->zStride());
90 internal_2cb10021ef8890fa965cb94996ae1510::reactionkernelbulk_2_double_precision_reactionkernelbulk_2_double_precision(_data_rho_0, _data_rho_1, _size_rho_0_0, _size_rho_0_1, _size_rho_0_2, _stride_rho_0_0, _stride_rho_0_1, _stride_rho_0_2, _stride_rho_1_0, _stride_rho_1_1, _stride_rho_1_2, order_0, order_1, rate_coefficient, stoech_0, stoech_1);
95 CellInterval ci = globalCellInterval;
96 CellInterval blockBB = blocks->getBlockCellBB(*
block);
97 blockBB.expand(ghostLayers);
98 ci.intersect(blockBB);
99 blocks->transformGlobalToBlockLocalCellInterval(ci, *
block);
103 auto rho_0 =
block->getData<field::GhostLayerField<double, 1>>(
rho_0ID);
104 auto rho_1 =
block->getData<field::GhostLayerField<double, 1>>(
rho_1ID);
111 WALBERLA_ASSERT_GREATER_EQUAL(ci.xMin(), -int_c(rho_0->nrOfGhostLayers()))
112 WALBERLA_ASSERT_GREATER_EQUAL(ci.yMin(), -int_c(rho_0->nrOfGhostLayers()))
113 WALBERLA_ASSERT_GREATER_EQUAL(ci.zMin(), -int_c(rho_0->nrOfGhostLayers()))
114 double *
RESTRICT _data_rho_0 = rho_0->dataAt(ci.xMin(), ci.yMin(), ci.zMin(), 0);
115 WALBERLA_ASSERT_GREATER_EQUAL(ci.xMin(), -int_c(rho_1->nrOfGhostLayers()))
116 WALBERLA_ASSERT_GREATER_EQUAL(ci.yMin(), -int_c(rho_1->nrOfGhostLayers()))
117 WALBERLA_ASSERT_GREATER_EQUAL(ci.zMin(), -int_c(rho_1->nrOfGhostLayers()))
118 double *
RESTRICT _data_rho_1 = rho_1->dataAt(ci.xMin(), ci.yMin(), ci.zMin(), 0);
119 WALBERLA_ASSERT_GREATER_EQUAL(rho_0->xSizeWithGhostLayer(), int64_t(int64_c(ci.xSize()) + 0))
120 const int64_t _size_rho_0_0 = int64_t(int64_c(ci.xSize()) + 0);
121 WALBERLA_ASSERT_GREATER_EQUAL(rho_0->ySizeWithGhostLayer(), int64_t(int64_c(ci.ySize()) + 0))
122 const int64_t _size_rho_0_1 = int64_t(int64_c(ci.ySize()) + 0);
123 WALBERLA_ASSERT_GREATER_EQUAL(rho_0->zSizeWithGhostLayer(), int64_t(int64_c(ci.zSize()) + 0))
124 const int64_t _size_rho_0_2 = int64_t(int64_c(ci.zSize()) + 0);
125 const int64_t _stride_rho_0_0 = int64_t(rho_0->xStride());
126 const int64_t _stride_rho_0_1 = int64_t(rho_0->yStride());
127 const int64_t _stride_rho_0_2 = int64_t(rho_0->zStride());
128 const int64_t _stride_rho_1_0 = int64_t(rho_1->xStride());
129 const int64_t _stride_rho_1_1 = int64_t(rho_1->yStride());
130 const int64_t _stride_rho_1_2 = int64_t(rho_1->zStride());
131 internal_2cb10021ef8890fa965cb94996ae1510::reactionkernelbulk_2_double_precision_reactionkernelbulk_2_double_precision(_data_rho_0, _data_rho_1, _size_rho_0_0, _size_rho_0_1, _size_rho_0_2, _stride_rho_0_0, _stride_rho_0_1, _stride_rho_0_2, _stride_rho_1_0, _stride_rho_1_1, _stride_rho_1_2, order_0, order_1, rate_coefficient, stoech_0, stoech_1);
static FUNC_PREFIX void reactionkernelbulk_2_double_precision_reactionkernelbulk_2_double_precision(double *RESTRICT _data_rho_0, double *RESTRICT _data_rho_1, int64_t const _size_rho_0_0, int64_t const _size_rho_0_1, int64_t const _size_rho_0_2, int64_t const _stride_rho_0_0, int64_t const _stride_rho_0_1, int64_t const _stride_rho_0_2, int64_t const _stride_rho_1_0, int64_t const _stride_rho_1_1, int64_t const _stride_rho_1_2, double order_0, double order_1, double rate_coefficient, double stoech_0, double stoech_1)