olb::ForceCoupling3D Struct Reference

#include <freeEnergyCoupling3D.h>

Collaboration diagram for olb::ForceCoupling3D:

Public Member Functions
template<typename CELLS >
void	apply (CELLS &cells) any_platform

Static Public Attributes
static constexpr OperatorScope	scope = OperatorScope::PerCell

Detailed Description

Definition at line 114 of file freeEnergyCoupling3D.h.

Member Function Documentation

apply()

template<typename CELLS >

void olb::ForceCoupling3D::apply ( CELLS & cells )

inline

Definition at line 120 of file freeEnergyCoupling3D.h.

                                        {
 
    using V = typename CELLS::template value_t<names::A>::value_t;
    using DESCRIPTOR = typename CELLS::template value_t<names::A>::descriptor_t;
 
    // Get the cell of the first lattice
    auto& cellA = cells.template get<names::A>();
 
    // Get the cell of the second lattice
    auto& cellB = cells.template get<names::B>();
 
    V phi = cellA.computeRho();
    V rho = cellB.computeRho();
 
    // Calculatations for lattice A.
    auto cellA1 = cellA.neighbor(descriptors::c<DESCRIPTOR>(1));
    auto cellA2 = cellA.neighbor(descriptors::c<DESCRIPTOR>(2));
    auto cellA3 = cellA.neighbor(descriptors::c<DESCRIPTOR>(3));
    auto cellA4 = cellA.neighbor(descriptors::c<DESCRIPTOR>(4));
    auto cellA5 = cellA.neighbor(descriptors::c<DESCRIPTOR>(5));
    auto cellA6 = cellA.neighbor(descriptors::c<DESCRIPTOR>(6));
    auto cellA7 = cellA.neighbor(descriptors::c<DESCRIPTOR>(7));
    auto cellA8 = cellA.neighbor(descriptors::c<DESCRIPTOR>(8));
    auto cellA9 = cellA.neighbor(descriptors::c<DESCRIPTOR>(9));
    auto cellA10 = cellA.neighbor(descriptors::c<DESCRIPTOR>(10));
    auto cellA11 = cellA.neighbor(descriptors::c<DESCRIPTOR>(11));
    auto cellA12 = cellA.neighbor(descriptors::c<DESCRIPTOR>(12));
    auto cellA13 = cellA.neighbor(descriptors::c<DESCRIPTOR>(13));
    auto cellA14 = cellA.neighbor(descriptors::c<DESCRIPTOR>(14));
    auto cellA15 = cellA.neighbor(descriptors::c<DESCRIPTOR>(15));
    auto cellA16 = cellA.neighbor(descriptors::c<DESCRIPTOR>(16));
    auto cellA17 = cellA.neighbor(descriptors::c<DESCRIPTOR>(17));
    auto cellA18 = cellA.neighbor(descriptors::c<DESCRIPTOR>(18));
 
    V gradMuPhiX = 1./12. * (
    -     cellA6.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellA7.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellA4.template getField<descriptors::CHEM_POTENTIAL>()
    - 2 * cellA1.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellA5.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellA14.template getField<descriptors::CHEM_POTENTIAL>()
    + 2 * cellA10.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellA13.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellA16.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellA15.template getField<descriptors::CHEM_POTENTIAL>());
 
    V gradMuPhiY = 1./12. * (
    -     cellA4.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellA14.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellA8.template getField<descriptors::CHEM_POTENTIAL>()
    - 2 * cellA2.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellA9.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellA18.template getField<descriptors::CHEM_POTENTIAL>()
    + 2 * cellA11.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellA17.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellA5.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellA13.template getField<descriptors::CHEM_POTENTIAL>());
 
    V gradMuPhiZ = 1./12. * (
    -     cellA8.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellA18.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellA6.template getField<descriptors::CHEM_POTENTIAL>()
    - 2 * cellA3.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellA16.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellA7.template getField<descriptors::CHEM_POTENTIAL>()
    + 2 * cellA12.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellA15.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellA9.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellA17.template getField<descriptors::CHEM_POTENTIAL>());
 
 
    // Calculatations for lattice B.
    auto cellB1 = cellB.neighbor(descriptors::c<DESCRIPTOR>(1));
    auto cellB2 = cellB.neighbor(descriptors::c<DESCRIPTOR>(2));
    auto cellB3 = cellB.neighbor(descriptors::c<DESCRIPTOR>(3));
    auto cellB4 = cellB.neighbor(descriptors::c<DESCRIPTOR>(4));
    auto cellB5 = cellB.neighbor(descriptors::c<DESCRIPTOR>(5));
    auto cellB6 = cellB.neighbor(descriptors::c<DESCRIPTOR>(6));
    auto cellB7 = cellB.neighbor(descriptors::c<DESCRIPTOR>(7));
    auto cellB8 = cellB.neighbor(descriptors::c<DESCRIPTOR>(8));
    auto cellB9 = cellB.neighbor(descriptors::c<DESCRIPTOR>(9));
    auto cellB10 = cellB.neighbor(descriptors::c<DESCRIPTOR>(10));
    auto cellB11 = cellB.neighbor(descriptors::c<DESCRIPTOR>(11));
    auto cellB12 = cellB.neighbor(descriptors::c<DESCRIPTOR>(12));
    auto cellB13 = cellB.neighbor(descriptors::c<DESCRIPTOR>(13));
    auto cellB14 = cellB.neighbor(descriptors::c<DESCRIPTOR>(14));
    auto cellB15 = cellB.neighbor(descriptors::c<DESCRIPTOR>(15));
    auto cellB16 = cellB.neighbor(descriptors::c<DESCRIPTOR>(16));
    auto cellB17 = cellB.neighbor(descriptors::c<DESCRIPTOR>(17));
    auto cellB18 = cellB.neighbor(descriptors::c<DESCRIPTOR>(18));
 
    V gradMuRhoX = 1./12. * (
    -     cellB6.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellB7.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellB4.template getField<descriptors::CHEM_POTENTIAL>()
    - 2 * cellB1.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellB5.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellB14.template getField<descriptors::CHEM_POTENTIAL>()
    + 2 * cellB10.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellB13.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellB16.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellB15.template getField<descriptors::CHEM_POTENTIAL>());
 
    V gradMuRhoY = 1./12. * (
    -     cellB4.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellB14.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellB8.template getField<descriptors::CHEM_POTENTIAL>()
    - 2 * cellB2.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellB9.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellB18.template getField<descriptors::CHEM_POTENTIAL>()
    + 2 * cellB11.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellB17.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellB5.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellB13.template getField<descriptors::CHEM_POTENTIAL>());
 
    V gradMuRhoZ = 1./12. * (
    -     cellB8.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellB18.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellB6.template getField<descriptors::CHEM_POTENTIAL>()
    - 2 * cellB3.template getField<descriptors::CHEM_POTENTIAL>()
    -     cellB16.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellB7.template getField<descriptors::CHEM_POTENTIAL>()
    + 2 * cellB12.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellB15.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellB9.template getField<descriptors::CHEM_POTENTIAL>()
    +     cellB17.template getField<descriptors::CHEM_POTENTIAL>());
 
    V psi = 0.;
    V gradMuPsiX = 0.;
    V gradMuPsiY = 0.;
    V gradMuPsiZ = 0.;
 
    V forceX = - rho*gradMuRhoX - phi*gradMuPhiX - psi*gradMuPsiX;
    V forceY = - rho*gradMuRhoY - phi*gradMuPhiY - psi*gradMuPsiY;
    V forceZ = - rho*gradMuRhoZ - phi*gradMuPhiZ - psi*gradMuPsiZ;
 
    cellB.template setField<descriptors::FORCE>({forceX, forceY, forceZ});
 
    V u[3];
    cellB.computeU(u);
 
    cellA.template setField<descriptors::FORCE>(u);
 
  }

Member Data Documentation

scope

constexpr OperatorScope olb::ForceCoupling3D::scope = OperatorScope::PerCell

staticconstexpr

Definition at line 117 of file freeEnergyCoupling3D.h.

---

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

• src/dynamics/freeEnergyCoupling3D.h