40 {
41auto x9 = cell.template getFieldComponent<descriptors::VELOCITY>(2);
42auto x10 = cell.template getFieldComponent<momenta::FixedDensity::RHO>(0);
43auto x7 = cell.template getFieldComponent<descriptors::VELOCITY>(0);
44auto x14 = parameters.template get<descriptors::OMEGA>();
45auto x8 = cell.template getFieldComponent<descriptors::VELOCITY>(1);
46auto x11 = x14 + V{-1};
47auto x12 = V{4}*x7;
48auto x13 = x12 + V{1};
49auto x15 = V{0.0625}*x10;
50auto x16 = x12 + V{-1};
51auto x17 = x11*(V{0.5}*cell[1] - V{0.5}*cell[4] + x13*x15 + x15*x16);
52auto x18 = V{0.125}*x10;
53auto x19 = V{4}*x8;
54auto x20 = x18*(x19 + V{-1}) + V{0.125};
55auto x21 = x11*(V{1}*cell[2] + x20);
56auto x22 = V{4}*x9;
57auto x23 = x18*(x22 + V{-1}) + V{0.125};
58auto x24 = x11*(V{1}*cell[3] + x23);
59cell[0] = V{0.25}*x10 + V{-0.25};
60cell[1] = -x16*x18 - x17 + V{-0.125};
61cell[2] = -x20 - x21;
62cell[3] = -x23 - x24;
63cell[4] = x13*x18 + x17 + V{-0.125};
64cell[5] = x18*(x19 + V{1}) + x21 + V{-0.125};
65cell[6] = x18*(x22 + V{1}) + x24 + V{-0.125};
66return { x10, x7*x7 + x8*x8 + x9*x9 };
67}