olb::dynamics::CSE< dynamics::Tuple< T, descriptors::D2Q9< FIELDS... >, momenta::Tuple< momenta::BulkDensity, momenta::BulkMomentum, momenta::BulkStress, momenta::DefineToNEq >, equilibria::SecondOrder, collision::ConstRhoBGK, dynamics::DefaultCombination > > Struct Template Reference

#include <624cb65402c15187e66409fd86ad9d5849ee280b2bfb2a16f451c6ce8ec0d9ae.cse.h>

Collaboration diagram for olb::dynamics::CSE< dynamics::Tuple< T, descriptors::D2Q9< FIELDS... >, momenta::Tuple< momenta::BulkDensity, momenta::BulkMomentum, momenta::BulkStress, momenta::DefineToNEq >, equilibria::SecondOrder, collision::ConstRhoBGK, 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::D2Q9< FIELDS... >, momenta::Tuple< momenta::BulkDensity, momenta::BulkMomentum, momenta::BulkStress, momenta::DefineToNEq >, equilibria::SecondOrder, collision::ConstRhoBGK, dynamics::DefaultCombination > >

Definition at line 38 of file 624cb65402c15187e66409fd86ad9d5849ee280b2bfb2a16f451c6ce8ec0d9ae.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::D2Q9< FIELDS... >, momenta::Tuple< momenta::BulkDensity, momenta::BulkMomentum, momenta::BulkStress, momenta::DefineToNEq >, equilibria::SecondOrder, collision::ConstRhoBGK, dynamics::DefaultCombination > >::collide ( CELL & cell, PARAMETERS & parameters )

inline

Definition at line 40 of file 624cb65402c15187e66409fd86ad9d5849ee280b2bfb2a16f451c6ce8ec0d9ae.cse.h.

                                                                          {
auto x9 = parameters.template get<descriptors::OMEGA>();
auto x10 = parameters.template get<statistics::AVERAGE_RHO>();
auto x11 = x9 + V{-1};
auto x12 = cell[0] + cell[1] + cell[2] + cell[3] + cell[4] + cell[5] + cell[6] + cell[7] + cell[8];
auto x13 = x12 + V{1};
auto x14 = V{1} / ((x13)*(x13));
auto x15 = V{1.5}*x14;
auto x16 = cell[1] - cell[5];
auto x17 = -cell[4] + cell[8];
auto x18 = -cell[3] + cell[7] + x16 + x17;
auto x19 = x18*x18;
auto x20 = x15*x19;
auto x21 = cell[2] - cell[6];
auto x22 = cell[3] - cell[7] + x16 + x21;
auto x23 = -x22;
auto x24 = x15*(x23*x23) + V{-1};
auto x25 = x20 + x24;
auto x26 = V{1} / (x13);
auto x27 = -x26*(x10 + V{-1}) + V{1};
auto x28 = x12 + V{1};
auto x29 = x27*x28;
auto x30 = V{0.0277777777777778}*cell[0] + V{0.0277777777777778}*cell[1] + V{0.0277777777777778}*cell[2] + V{0.0277777777777778}*cell[3] + V{0.0277777777777778}*cell[4] + V{0.0277777777777778}*cell[5] + V{0.0277777777777778}*cell[6] + V{0.0277777777777778}*cell[7] + V{0.0277777777777778}*cell[8] + V{0.0277777777777778};
auto x31 = V{3}*cell[3];
auto x32 = V{3}*cell[7];
auto x33 = V{3}*cell[1] - V{3}*cell[5];
auto x34 = V{3}*cell[2] - V{3}*cell[6] + x31 - x32 + x33;
auto x35 = x26*x34;
auto x36 = V{4.5}*x14;
auto x37 = V{2}*cell[1] - V{2}*cell[5] + x17 + x21;
auto x38 = -x37;
auto x39 = V{1} - x20;
auto x40 = x26*(-V{3}*cell[4] + V{3}*cell[8] - x31 + x32 + x33);
auto x41 = x22*x22;
auto x42 = x15*x41;
auto x43 = x40 - x42;
auto x44 = x39 + x43;
auto x45 = x35 + x36*(x38*x38) + x44;
auto x46 = V{0.111111111111111}*cell[0] + V{0.111111111111111}*cell[1] + V{0.111111111111111}*cell[2] + V{0.111111111111111}*cell[3] + V{0.111111111111111}*cell[4] + V{0.111111111111111}*cell[5] + V{0.111111111111111}*cell[6] + V{0.111111111111111}*cell[7] + V{0.111111111111111}*cell[8] + V{0.111111111111111};
auto x47 = V{3}*x14;
auto x48 = x39 + x41*x47;
auto x49 = x35 + x48;
auto x50 = -V{2}*cell[3] - cell[4] + V{2}*cell[7] + cell[8] - x21;
auto x51 = x36*(x50*x50);
auto x52 = x25 - x26*x34 + x40 - x51;
auto x53 = x19*x47;
auto x54 = x24 + x40 - x53;
auto x55 = -(-V{4.5}*x14*x37*x37 + x20 + x35 + x40 + x42 + V{-1});
auto x56 = -x35;
auto x57 = x48 + x56;
auto x58 = x44 + x51 + x56;
auto x59 = x43 + x53 + V{1};
auto x0 = -x11*(cell[0] + x25*(V{0.444444444444444}*cell[0] + V{0.444444444444444}*cell[1] + V{0.444444444444444}*cell[2] + V{0.444444444444444}*cell[3] + V{0.444444444444444}*cell[4] + V{0.444444444444444}*cell[5] + V{0.444444444444444}*cell[6] + V{0.444444444444444}*cell[7] + V{0.444444444444444}*cell[8] + V{0.444444444444444}) + V{0.444444444444444}) - V{0.444444444444444}*x25*x29 + V{-0.444444444444444};
auto x1 = -x11*(cell[1] - x30*x45 + V{0.0277777777777778}) + V{0.0277777777777778}*x27*x28*x45 + V{-0.0277777777777778};
auto x2 = -x11*(cell[2] - x46*x49 + V{0.111111111111111}) + V{0.111111111111111}*x27*x28*x49 + V{-0.111111111111111};
auto x3 = -x11*(cell[3] + x30*x52 + V{0.0277777777777778}) - V{0.0277777777777778}*x29*x52 + V{-0.0277777777777778};
auto x4 = -x11*(cell[4] + x46*x54 + V{0.111111111111111}) - V{0.111111111111111}*x29*x54 + V{-0.111111111111111};
auto x5 = -x11*(cell[5] - x30*x55 + V{0.0277777777777778}) + V{0.0277777777777778}*x27*x28*x55 + V{-0.0277777777777778};
auto x6 = -x11*(cell[6] - x46*x57 + V{0.111111111111111}) + V{0.111111111111111}*x27*x28*x57 + V{-0.111111111111111};
auto x7 = -x11*(cell[7] - x30*x58 + V{0.0277777777777778}) + V{0.0277777777777778}*x27*x28*x58 + V{-0.0277777777777778};
auto x8 = -x11*(cell[8] - x46*x59 + V{0.111111111111111}) + V{0.111111111111111}*x27*x28*x59 + V{-0.111111111111111};
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;
return { -x10 + x12 + V{2}, V{1}*x14*(x19 + x41) };
}

---

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

• src/cse/dynamics/624cb65402c15187e66409fd86ad9d5849ee280b2bfb2a16f451c6ce8ec0d9ae.cse.h