40 {
41auto x12 = parameters.template get<descriptors::OMEGA>();
42auto x9 = cell.template getFieldComponent<momenta::FixedDensity::RHO>(0);
43auto x10 = cell.template getFieldComponent<momenta::FixedPressureMomentum<1, -1>::VELOCITY>(0);
44auto x11 = x12 + V{-1};
45auto x13 = x10*x10;
46auto x14 = V{1.5}*x13;
47auto x15 = V{1} / (x9);
48auto x16 = x15*(cell[0] + cell[2] + V{2}*cell[3] + V{2}*cell[4] + V{2}*cell[5] + cell[6] + V{1});
49auto x17 = V{1} - x16;
50auto x18 = x17*x17;
51auto x19 = V{1.5}*x18;
52auto x20 = x14 + x19;
53auto x21 = V{0.0277777777777778}*x12;
54auto x22 = V{3}*x10;
55auto x23 = x10 + x16 + V{-1};
56auto x24 = -x23;
57auto x25 = x15*(V{3}*cell[0] + V{3}*cell[2] + V{6}*cell[3] + V{6}*cell[4] + V{6}*cell[5] + V{3}*cell[6] + V{3});
58auto x26 = x25 + V{-4};
59auto x27 = x20 + x26;
60auto x28 = V{0.111111111111111}*x12;
61auto x29 = -x22;
62auto x30 = V{3}*x13 - x19 + V{1};
63auto x31 = x10 + x17;
64auto x32 = -x31;
65auto x33 = x20 - x25 + V{2};
66auto x34 = -x17;
67auto x0 = -cell[0]*x11 - V{0.444444444444444}*x12*(x9*(x20 + V{-1}) + V{1});
68auto x1 = -(cell[1]*x11 + x21*(x9*(x22 + x27 - V{4.5}*x24*x24) + V{1}));
69auto x2 = -cell[2]*x11 + x28*(x9*(x29 + x30) + V{-1});
70auto x3 = -(cell[3]*x11 + x21*(x9*(x22 + x33 - V{4.5}*x32*x32) + V{1}));
71auto x4 = -(cell[4]*x11 + x28*(x9*(x33 - V{4.5}*x34*x34) + V{1}));
72auto x5 = -(cell[5]*x11 + x21*(x9*(x29 + x33 - V{4.5}*x23*x23) + V{1}));
73auto x6 = -cell[6]*x11 + x28*(x9*(x22 + x30) + V{-1});
74auto x7 = -(cell[7]*x11 + x21*(x9*(x27 + x29 - V{4.5}*x31*x31) + V{1}));
75auto x8 = -cell[8]*x11 - x28*(x9*(x14 - V{3}*x18 + x26) + V{1});
76cell[0] = x0;
77cell[1] = x1;
78cell[2] = x2;
79cell[3] = x3;
80cell[4] = x4;
81cell[5] = x5;
82cell[6] = x6;
83cell[7] = x7;
84cell[8] = x8;
85return { x9, x13 + V{1}*x18 };
86}