38apply(T& f0, T f[], T F[], T u[], PARAMETERS& params)
42 for (
unsigned iD=0; iD<D; ++iD) {
43 fNew += F[getArrayPos<1>(0,iD)] + f[getArrayPos<1>(0,iD)];
44 fCorr += (F[getArrayPos<1>(0,iD)] + f[getArrayPos<1>(0,iD)] - 2.*f0) *
abs(u[iD]);
47 return fNew * params.template get<fdParams::Diffusivity>()
48 - 0.5 * fCorr * params.template get<fdParams::AntiDiffusivityTuning>();