40 {
41auto x7 = cell.template getFieldComponent<momenta::FixedDensity::RHO>(0);
42auto x8 = parameters.template get<descriptors::OMEGA>();
43auto x5 = cell.template getFieldComponent<descriptors::VELOCITY>(0);
44auto x6 = cell.template getFieldComponent<descriptors::VELOCITY>(1);
45auto x9 = V{3}*x5;
46auto x10 = x9 + V{-1};
47auto x11 = V{0.0555555555555556}*x7;
48auto x12 = V{3}*x6;
49auto x13 = x12 + V{-1};
50auto x14 = x9 + V{1};
51auto x15 = x12 + V{1};
52auto x16 = V{0.0277777777777778}*x7;
53auto x17 = x11 + x14*x16 + x15*x16 + V{1.85037170770859e-17};
54auto x18 = -x10*x16 - x13*x16 + x17;
55auto x19 = x10*x18;
56auto x20 = x8 + V{-1};
57auto x21 = -x10;
58auto x22 = -V{0.5}*x13*x16 + V{0.5}*x16*x21 + V{0.5}*x17;
59auto x23 = V{0.0833333333333333}*x7;
60auto x24 = V{1}*cell[3] - x14*x23 + V{0.166666666666667};
61auto x25 = x13*x18;
62auto x26 = x15*x18;
63auto x27 = V{0.5}*x20*(cell[2] - cell[4] + x25 + x26);
64auto x28 = V{0.166666666666667}*x7;
65cell[0] = -x10*x11 - x11*x13 + V{0.0555555555555556}*x14*x7 + V{0.0555555555555556}*x15*x7 + V{0.111111111111111}*x7 + V{-0.333333333333333};
66cell[1] = -x19 + x20*(-x14*x22 + x21*x22 - x21*x23 + x24) + V{-0.166666666666667};
67cell[2] = -x25 - x27 + V{-0.166666666666667};
68cell[3] = x14*x18 - x20*(x10*x23 - V{0.5}*x14*x18 - V{0.5}*x19 + x24) + V{-0.166666666666667};
69cell[4] = x26 + x27 + V{-0.166666666666667};
70return { -x10*x28 - x13*x28 + x14*x28 + x15*x28 + V{0.333333333333333}*x7 + V{1.11022302462516e-16}, x5*x5 + x6*x6 };
71}