38apply(T& f0, T f[], T F[], T u[], PARAMETERS& params)
41 for (
unsigned iD=0; iD<D; ++iD) {
42 fNew += u[iD] * ( u[iD]==0. ? T(0.)
43 : ( u[iD] >0. ? 3.*f0 - 4.*f[getArrayPos<2>(0,iD)] + f[getArrayPos<2>(1,iD)]
44 : -F[getArrayPos<2>(1,iD)] + 4.*f[getArrayPos<2>(0,iD)] - 3.*f0 ) );