64 auto rho_0 =
block->getData<field::GhostLayerField<double, 1>>(
rho_0ID);
69 WALBERLA_ASSERT_GREATER_EQUAL(0, -int_c(rho_0->nrOfGhostLayers()))
70 double *
RESTRICT _data_rho_0 = rho_0->dataAt(0, 0, 0, 0);
71 WALBERLA_ASSERT_GREATER_EQUAL(rho_0->xSizeWithGhostLayer(), int64_t(int64_c(rho_0->xSize()) + 0))
72 const int64_t _size_rho_0_0 = int64_t(int64_c(rho_0->xSize()) + 0);
73 WALBERLA_ASSERT_GREATER_EQUAL(rho_0->ySizeWithGhostLayer(), int64_t(int64_c(rho_0->ySize()) + 0))
74 const int64_t _size_rho_0_1 = int64_t(int64_c(rho_0->ySize()) + 0);
75 WALBERLA_ASSERT_GREATER_EQUAL(rho_0->zSizeWithGhostLayer(), int64_t(int64_c(rho_0->zSize()) + 0))
76 const int64_t _size_rho_0_2 = int64_t(int64_c(rho_0->zSize()) + 0);
77 const int64_t _stride_rho_0_0 = int64_t(rho_0->xStride());
78 const int64_t _stride_rho_0_1 = int64_t(rho_0->yStride());
79 const int64_t _stride_rho_0_2 = int64_t(rho_0->zStride());
85 CellInterval ci = globalCellInterval;
86 CellInterval blockBB = blocks->getBlockCellBB(*
block);
87 blockBB.expand(ghostLayers);
88 ci.intersect(blockBB);
89 blocks->transformGlobalToBlockLocalCellInterval(ci, *
block);
93 auto rho_0 =
block->getData<field::GhostLayerField<double, 1>>(
rho_0ID);
98 WALBERLA_ASSERT_GREATER_EQUAL(ci.xMin(), -int_c(rho_0->nrOfGhostLayers()))
99 WALBERLA_ASSERT_GREATER_EQUAL(ci.yMin(), -int_c(rho_0->nrOfGhostLayers()))
100 WALBERLA_ASSERT_GREATER_EQUAL(ci.zMin(), -int_c(rho_0->nrOfGhostLayers()))
101 double *
RESTRICT _data_rho_0 = rho_0->dataAt(ci.xMin(), ci.yMin(), ci.zMin(), 0);
102 WALBERLA_ASSERT_GREATER_EQUAL(rho_0->xSizeWithGhostLayer(), int64_t(int64_c(ci.xSize()) + 0))
103 const int64_t _size_rho_0_0 = int64_t(int64_c(ci.xSize()) + 0);
104 WALBERLA_ASSERT_GREATER_EQUAL(rho_0->ySizeWithGhostLayer(), int64_t(int64_c(ci.ySize()) + 0))
105 const int64_t _size_rho_0_1 = int64_t(int64_c(ci.ySize()) + 0);
106 WALBERLA_ASSERT_GREATER_EQUAL(rho_0->zSizeWithGhostLayer(), int64_t(int64_c(ci.zSize()) + 0))
107 const int64_t _size_rho_0_2 = int64_t(int64_c(ci.zSize()) + 0);
108 const int64_t _stride_rho_0_0 = int64_t(rho_0->xStride());
109 const int64_t _stride_rho_0_1 = int64_t(rho_0->yStride());
110 const int64_t _stride_rho_0_2 = int64_t(rho_0->zStride());
static FUNC_PREFIX void reactionkernelbulk_1_double_precision_reactionkernelbulk_1_double_precision(double *RESTRICT _data_rho_0, 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, double order_0, double rate_coefficient, double stoech_0)