olb::dynamics::CSE< dynamics::Tuple< T, descriptors::D3Q19< FIELDS... >, momenta::Tuple< momenta::FixedDensity, momenta::FixedPressureMomentum< 2, 1 >, momenta::BulkStress, momenta::DefineSeparately >, equilibria::SecondOrder, collision::BGK, dynamics::DefaultCombination > > Struct Template Reference

#include <3527b2ad793d4247ffa10d96c60c6548b8f6b60401d2a4b74e03ea6d6ec1ade4.cse.h>

Collaboration diagram for olb::dynamics::CSE< dynamics::Tuple< T, descriptors::D3Q19< FIELDS... >, momenta::Tuple< momenta::FixedDensity, momenta::FixedPressureMomentum< 2, 1 >, momenta::BulkStress, momenta::DefineSeparately >, equilibria::SecondOrder, collision::BGK, dynamics::DefaultCombination > >:

Public Member Functions
template<concepts::Cell CELL, concepts::Parameters PARAMETERS, concepts::BaseType V = typename CELL::value_t>
CellStatistic< V >	collide (CELL &cell, PARAMETERS &parameters) any_platform

Detailed Description

template<typename T, typename... FIELDS>
struct olb::dynamics::CSE< dynamics::Tuple< T, descriptors::D3Q19< FIELDS... >, momenta::Tuple< momenta::FixedDensity, momenta::FixedPressureMomentum< 2, 1 >, momenta::BulkStress, momenta::DefineSeparately >, equilibria::SecondOrder, collision::BGK, dynamics::DefaultCombination > >

Definition at line 38 of file 3527b2ad793d4247ffa10d96c60c6548b8f6b60401d2a4b74e03ea6d6ec1ade4.cse.h.

Member Function Documentation

collide()

template<typename T , typename... FIELDS>

template<concepts::Cell CELL, concepts::Parameters PARAMETERS, concepts::BaseType V = typename CELL::value_t>

CellStatistic< V > olb::dynamics::CSE< dynamics::Tuple< T, descriptors::D3Q19< FIELDS... >, momenta::Tuple< momenta::FixedDensity, momenta::FixedPressureMomentum< 2, 1 >, momenta::BulkStress, momenta::DefineSeparately >, equilibria::SecondOrder, collision::BGK, dynamics::DefaultCombination > >::collide ( CELL & cell, PARAMETERS & parameters )

inline

