olb::collision::TRT::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM > Struct Template Reference

#include <collision.h>

Collaboration diagram for olb::collision::TRT::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >:

Public Types
using	MomentaF = typename MOMENTA::template type<DESCRIPTOR>
using	EquilibriumF = typename EQUILIBRIUM::template type<DESCRIPTOR,MOMENTA>

Public Member Functions
template<CONCEPT(MinimalCell) CELL, typename PARAMETERS , typename V = typename CELL::value_t>
CellStatistic< V >	apply (CELL &cell, PARAMETERS &parameters) any_platform

Detailed Description

template<typename DESCRIPTOR, typename MOMENTA, typename EQUILIBRIUM>
struct olb::collision::TRT::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >

Definition at line 215 of file collision.h.

Member Typedef Documentation

EquilibriumF

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

using olb::collision::TRT::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >::EquilibriumF = typename EQUILIBRIUM::template type<DESCRIPTOR,MOMENTA>

Definition at line 217 of file collision.h.

MomentaF

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

using olb::collision::TRT::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >::MomentaF = typename MOMENTA::template type<DESCRIPTOR>

Definition at line 216 of file collision.h.

Member Function Documentation

apply()

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

template<CONCEPT(MinimalCell) CELL, typename PARAMETERS , typename V = typename CELL::value_t>

CellStatistic< V > olb::collision::TRT::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >::apply ( CELL & cell, PARAMETERS & parameters )

inline

Definition at line 220 of file collision.h.

                                                                            {
      V fPlus[DESCRIPTOR::q], fMinus[DESCRIPTOR::q];
      V fEq[DESCRIPTOR::q], fEqPlus[DESCRIPTOR::q], fEqMinus[DESCRIPTOR::q];
      const auto statistic = EquilibriumF().compute(cell, parameters, fEq);
      for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
        fPlus[iPop]  = 0.5 * (cell[iPop] + cell[descriptors::opposite<DESCRIPTOR>(iPop)]);
        fMinus[iPop] = 0.5 * (cell[iPop] - cell[descriptors::opposite<DESCRIPTOR>(iPop)]);
      }
      for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
        fEqPlus[iPop]  = 0.5 * (fEq[iPop] + fEq[descriptors::opposite<DESCRIPTOR>(iPop)]);
        fEqMinus[iPop] = 0.5 * (fEq[iPop] - fEq[descriptors::opposite<DESCRIPTOR>(iPop)]);
      }
 
      const V omega = parameters.template get<descriptors::OMEGA>();
      const V magic = parameters.template get<MAGIC>();
      const V omega2 = V{1} / (magic / (V{1} / omega - 0.5) + V{0.5});
 
      for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
        cell[iPop] -= omega * (fPlus[iPop] - fEqPlus[iPop]) + omega2 * (fMinus[iPop] - fEqMinus[iPop]);
      }
      return statistic;
    };

---

The documentation for this struct was generated from the following file:

• src/dynamics/collision.h