25#include "core/DataTypes.h"
26#include "core/logging/Logging.h"
28#include "domain_decomposition/BlockDataID.h"
29#include "domain_decomposition/IBlock.h"
30#include "domain_decomposition/StructuredBlockStorage.h"
31#include "field/GhostLayerField.h"
32#include "field/SwapableCompare.h"
36#define RESTRICT __restrict__
38#define RESTRICT __restrict
43#if (defined WALBERLA_CXX_COMPILER_IS_GNU) || \
44 (defined WALBERLA_CXX_COMPILER_IS_CLANG)
45#pragma GCC diagnostic push
46#pragma GCC diagnostic ignored "-Wunused-parameter"
47#pragma GCC diagnostic ignored "-Wreorder"
56 BlockDataID rho_0ID_, BlockDataID rho_1ID_, BlockDataID rho_2ID_,
57 BlockDataID rho_3ID_, BlockDataID rho_4ID_,
double order_0,
58 double order_1,
double order_2,
double order_3,
double order_4,
59 double rate_coefficient,
double stoech_0,
double stoech_1,
60 double stoech_2,
double stoech_3,
double stoech_4)
71 const CellInterval &globalCellInterval,
72 cell_idx_t ghostLayers, IBlock *
block);
76 static std::function<void(IBlock *)>
77 getSweep(
const shared_ptr<ReactionKernelBulk_5_double_precision> &kernel) {
78 return [kernel](IBlock *b) { kernel->run(b); };
82 const shared_ptr<ReactionKernelBulk_5_double_precision> &kernel,
83 const shared_ptr<StructuredBlockStorage> &blocks,
84 const CellInterval &globalCellInterval, cell_idx_t ghostLayers = 1) {
85 return [kernel, blocks, globalCellInterval, ghostLayers](IBlock *b) {
86 kernel->runOnCellInterval(blocks, globalCellInterval, ghostLayers, b);
91 return [
this](IBlock *b) { this->
run(b); };
94 std::function<void(IBlock *)>
96 const CellInterval &globalCellInterval,
97 cell_idx_t ghostLayers = 1) {
98 return [
this, blocks, globalCellInterval, ghostLayers](IBlock *b) {
103 void configure(
const shared_ptr<StructuredBlockStorage> &blocks,
127#if (defined WALBERLA_CXX_COMPILER_IS_GNU) || \
128 (defined WALBERLA_CXX_COMPILER_IS_CLANG)
129#pragma GCC diagnostic pop
void operator()(IBlock *block)
static std::function< void(IBlock *)> getSweep(const shared_ptr< ReactionKernelBulk_5_double_precision > &kernel)
ReactionKernelBulk_5_double_precision(BlockDataID rho_0ID_, BlockDataID rho_1ID_, BlockDataID rho_2ID_, BlockDataID rho_3ID_, BlockDataID rho_4ID_, double order_0, double order_1, double order_2, double order_3, double order_4, double rate_coefficient, double stoech_0, double stoech_1, double stoech_2, double stoech_3, double stoech_4)
std::function< void(IBlock *)> getSweep()
void configure(const shared_ptr< StructuredBlockStorage > &blocks, IBlock *block)
static std::function< void(IBlock *)> getSweepOnCellInterval(const shared_ptr< ReactionKernelBulk_5_double_precision > &kernel, const shared_ptr< StructuredBlockStorage > &blocks, const CellInterval &globalCellInterval, cell_idx_t ghostLayers=1)
std::function< void(IBlock *)> getSweepOnCellInterval(const shared_ptr< StructuredBlockStorage > &blocks, const CellInterval &globalCellInterval, cell_idx_t ghostLayers=1)
void runOnCellInterval(const shared_ptr< StructuredBlockStorage > &blocks, const CellInterval &globalCellInterval, cell_idx_t ghostLayers, IBlock *block)
static double * block(double *p, std::size_t index, std::size_t size)
\file PackInfoPdfDoublePrecision.cpp \author pystencils