olb::equilibria::Chopard::type< DESCRIPTOR, MOMENTA > Struct Template Reference

#include <equilibrium.h>

Collaboration diagram for olb::equilibria::Chopard::type< DESCRIPTOR, MOMENTA >:

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

Public Member Functions
template<typename CELL , typename RHO , typename U , typename FEQ , typename VS2 , typename V = typename CELL::value_t>
CellStatistic< V >	compute (CELL &cell, RHO &rho, U &u, FEQ &fEq, VS2 vs2=V{1}/descriptors::invCs2< V, DESCRIPTOR >()) any_platform
template<typename CELL , typename PARAMETERS , typename FEQ , typename V = typename CELL::value_t>
CellStatistic< V >	compute (CELL &cell, PARAMETERS &parameters, FEQ &fEq) any_platform

Detailed Description

template<typename DESCRIPTOR, typename MOMENTA>
struct olb::equilibria::Chopard::type< DESCRIPTOR, MOMENTA >

Definition at line 306 of file equilibrium.h.

Member Typedef Documentation

MomentaF

template<typename DESCRIPTOR , typename MOMENTA >

using olb::equilibria::Chopard::type< DESCRIPTOR, MOMENTA >::MomentaF = typename MOMENTA::template type<DESCRIPTOR>

Definition at line 307 of file equilibrium.h.

Member Function Documentation

compute() [1/2]

template<typename DESCRIPTOR , typename MOMENTA >

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

CellStatistic< V > olb::equilibria::Chopard::type< DESCRIPTOR, MOMENTA >::compute ( CELL & cell, PARAMETERS & parameters, FEQ & fEq )

inline

Definition at line 334 of file equilibrium.h.

                                                                                        {
      V rho, u[DESCRIPTOR::d];
      MomentaF().computeRhoU(cell, rho, u);
      const V vs2 = parameters.template get<SPEED_OF_SOUND>();
      const V uSqr = util::normSqr<V,DESCRIPTOR::d>(u);
      for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
        if (iPop==0) {
          return rho * (V{1} - vs2 * descriptors::invCs2<V,DESCRIPTOR>()*(V{1}-descriptors::t<V,DESCRIPTOR>(0))
                              - descriptors::t<V,DESCRIPTOR>(0)/V{2}*descriptors::invCs2<V,DESCRIPTOR>()*uSqr)
                  - descriptors::t<V,DESCRIPTOR>(0);
        }
        else {
          V c_u{};
          for (int iD=0; iD < DESCRIPTOR::d; ++iD) {
            c_u += descriptors::c<DESCRIPTOR>(iPop,iD)*u[iD];
          }
          return rho * descriptors::t<V,DESCRIPTOR>(iPop) * descriptors::invCs2<V,DESCRIPTOR>()
                      * (vs2 + c_u
                        + descriptors::invCs2<V,DESCRIPTOR>() / V{2} * c_u*c_u
                        - uSqr / V{2})
            - descriptors::t<V,DESCRIPTOR>(iPop);
        }
      }
      return {rho, uSqr};
    };

References olb::descriptors::c(), olb::descriptors::invCs2(), olb::util::normSqr(), and olb::descriptors::t().

Here is the call graph for this function:

compute() [2/2]

template<typename DESCRIPTOR , typename MOMENTA >

template<typename CELL , typename RHO , typename U , typename FEQ , typename VS2 , typename V = typename CELL::value_t>

CellStatistic< V > olb::equilibria::Chopard::type< DESCRIPTOR, MOMENTA >::compute ( CELL & cell, RHO & rho, U & u, FEQ & fEq, VS2 vs2 = V{1} / descriptors::invCs2<V,DESCRIPTOR>() )

inline

Definition at line 310 of file equilibrium.h.

                                                                              {1} / descriptors::invCs2<V,DESCRIPTOR>()) any_platform {
      const V uSqr = util::normSqr<V,DESCRIPTOR::d>(u);
      for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
        if (iPop==0) {
          return rho * (V{1} - vs2 * descriptors::invCs2<V,DESCRIPTOR>()*(V{1}-descriptors::t<V,DESCRIPTOR>(0))
                              - descriptors::t<V,DESCRIPTOR>(0)/V{2}*descriptors::invCs2<V,DESCRIPTOR>()*uSqr)
                  - descriptors::t<V,DESCRIPTOR>(0);
        }
        else {
          V c_u{};
          for (int iD=0; iD < DESCRIPTOR::d; ++iD) {
            c_u += descriptors::c<DESCRIPTOR>(iPop,iD)*u[iD];
          }
          return rho * descriptors::t<V,DESCRIPTOR>(iPop) * descriptors::invCs2<V,DESCRIPTOR>()
                      * (vs2 + c_u
                        + descriptors::invCs2<V,DESCRIPTOR>() / V{2} * c_u*c_u
                        - uSqr / V{2})
            - descriptors::t<V,DESCRIPTOR>(iPop);
        }
      }
      return {rho, uSqr};
    };

---

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

• src/dynamics/equilibrium.h