olb::forcing::LaddVerberg::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION > Struct Template Reference

#include <collisionMRT.h>

Collaboration diagram for olb::forcing::LaddVerberg::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >:

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

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

Detailed Description

template<typename DESCRIPTOR, typename MOMENTA, typename EQUILIBRIUM, typename COLLISION>
struct olb::forcing::LaddVerberg::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >

Definition at line 101 of file collisionMRT.h.

Member Typedef Documentation

CollisionO

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM , typename COLLISION >

using olb::forcing::LaddVerberg::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >::CollisionO = typename COLLISION::template type<DESCRIPTOR,MOMENTA,EQUILIBRIUM>

Definition at line 103 of file collisionMRT.h.

MomentaF

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM , typename COLLISION >

using olb::forcing::LaddVerberg::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >::MomentaF = typename MOMENTA::template type<DESCRIPTOR>

Definition at line 102 of file collisionMRT.h.

Member Function Documentation

apply()

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM , typename COLLISION >

template<typename CELL , typename PARAMETERS , typename V = typename CELL::value_t>

CellStatistic< V > olb::forcing::LaddVerberg::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >::apply ( CELL & cell, PARAMETERS & parameters )

inline

Definition at line 106 of file collisionMRT.h.

                                                               {
      V rho, u[DESCRIPTOR::d];
      MomentaF().computeRhoU(cell, rho, u);
      const auto statistic = CollisionO().apply(cell, parameters);
      const V omega = parameters.template get<descriptors::OMEGA>();
      // While this duplication can be resolved using CSE it should be extracted into a helper
      V rt[DESCRIPTOR::q] { }; // relaxation times vector.
      for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
        rt[iPop] = descriptors::s<V,DESCRIPTOR>(iPop);
      }
      for (int iPop=0; iPop < descriptors::shearIndexes<DESCRIPTOR>(); ++iPop) {
        rt[descriptors::shearViscIndexes<DESCRIPTOR>(iPop)] = omega;
      }
      V invM_S[DESCRIPTOR::q][DESCRIPTOR::q]; // relaxation times matrix
      for (int iPop = 0; iPop < DESCRIPTOR::q; ++iPop) {
        for (int jPop = 0; jPop < DESCRIPTOR::q; ++jPop) {
          invM_S[iPop][jPop] = V{};
          for (int kPop = 0; kPop < DESCRIPTOR::q; ++kPop) {
            if (kPop == jPop) {
              invM_S[iPop][jPop] += descriptors::invM<V,DESCRIPTOR>(iPop,kPop) * rt[kPop];
            }
          }
        }
      }
      const auto force = cell.template getField<descriptors::FORCE>();
      mrt<DESCRIPTOR>::addExternalForce(cell, rho, u, invM_S, force);
      return statistic;
    };

References olb::mrt< DESCRIPTOR >::addExternalForce(), olb::descriptors::invM(), olb::descriptors::s(), and olb::descriptors::shearViscIndexes().

Here is the call graph for this function:

---

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

• src/dynamics/collisionMRT.h