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