62 auto rho =
block->getData<field::GhostLayerField<float, 1>>(
rhoID);
63 auto j =
block->getData<field::GhostLayerField<float, 13>>(
jID);
65 WALBERLA_ASSERT_GREATER_EQUAL(-1, -int_c(j->nrOfGhostLayers()))
66 float *
RESTRICT const _data_j = j->dataAt(-1, -1, -1, 0);
67 WALBERLA_ASSERT_GREATER_EQUAL(-1, -int_c(rho->nrOfGhostLayers()))
68 float *
RESTRICT _data_rho = rho->dataAt(-1, -1, -1, 0);
69 WALBERLA_ASSERT_GREATER_EQUAL(j->xSizeWithGhostLayer(), int64_t(int64_c(j->xSize()) + 2))
70 const int64_t _size_j_0 = int64_t(int64_c(j->xSize()) + 2);
71 WALBERLA_ASSERT_GREATER_EQUAL(j->ySizeWithGhostLayer(), int64_t(int64_c(j->ySize()) + 2))
72 const int64_t _size_j_1 = int64_t(int64_c(j->ySize()) + 2);
73 WALBERLA_ASSERT_GREATER_EQUAL(j->zSizeWithGhostLayer(), int64_t(int64_c(j->zSize()) + 2))
74 const int64_t _size_j_2 = int64_t(int64_c(j->zSize()) + 2);
75 const int64_t _stride_j_0 = int64_t(j->xStride());
76 const int64_t _stride_j_1 = int64_t(j->yStride());
77 const int64_t _stride_j_2 = int64_t(j->zStride());
78 const int64_t _stride_j_3 = int64_t(1 * int64_t(j->fStride()));
79 const int64_t _stride_rho_0 = int64_t(rho->xStride());
80 const int64_t _stride_rho_1 = int64_t(rho->yStride());
81 const int64_t _stride_rho_2 = int64_t(rho->zStride());
82 internal_990034b4e4dd57d2802b4bcb5f716e46::continuitykernel_single_precision_continuitykernel_single_precision(_data_j, _data_rho, _size_j_0, _size_j_1, _size_j_2, _stride_j_0, _stride_j_1, _stride_j_2, _stride_j_3, _stride_rho_0, _stride_rho_1, _stride_rho_2);
87 CellInterval ci = globalCellInterval;
88 CellInterval blockBB = blocks->getBlockCellBB(*
block);
89 blockBB.expand(ghostLayers);
90 ci.intersect(blockBB);
91 blocks->transformGlobalToBlockLocalCellInterval(ci, *
block);
95 auto rho =
block->getData<field::GhostLayerField<float, 1>>(
rhoID);
96 auto j =
block->getData<field::GhostLayerField<float, 13>>(
jID);
98 WALBERLA_ASSERT_GREATER_EQUAL(ci.xMin() - 1, -int_c(j->nrOfGhostLayers()))
99 WALBERLA_ASSERT_GREATER_EQUAL(ci.yMin() - 1, -int_c(j->nrOfGhostLayers()))
100 WALBERLA_ASSERT_GREATER_EQUAL(ci.zMin() - 1, -int_c(j->nrOfGhostLayers()))
101 float *
RESTRICT const _data_j = j->dataAt(ci.xMin() - 1, ci.yMin() - 1, ci.zMin() - 1, 0);
102 WALBERLA_ASSERT_GREATER_EQUAL(ci.xMin() - 1, -int_c(rho->nrOfGhostLayers()))
103 WALBERLA_ASSERT_GREATER_EQUAL(ci.yMin() - 1, -int_c(rho->nrOfGhostLayers()))
104 WALBERLA_ASSERT_GREATER_EQUAL(ci.zMin() - 1, -int_c(rho->nrOfGhostLayers()))
105 float *
RESTRICT _data_rho = rho->dataAt(ci.xMin() - 1, ci.yMin() - 1, ci.zMin() - 1, 0);
106 WALBERLA_ASSERT_GREATER_EQUAL(j->xSizeWithGhostLayer(), int64_t(int64_c(ci.xSize()) + 2))
107 const int64_t _size_j_0 = int64_t(int64_c(ci.xSize()) + 2);
108 WALBERLA_ASSERT_GREATER_EQUAL(j->ySizeWithGhostLayer(), int64_t(int64_c(ci.ySize()) + 2))
109 const int64_t _size_j_1 = int64_t(int64_c(ci.ySize()) + 2);
110 WALBERLA_ASSERT_GREATER_EQUAL(j->zSizeWithGhostLayer(), int64_t(int64_c(ci.zSize()) + 2))
111 const int64_t _size_j_2 = int64_t(int64_c(ci.zSize()) + 2);
112 const int64_t _stride_j_0 = int64_t(j->xStride());
113 const int64_t _stride_j_1 = int64_t(j->yStride());
114 const int64_t _stride_j_2 = int64_t(j->zStride());
115 const int64_t _stride_j_3 = int64_t(1 * int64_t(j->fStride()));
116 const int64_t _stride_rho_0 = int64_t(rho->xStride());
117 const int64_t _stride_rho_1 = int64_t(rho->yStride());
118 const int64_t _stride_rho_2 = int64_t(rho->zStride());
119 internal_990034b4e4dd57d2802b4bcb5f716e46::continuitykernel_single_precision_continuitykernel_single_precision(_data_j, _data_rho, _size_j_0, _size_j_1, _size_j_2, _stride_j_0, _stride_j_1, _stride_j_2, _stride_j_3, _stride_rho_0, _stride_rho_1, _stride_rho_2);
static FUNC_PREFIX void continuitykernel_single_precision_continuitykernel_single_precision(float *RESTRICT const _data_j, float *RESTRICT _data_rho, int64_t const _size_j_0, int64_t const _size_j_1, int64_t const _size_j_2, int64_t const _stride_j_0, int64_t const _stride_j_1, int64_t const _stride_j_2, int64_t const _stride_j_3, int64_t const _stride_rho_0, int64_t const _stride_rho_1, int64_t const _stride_rho_2)