50   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in ChildStatsCollection()  local
51     const CeedScalar wdetJ = q_data[0][i] * delta_t;  in ChildStatsCollection()
53     const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]};  in ChildStatsCollection()
56     v[TURB_MEAN_DENSITY][i]                    = wdetJ * s.U.density;  in ChildStatsCollection()
57     v[TURB_MEAN_PRESSURE][i]                   = wdetJ * s.Y.pressure;  in ChildStatsCollection()
58     v[TURB_MEAN_PRESSURE_SQUARED][i]           = wdetJ * Square(s.Y.pressure);  in ChildStatsCollection()
59     v[TURB_MEAN_PRESSURE_VELOCITY_X][i]        = wdetJ * s.Y.pressure * s.Y.velocity[0];  in ChildStatsCollection()
60     v[TURB_MEAN_PRESSURE_VELOCITY_Y][i]        = wdetJ * s.Y.pressure * s.Y.velocity[1];  in ChildStatsCollection()
61     v[TURB_MEAN_PRESSURE_VELOCITY_Z][i]        = wdetJ * s.Y.pressure * s.Y.velocity[2];  in ChildStatsCollection()
62     v[TURB_MEAN_DENSITY_TEMPERATURE][i]        = wdetJ * s.U.density * s.Y.temperature;  in ChildStatsCollection()
63 …v[TURB_MEAN_DENSITY_TEMPERATURE_FLUX_X][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[…  in ChildStatsCollection()
64 …v[TURB_MEAN_DENSITY_TEMPERATURE_FLUX_Y][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[…  in ChildStatsCollection()
65 …v[TURB_MEAN_DENSITY_TEMPERATURE_FLUX_Z][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[…  in ChildStatsCollection()
66     v[TURB_MEAN_MOMENTUM_X][i]                 = wdetJ * s.U.momentum[0];  in ChildStatsCollection()
67     v[TURB_MEAN_MOMENTUM_Y][i]                 = wdetJ * s.U.momentum[1];  in ChildStatsCollection()
68     v[TURB_MEAN_MOMENTUM_Z][i]                 = wdetJ * s.U.momentum[2];  in ChildStatsCollection()
69     v[TURB_MEAN_MOMENTUMFLUX_XX][i]            = wdetJ * s.U.momentum[0] * s.Y.velocity[0];  in ChildStatsCollection()
70     v[TURB_MEAN_MOMENTUMFLUX_YY][i]            = wdetJ * s.U.momentum[1] * s.Y.velocity[1];  in ChildStatsCollection()
71     v[TURB_MEAN_MOMENTUMFLUX_ZZ][i]            = wdetJ * s.U.momentum[2] * s.Y.velocity[2];  in ChildStatsCollection()
72     v[TURB_MEAN_MOMENTUMFLUX_YZ][i]            = wdetJ * s.U.momentum[1] * s.Y.velocity[2];  in ChildStatsCollection()
73     v[TURB_MEAN_MOMENTUMFLUX_XZ][i]            = wdetJ * s.U.momentum[0] * s.Y.velocity[2];  in ChildStatsCollection()
74     v[TURB_MEAN_MOMENTUMFLUX_XY][i]            = wdetJ * s.U.momentum[0] * s.Y.velocity[1];  in ChildStatsCollection()
75     v[TURB_MEAN_VELOCITY_X][i]                 = wdetJ * s.Y.velocity[0];  in ChildStatsCollection()
76     v[TURB_MEAN_VELOCITY_Y][i]                 = wdetJ * s.Y.velocity[1];  in ChildStatsCollection()
77     v[TURB_MEAN_VELOCITY_Z][i]                 = wdetJ * s.Y.velocity[2];  in ChildStatsCollection()
105   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in ChildStatsCollectionMMSTest()  local
106     const CeedScalar wdetJ  = q_data[0][i];  in ChildStatsCollectionMMSTest()
107     const CeedScalar x_i[3] = {x[0][i], x[1][i], x[2][i]};  in ChildStatsCollectionMMSTest()
110     v[0][i] = wdetJ * (ChildStatsCollectionTest_Exact(x_i) + t - 0.5) * 4 * Cube(x_i[2]);  in ChildStatsCollectionMMSTest()
111     for (int j = 1; j < 22; j++) v[j][i] = 0;  in ChildStatsCollectionMMSTest()
122   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in ChildStatsCollectionMMSTest_Error()  local
123     const CeedScalar wdetJ  = q_data[0][i];  in ChildStatsCollectionMMSTest_Error()
124     const CeedScalar x_i[3] = {x[0][i], x[1][i], x[2][i]};  in ChildStatsCollectionMMSTest_Error()
126     v[0][i] = wdetJ * Square(ChildStatsCollectionTest_Exact(x_i) - q[0][i]);  in ChildStatsCollectionMMSTest_Error()