28#ifndef FD_SCHEMES_CENTRAL_HH
29#define FD_SCHEMES_CENTRAL_HH
36template <
unsigned D,
typename T>
37template <
typename PARAMETERS>
39apply(T& f0, T f[], T F[], T u[], PARAMETERS& params)
42 for (
unsigned iD=0; iD<D; ++iD) {
43 fNew += u[iD]*(F[getArrayPos<1>(0,iD)] - f[getArrayPos<1>(0,iD)]);
48template <
unsigned D,
typename T>
49template <
typename PARAMETERS>
51apply(T& f0, T f[], T F[], T u[], PARAMETERS& params)
54 for (
unsigned iD=0; iD<D; ++iD) {
55 fNew += F[getArrayPos<1>(0,iD)] + f[getArrayPos<1>(0,iD)];
58 return fNew * params.template get<fdParams::Diffusivity>();
61template <
unsigned D,
typename T>
67template <
unsigned D,
typename T>
68template <
typename PARAMETERS>
70apply(T fOut[], T& f0, T fIn[],
int normal[], T u[], PARAMETERS& params)
72 for (
unsigned iD=0; iD<D; ++iD) {
73 fOut[getArrayPos<1>(0,iD)] = fIn[getArrayPos<1>(0,iD)];
Top level namespace for all of OpenLB.
static constexpr int getExtraExtent()
static void apply(T fOut[], T &f0, T fIn[], int normal[], T u[], PARAMETERS ¶ms)
static T apply(T &f0, T f[], T F[], T u[], PARAMETERS ¶ms)
static T apply(T &f0, T f[], T F[], T u[], PARAMETERS ¶ms)