65   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in ICsChannel()  local
66     const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]};  in ICsChannel()
70     for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j];  in ICsChannel()
89   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in Channel_Inflow()  local
91     QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, NULL, norm);  in Channel_Inflow()
99     const CeedScalar x[3] = {0, X[1][i], X[2][i]};  in Channel_Inflow()
108     for (CeedInt j = 0; j < 5; j++) q_inside[j] = q[j][i];  in Channel_Inflow()
120     for (CeedInt j = 0; j < 5; j++) v[j][i] = 0.;  in Channel_Inflow()
126     v[0][i] -= wdetJb * rho_in * u_normal;  in Channel_Inflow()
129 …for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho_in * u_normal * s_exact.Y.velocity[j]…  in Channel_Inflow()
132     v[4][i] -= wdetJb * u_normal * (E + P);  in Channel_Inflow()
148   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in Channel_Outflow()  local
150     QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, NULL, norm);  in Channel_Outflow()
153     const CeedScalar rho  = q[0][i];  in Channel_Outflow()
154     const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho};  in Channel_Outflow()
155     const CeedScalar E    = q[4][i];  in Channel_Outflow()
159     for (CeedInt j = 0; j < 5; j++) v[j][i] = 0.;  in Channel_Outflow()
166     v[0][i] -= wdetJb * rho * u_normal;  in Channel_Outflow()
169     for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho * u_normal * u[j] + norm[j] * P);  in Channel_Outflow()
172     v[4][i] -= wdetJb * u_normal * (E + P);  in Channel_Outflow()