Definition at line 40 of file 3527b2ad793d4247ffa10d96c60c6548b8f6b60401d2a4b74e03ea6d6ec1ade4.cse.h.

                                                                          {
auto x21 = cell.template getFieldComponent<momenta::FixedPressureMomentum<2, 1>::VELOCITY>(1);
auto x20 = cell.template getFieldComponent<momenta::FixedPressureMomentum<2, 1>::VELOCITY>(0);
auto x23 = parameters.template get<descriptors::OMEGA>();
auto x19 = cell.template getFieldComponent<momenta::FixedDensity::RHO>(0);
auto x22 = x23 + V{-1};
auto x24 = x20*x20;
auto x25 = V{1.5}*x24;
auto x26 = x21*x21;
auto x27 = V{1.5}*x26;
auto x28 = x25 + x27;
auto x29 = V{1} / (x19);
auto x30 = x29*(cell[0] + cell[10] + cell[11] + V{2}*cell[12] + cell[13] + cell[14] + V{2}*cell[15] + V{2}*cell[17] + cell[1] + cell[2] + cell[4] + cell[5] + V{2}*cell[7] + V{2}*cell[9] + V{1});
auto x31 = V{1} - x30;
auto x32 = -x31;
auto x33 = x32*x32;
auto x34 = V{1.5}*x33;
auto x35 = x34 + V{-1};
auto x36 = x28 + x35;
auto x37 = V{0.0555555555555556}*x23;
auto x38 = V{3}*x20;
auto x39 = V{3}*x24;
auto x40 = V{3}*x21;
auto x41 = V{3}*x26;
auto x42 = x31*x31;
auto x43 = x29*(V{3}*cell[0] + V{3}*cell[10] + V{3}*cell[11] + V{6}*cell[12] + V{3}*cell[13] + V{3}*cell[14] + V{6}*cell[15] + V{6}*cell[17] + V{3}*cell[1] + V{3}*cell[2] + V{3}*cell[4] + V{3}*cell[5] + V{6}*cell[7] + V{6}*cell[9] + V{3});
auto x44 = x28 + x34 + x43 + V{-4};
auto x45 = V{0.0277777777777778}*x23;
auto x46 = x20 + x21;
auto x47 = V{4.5}*(x46*x46);
auto x48 = x36 + x38;
auto x49 = -x40;
auto x50 = x20 - x21;
auto x51 = -x50;
auto x52 = x30 + V{-1};
auto x53 = x20 + x52;
auto x54 = -x53;
auto x55 = x20 + x31;
auto x56 = -x55;
auto x57 = x28 - x43 + V{2};
auto x58 = x34 + x57;
auto x59 = x21 + x52;
auto x60 = -x59;
auto x61 = x21 + x31;
auto x62 = -x61;
auto x63 = V{1} - V{1.5}*x42;
auto x64 = -x27 + x38 + x63;
auto x65 = -x25 + x40;
auto x66 = -x38;
auto x0 = -cell[0]*x22 - V{0.333333333333333}*x23*(x19*x36 + V{1});
auto x1 = -cell[1]*x22 - x37*(x19*(x27 + x35 + x38 - x39) + V{1});
auto x2 = -cell[2]*x22 - x37*(x19*(x25 + x35 + x40 - x41) + V{1});
auto x3 = -cell[3]*x22 - x37*(x19*(-V{4.5}*x42 + x44) + V{1});
auto x4 = -cell[4]*x22 - x45*(x19*(x40 - x47 + x48) + V{1});
auto x5 = -(cell[5]*x22 + x45*(x19*(x48 + x49 - V{4.5}*x51*x51) + V{1}));
auto x6 = -(cell[6]*x22 + x45*(x19*(x38 + x44 - V{4.5}*x54*x54) + V{1}));
auto x7 = -(cell[7]*x22 + x45*(x19*(x38 + x58 - V{4.5}*x56*x56) + V{1}));
auto x8 = -(cell[8]*x22 + x45*(x19*(x40 + x44 - V{4.5}*x60*x60) + V{1}));
auto x9 = -(cell[9]*x22 + x45*(x19*(x40 + x58 - V{4.5}*x62*x62) + V{1}));
auto x10 = -cell[10]*x22 + x37*(x19*(x39 + x64) + V{-1});
auto x11 = -cell[11]*x22 + x37*(x19*(x41 + x63 + x65) + V{-1});
auto x12 = -cell[12]*x22 - x37*(x19*(-V{3}*x33 + x57) + V{1});
auto x13 = -cell[13]*x22 + x45*(x19*(x47 + x64 + x65) + V{-1});
auto x14 = -(cell[14]*x22 + x45*(x19*(x36 + x40 + x66 - V{4.5}*x50*x50) + V{1}));
auto x15 = -(cell[15]*x22 + x45*(x19*(x58 + x66 - V{4.5}*x53*x53) + V{1}));
auto x16 = -(cell[16]*x22 + x45*(x19*(x44 + x66 - V{4.5}*x55*x55) + V{1}));
auto x17 = -(cell[17]*x22 + x45*(x19*(x49 + x58 - V{4.5}*x59*x59) + V{1}));
auto x18 = -(cell[18]*x22 + x45*(x19*(x44 + x49 - V{4.5}*x61*x61) + V{1}));
cell[0] = x0;
cell[1] = x1;
cell[2] = x2;
cell[3] = x3;
cell[4] = x4;
cell[5] = x5;
cell[6] = x6;
cell[7] = x7;
cell[8] = x8;
cell[9] = x9;
cell[10] = x10;
cell[11] = x11;
cell[12] = x12;
cell[13] = x13;
cell[14] = x14;
cell[15] = x15;
cell[16] = x16;
cell[17] = x17;
cell[18] = x18;
return { x19, x24 + x26 + V{1}*x42 };
}

---

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

• src/cse/dynamics/3527b2ad793d4247ffa10d96c60c6548b8f6b60401d2a4b74e03ea6d6ec1ade4.cse.h