OpenLB 1.7
Loading...
Searching...
No Matches
fdPostProcessor3D.hh
Go to the documentation of this file.
1/* This file is part of the OpenLB library
2 *
3 * Copyright (C) 2020 Davide Dapelo
4 * E-mail contact: info@openlb.net
5 * The most recent release of OpenLB can be downloaded at
6 * <http://www.openlb.net/>
7 *
8 * This program is free software; you can redistribute it and/or
9 * modify it under the terms of the GNU General Public License
10 * as published by the Free Software Foundation; either version 2
11 * of the License, or (at your option) any later version.
12 *
13 * This program is distributed in the hope that it will be useful,
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16 * GNU General Public License for more details.
17 *
18 * You should have received a copy of the GNU General Public
19 * License along with this program; if not, write to the Free
20 * Software Foundation, Inc., 51 Franklin Street, Fifth Floor,
21 * Boston, MA 02110-1301, USA.
22*/
23
28#ifndef AD_POST_PROCESSOR_3D_DEV03_HH
29#define AD_POST_PROCESSOR_3D_DEV03_HH
30
31#include <string>
32namespace olb {
33
34
36
37template<typename T, typename DESCRIPTOR, typename FIELD, typename SOURCE>
39{
40 static_assert(DESCRIPTOR::template size<FIELD>() == 2, "FIELD must have size 2." );
41}
43template<typename T, typename DESCRIPTOR, typename FIELD, typename SOURCE>
44template <typename CELL>
46{
47 if constexpr (! std::is_void<SOURCE>::value) {
48 *fNew += cell.template getFieldPointer<SOURCE>()[0];
49 }
50}
51
52
54
55template<typename T, typename DESCRIPTOR, typename MODEL, typename PARAMS, typename FIELD, typename SOURCE>
60
61template<typename T, typename DESCRIPTOR, typename MODEL, typename PARAMS, typename FIELD, typename SOURCE>
66
67template<typename T, typename DESCRIPTOR, typename MODEL, typename PARAMS, typename FIELD, typename SOURCE>
68template <typename CELL, typename PARAMETERS>
70{
71 std::size_t iT = vars.template get<fd::fdParams::Timestep>();
72 T f[MODEL::extent()*DESCRIPTOR::d], F[MODEL::extent()*DESCRIPTOR::d];
73 T u[DESCRIPTOR::d];
74 T* fNew = fd::accessNew<T,FIELD>(cell, iT);
75 T* f0 = fd::accessOld<T,FIELD>(cell, iT);
76 for (int iN=1; iN<=MODEL::extent(); ++iN) {
77 f[fd::getArrayPos<MODEL::extent()>(iN-1,0)] = *fd::accessOld<T,FIELD>( cell.neighbor({-iN, 0, 0}), iT );
78 F[fd::getArrayPos<MODEL::extent()>(iN-1,0)] = *fd::accessOld<T,FIELD>( cell.neighbor({+iN, 0, 0}), iT );
79 f[fd::getArrayPos<MODEL::extent()>(iN-1,1)] = *fd::accessOld<T,FIELD>( cell.neighbor({ 0, -iN, 0}), iT );
80 F[fd::getArrayPos<MODEL::extent()>(iN-1,1)] = *fd::accessOld<T,FIELD>( cell.neighbor({ 0, +iN, 0}), iT );
81 f[fd::getArrayPos<MODEL::extent()>(iN-1,2)] = *fd::accessOld<T,FIELD>( cell.neighbor({ 0, 0, -iN}), iT );
82 F[fd::getArrayPos<MODEL::extent()>(iN-1,2)] = *fd::accessOld<T,FIELD>( cell.neighbor({ 0, 0, +iN}), iT );
83 }
84 cell.computeU(u);
85 MODEL::apply(fNew, f0, f, F, u, vars);
86 this->applySourceTerm(fNew, cell);
87}
88
89} // namespace olb
90
91#endif
void applySourceTerm(T *fNew, CELL &cell) any_platform
void apply(CELL &cell, PARAMETERS &vars) any_platform
Top level namespace for all of OpenLB.