ESPResSo
Extensible Simulation Package for Research on Soft Matter Systems
Loading...
Searching...
No Matches
ReactionKernelBulk_5_double_precision.h
Go to the documentation of this file.
1//======================================================================================================================
2//
3// This file is part of waLBerla. waLBerla is free software: you can
4// redistribute it and/or modify it under the terms of the GNU General Public
5// License as published by the Free Software Foundation, either version 3 of
6// the License, or (at your option) any later version.
7//
8// waLBerla is distributed in the hope that it will be useful, but WITHOUT
9// ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
10// FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
11// for more details.
12//
13// You should have received a copy of the GNU General Public License along
14// with waLBerla (see COPYING.txt). If not, see <http://www.gnu.org/licenses/>.
15//
16//! \\file ReactionKernelBulk_5_double_precision.h
17//! \\author pystencils
18//======================================================================================================================
19
20// kernel generated with pystencils v1.3.3, lbmpy v1.3.3,
21// lbmpy_walberla/pystencils_walberla from waLBerla commit
22// b0842e1a493ce19ef1bbb8d2cf382fc343970a7f
23
24#pragma once
25#include "core/DataTypes.h"
26#include "core/logging/Logging.h"
27
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"
33#include <set>
34
35#ifdef __GNUC__
36#define RESTRICT __restrict__
37#elif _MSC_VER
38#define RESTRICT __restrict
39#else
40#define RESTRICT
41#endif
42
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"
48#endif
49
50namespace walberla {
51namespace pystencils {
52
54public:
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)
61 : rho_0ID(rho_0ID_), rho_1ID(rho_1ID_), rho_2ID(rho_2ID_),
62 rho_3ID(rho_3ID_), rho_4ID(rho_4ID_), order_0_(order_0),
63 order_1_(order_1), order_2_(order_2), order_3_(order_3),
64 order_4_(order_4), rate_coefficient_(rate_coefficient),
65 stoech_0_(stoech_0), stoech_1_(stoech_1), stoech_2_(stoech_2),
66 stoech_3_(stoech_3), stoech_4_(stoech_4){};
67
68 void run(IBlock *block);
69
70 void runOnCellInterval(const shared_ptr<StructuredBlockStorage> &blocks,
71 const CellInterval &globalCellInterval,
72 cell_idx_t ghostLayers, IBlock *block);
73
74 void operator()(IBlock *block) { run(block); }
75
76 static std::function<void(IBlock *)>
77 getSweep(const shared_ptr<ReactionKernelBulk_5_double_precision> &kernel) {
78 return [kernel](IBlock *b) { kernel->run(b); };
79 }
80
81 static std::function<void(IBlock *)> getSweepOnCellInterval(
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);
87 };
88 }
89
90 std::function<void(IBlock *)> getSweep() {
91 return [this](IBlock *b) { this->run(b); };
92 }
93
94 std::function<void(IBlock *)>
95 getSweepOnCellInterval(const shared_ptr<StructuredBlockStorage> &blocks,
96 const CellInterval &globalCellInterval,
97 cell_idx_t ghostLayers = 1) {
98 return [this, blocks, globalCellInterval, ghostLayers](IBlock *b) {
99 this->runOnCellInterval(blocks, globalCellInterval, ghostLayers, b);
100 };
101 }
102
103 void configure(const shared_ptr<StructuredBlockStorage> &blocks,
104 IBlock *block) {}
105
106 BlockDataID rho_0ID;
107 BlockDataID rho_1ID;
108 BlockDataID rho_2ID;
109 BlockDataID rho_3ID;
110 BlockDataID rho_4ID;
111 double order_0_;
112 double order_1_;
113 double order_2_;
114 double order_3_;
115 double order_4_;
117 double stoech_0_;
118 double stoech_1_;
119 double stoech_2_;
120 double stoech_3_;
121 double stoech_4_;
122};
123
124} // namespace pystencils
125} // namespace walberla
126
127#if (defined WALBERLA_CXX_COMPILER_IS_GNU) || \
128 (defined WALBERLA_CXX_COMPILER_IS_CLANG)
129#pragma GCC diagnostic pop
130#endif
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)
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)
Definition elc.cpp:172
\file PackInfoPdfDoublePrecision.cpp \author pystencils