ESPResSo
Extensible Simulation Package for Research on Soft Matter Systems
Loading...
Searching...
No Matches
DiffusiveFluxKernelWithElectrostaticThermalized_single_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 DiffusiveFluxKernelWithElectrostaticThermalized_single_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 jID_, BlockDataID phiID_, BlockDataID rhoID_, float D,
57 float f_ext_0, float f_ext_1, float f_ext_2, uint32_t field_size_0,
58 uint32_t field_size_1, uint32_t field_size_2, float kT, uint32_t seed,
59 uint32_t time_step, float z)
60 : jID(jID_), phiID(phiID_), rhoID(rhoID_), D_(D), f_ext_0_(f_ext_0),
61 f_ext_1_(f_ext_1), f_ext_2_(f_ext_2), field_size_0_(field_size_0),
62 field_size_1_(field_size_1), field_size_2_(field_size_2), kT_(kT),
63 seed_(seed), time_step_(time_step), z_(z), configured_(false){};
64
65 void run(IBlock *block);
66
67 void runOnCellInterval(const shared_ptr<StructuredBlockStorage> &blocks,
68 const CellInterval &globalCellInterval,
69 cell_idx_t ghostLayers, IBlock *block);
70
71 void operator()(IBlock *block) { run(block); }
72
73 static std::function<void(IBlock *)>
74 getSweep(const shared_ptr<
76 &kernel) {
77 return [kernel](IBlock *b) { kernel->run(b); };
78 }
79
80 static std::function<void(IBlock *)> getSweepOnCellInterval(
81 const shared_ptr<
83 &kernel,
84 const shared_ptr<StructuredBlockStorage> &blocks,
85 const CellInterval &globalCellInterval, cell_idx_t ghostLayers = 1) {
86 return [kernel, blocks, globalCellInterval, ghostLayers](IBlock *b) {
87 kernel->runOnCellInterval(blocks, globalCellInterval, ghostLayers, b);
88 };
89 }
90
91 std::function<void(IBlock *)> getSweep() {
92 return [this](IBlock *b) { this->run(b); };
93 }
94
95 std::function<void(IBlock *)>
96 getSweepOnCellInterval(const shared_ptr<StructuredBlockStorage> &blocks,
97 const CellInterval &globalCellInterval,
98 cell_idx_t ghostLayers = 1) {
99 return [this, blocks, globalCellInterval, ghostLayers](IBlock *b) {
100 this->runOnCellInterval(blocks, globalCellInterval, ghostLayers, b);
101 };
102 }
103
104 void configure(const shared_ptr<StructuredBlockStorage> &blocks,
105 IBlock *block) {
106 Cell BlockCellBB = blocks->getBlockCellBB(*block).min();
107 block_offset_0_ = uint32_t(BlockCellBB[0]);
108 block_offset_1_ = uint32_t(BlockCellBB[1]);
109 block_offset_2_ = uint32_t(BlockCellBB[2]);
110 configured_ = true;
111 }
112
113 BlockDataID jID;
114 BlockDataID phiID;
115 BlockDataID rhoID;
116 float D_;
120 float f_ext_0_;
121 float f_ext_1_;
122 float f_ext_2_;
126 float kT_;
127 uint32_t seed_;
128 uint32_t time_step_;
129 float z_;
131};
132
133} // namespace pystencils
134} // namespace walberla
135
136#if (defined WALBERLA_CXX_COMPILER_IS_GNU) || \
137 (defined WALBERLA_CXX_COMPILER_IS_CLANG)
138#pragma GCC diagnostic pop
139#endif
Definition Cell.hpp:97
std::function< void(IBlock *)> getSweepOnCellInterval(const shared_ptr< StructuredBlockStorage > &blocks, const CellInterval &globalCellInterval, cell_idx_t ghostLayers=1)
static std::function< void(IBlock *)> getSweepOnCellInterval(const shared_ptr< DiffusiveFluxKernelWithElectrostaticThermalized_single_precision > &kernel, const shared_ptr< StructuredBlockStorage > &blocks, const CellInterval &globalCellInterval, cell_idx_t ghostLayers=1)
DiffusiveFluxKernelWithElectrostaticThermalized_single_precision(BlockDataID jID_, BlockDataID phiID_, BlockDataID rhoID_, float D, float f_ext_0, float f_ext_1, float f_ext_2, uint32_t field_size_0, uint32_t field_size_1, uint32_t field_size_2, float kT, uint32_t seed, uint32_t time_step, float z)
void runOnCellInterval(const shared_ptr< StructuredBlockStorage > &blocks, const CellInterval &globalCellInterval, cell_idx_t ghostLayers, IBlock *block)
static std::function< void(IBlock *)> getSweep(const shared_ptr< DiffusiveFluxKernelWithElectrostaticThermalized_single_precision > &kernel)
static double * block(double *p, std::size_t index, std::size_t size)
Definition elc.cpp:172
\file PackInfoPdfDoublePrecision.cpp \author pystencils