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 = V{0.125}*x10;
49auto x15 = x13*(x12 + V{-1}) + V{0.125};
50auto x16 = x11*(V{1}*cell[1] + x15);
51auto x17 = V{4}*x8;
52auto x18 = x13*(x17 + V{-1}) + V{0.125};
53auto x19 = x11*(V{1}*cell[2] + x18);
54auto x20 = V{4}*x9;
55auto x21 = x20 + V{1};
56auto x22 = V{0.0625}*x10;
57auto x23 = x20 + V{-1};
58auto x24 = x11*(V{0.5}*cell[3] - V{0.5}*cell[6] + x21*x22 + x22*x23);
59cell[0] = V{0.25}*x10 + V{-0.25};
60cell[1] = -x15 - x16;
61cell[2] = -x18 - x19;
62cell[3] = -x13*x23 - x24 + V{-0.125};
63cell[4] = x13*(x12 + V{1}) + x16 + V{-0.125};
64cell[5] = x13*(x17 + V{1}) + x19 + V{-0.125};
65cell[6] = x13*x21 + x24 + V{-0.125};
66return { x10, x7*x7 + x8*x8 + x9*x9 };
67}