373 {
374
375 using V = typename CELLS::template value_t<names::A>::value_t;
376
377 auto outletRho =
parameters.template get<RHO>();
378
379
380 auto& cellA = cells.template get<names::A>();
381
382
383 auto& cellB = cells.template get<names::B>();
384
385 V rho0, phi;
386 rho0 = cellA.computeRho();
387 phi = cellB.computeRho();
388
389 cellA.defineRho(outletRho);
390
391 V temp = phi * outletRho / rho0;
392 cellB.defineRho(temp);
393
394
395
396 if constexpr (CELLS::map_t::keys_t::template contains<names::C>()){
397 auto& cellC = cells.template get<names::C>();
398
399 V psi = cellC.computeRho();
400 temp = psi * outletRho / rho0;
401 cellC.defineRho(temp);
402 }
403
404 }
meta::list< RHO > parameters