38   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in DifferentialFilter_RHS()  local
39     const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]};  in DifferentialFilter_RHS()
40     const CeedScalar wdetJ = q_data[0][i];  in DifferentialFilter_RHS()
43     v0[DIFF_FILTER_PRESSURE][i]            = wdetJ * s.Y.pressure;  in DifferentialFilter_RHS()
44     v0[DIFF_FILTER_VELOCITY_X][i]          = wdetJ * s.Y.velocity[0];  in DifferentialFilter_RHS()
45     v0[DIFF_FILTER_VELOCITY_Y][i]          = wdetJ * s.Y.velocity[1];  in DifferentialFilter_RHS()
46     v0[DIFF_FILTER_VELOCITY_Z][i]          = wdetJ * s.Y.velocity[2];  in DifferentialFilter_RHS()
47     v0[DIFF_FILTER_TEMPERATURE][i]         = wdetJ * s.Y.temperature;  in DifferentialFilter_RHS()
48     v1[DIFF_FILTER_VELOCITY_SQUARED_XX][i] = wdetJ * s.Y.velocity[0] * s.Y.velocity[0];  in DifferentialFilter_RHS()
49     v1[DIFF_FILTER_VELOCITY_SQUARED_YY][i] = wdetJ * s.Y.velocity[1] * s.Y.velocity[1];  in DifferentialFilter_RHS()
50     v1[DIFF_FILTER_VELOCITY_SQUARED_ZZ][i] = wdetJ * s.Y.velocity[2] * s.Y.velocity[2];  in DifferentialFilter_RHS()
51     v1[DIFF_FILTER_VELOCITY_SQUARED_YZ][i] = wdetJ * s.Y.velocity[1] * s.Y.velocity[2];  in DifferentialFilter_RHS()
52     v1[DIFF_FILTER_VELOCITY_SQUARED_XZ][i] = wdetJ * s.Y.velocity[0] * s.Y.velocity[2];  in DifferentialFilter_RHS()
53     v1[DIFF_FILTER_VELOCITY_SQUARED_XY][i] = wdetJ * s.Y.velocity[0] * s.Y.velocity[1];  in DifferentialFilter_RHS()
85   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in DifferentialFilter_LHS_N()  local
87       const CeedScalar x_i[3] = {x[0][i], x[1][i], x[2][i]};  in DifferentialFilter_LHS_N()
89       QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx);  in DifferentialFilter_LHS_N()
93 …   km_A_ij[6] = {A_ij_delta[0][i], A_ij_delta[1][i], A_ij_delta[2][i], A_ij_delta[3][i], A_ij_delt…  in DifferentialFilter_LHS_N()
94         const CeedScalar delta      = A_ij_delta[6][i];  in DifferentialFilter_LHS_N()
126       v[j][i] = wdetJ * q[j][i];  in DifferentialFilter_LHS_N()
129 …dq[k] = Grad_q[0 * N + j][i] * dXdx[0][k] + Grad_q[1 * N + j][i] * dXdx[1][k] + Grad_q[2 * N + j][…  in DifferentialFilter_LHS_N()
134         Grad_v[k * N + j][i] = wdetJ * dq_dXdx_a[k];  in DifferentialFilter_LHS_N()
166   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in DifferentialFilter_MMS_RHS()  local
167     const CeedScalar wdetJ = q_data[0][i];  in DifferentialFilter_MMS_RHS()
168     v[0][i]                = wdetJ * q[0][i];  in DifferentialFilter_MMS_RHS()
188   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in DifferentialFilter_MMS_IC()  local
189     const CeedScalar x_i[3] = {x[0][i], x[1][i], x[2][i]};  in DifferentialFilter_MMS_IC()
214     for (CeedInt j = 0; j < 5; j++) q0[j][i] = phi;  in DifferentialFilter_MMS_IC()