222 {
223 using V = typename CELL::value_t;
224 using DESCRIPTOR = typename CELL::descriptor_t;
225 V w =
parameters.template get<descriptors::INTERFACE_WIDTH>();
226 V psi = cell.template getField<descriptors::PSI>();
227 auto phi = cell.template getField<descriptors::STATISTIC>();
228 if (
fabs(psi) >= 1.5*w) {
229 if (phi[0] > 0.5 && phi[0] < 0.995) {
230 for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
232 }
233 } else if (phi[0] < 0.5 && phi[0] > 0.005) {
234 for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
236 }
237 }
238 phi[0] = cell.computeRho();
239 cell.template setField<descriptors::STATISTIC>(phi);
240 }
241 }
constexpr T t(unsigned iPop, tag::CUM) any_platform
meta::list< descriptors::INTERFACE_WIDTH > parameters