olb::opti::dualLbDynamicsHelpers< T, DESCRIPTOR > Struct Template Reference

All helper functions are inside this structure. More...

#include <dualLbHelpers.h>

Collaboration diagram for olb::opti::dualLbDynamicsHelpers< T, DESCRIPTOR >:

Static Public Member Functions
static T	equilibrium (int iPop, int jPop, T rho, const T u[DESCRIPTOR::d])
	Computation of adjoint equilibrium distribution operator dEq_i/dF_j -> Mathias J. Krause.
static T	equilibrium2 (int iPop, int jPop, T rho, const T u[DESCRIPTOR::d])
	Computation of adjoint equilibrium distribution operator dEq_i/dF_j (discrete approach) -> Geng Liu.
static void	computeRhoU (T const *cell, T &rho, T u[DESCRIPTOR::d])
	Computation of hydrodynamic variables.
static void	computeRhoJ (T const *cell, T &rho, T j[DESCRIPTOR::d])
	Computation of hydrodynamic variables.

Detailed Description

template<typename T, typename DESCRIPTOR>
struct olb::opti::dualLbDynamicsHelpers< T, DESCRIPTOR >

All helper functions are inside this structure.

Definition at line 148 of file dualLbHelpers.h.

Member Function Documentation

computeRhoJ()

template<typename T , typename DESCRIPTOR >

static void olb::opti::dualLbDynamicsHelpers< T, DESCRIPTOR >::computeRhoJ ( T const * cell, T & rho, T j[DESCRIPTOR::d] )

inlinestatic

Computation of hydrodynamic variables.

Definition at line 276 of file dualLbHelpers.h.

  {
 
    lbm<DESCRIPTOR>::computeRhoJ(&cell[0], rho, j);
    rho-=T(1);
  }

References olb::lbm< DESCRIPTOR >::computeRhoJ().

Here is the call graph for this function:

computeRhoU()

template<typename T , typename DESCRIPTOR >

static void olb::opti::dualLbDynamicsHelpers< T, DESCRIPTOR >::computeRhoU ( T const * cell, T & rho, T u[DESCRIPTOR::d] )

inlinestatic

Computation of hydrodynamic variables.

Definition at line 257 of file dualLbHelpers.h.

  {
 
    rho = T();
    for (int iD=0; iD < DESCRIPTOR::d; ++iD) {
      u[iD] = T();
    }
    for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
      rho += cell[iPop];
      for (int iD=0; iD < DESCRIPTOR::d; ++iD) {
        u[iD] += cell[iPop]*descriptors::c<DESCRIPTOR>(iPop,iD);
      }
    }
    for (int iD=0; iD < DESCRIPTOR::d; ++iD) {
      u[iD] /= rho;
    }
  }

Here is the caller graph for this function:

equilibrium()

template<typename T , typename DESCRIPTOR >

static T olb::opti::dualLbDynamicsHelpers< T, DESCRIPTOR >::equilibrium ( int iPop, int jPop, T rho, const T u[DESCRIPTOR::d] )

inlinestatic

Computation of adjoint equilibrium distribution operator dEq_i/dF_j -> Mathias J. Krause.

Definition at line 150 of file dualLbHelpers.h.

  {
    T eq = olb::equilibrium<DESCRIPTOR>::template secondOrder(iPop, rho, u) + descriptors::t<T,DESCRIPTOR>(iPop);
    T sum = T();
    for (int iD=0; iD < DESCRIPTOR::d; ++iD) {
      sum += (u[iD] - descriptors::c<DESCRIPTOR>(jPop,iD))*(u[iD] - descriptors::c<DESCRIPTOR>(iPop,iD));
    } //std::cout<<u[0]<<std::endl;
    return (descriptors::invCs2<T,DESCRIPTOR>()*sum + T(1))/rho*eq;
  }

equilibrium2()

template<typename T , typename DESCRIPTOR >

static T olb::opti::dualLbDynamicsHelpers< T, DESCRIPTOR >::equilibrium2 ( int iPop, int jPop, T rho, const T u[DESCRIPTOR::d] )

inlinestatic

Computation of adjoint equilibrium distribution operator dEq_i/dF_j (discrete approach) -> Geng Liu.

Definition at line 161 of file dualLbHelpers.h.

  {
    T u_c = T();
    for (int iD=0; iD < DESCRIPTOR::d; ++iD) {
      u_c += u[iD]*descriptors::c<DESCRIPTOR>(iPop,iD);
    }
    T sum = T();
    for (int iD=0; iD < DESCRIPTOR::d; ++iD) {
      sum += (u[iD]*u[iD] + 2.*descriptors::c<DESCRIPTOR>(jPop,iD)*(descriptors::c<DESCRIPTOR>(iPop,iD)-u[iD]) )
            * descriptors::invCs2<T,DESCRIPTOR>()*0.5
         + (u_c*descriptors::c<DESCRIPTOR>(iPop,iD)*(2.*descriptors::c<DESCRIPTOR>(jPop,iD)-u[iD]))
            * descriptors::invCs2<T,DESCRIPTOR>()*descriptors::invCs2<T,DESCRIPTOR>()*0.5;
    }
    return descriptors::t<T,DESCRIPTOR>(iPop)*(1.+sum);
  }

---

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

• src/optimization/dualLbHelpers.h