48   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in RiemannOutflow()  local
50     QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, dXdx, normal);  in RiemannOutflow()
52     const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]};  in RiemannOutflow()
66     StatePhysicalGradientFromReference_Boundary(Q, i, gas, s_int, state_var, Grad_q, dXdx, grad_s);  in RiemannOutflow()
79     for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j];  in RiemannOutflow()
83       StoredValuesPack(Q, i, 0, 5, qi, jac_data_sur);  in RiemannOutflow()
84       StoredValuesPack(Q, i, 5, 6, kmstress, jac_data_sur);  in RiemannOutflow()
117   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in RiemannOutflow_Jacobian()  local
119     QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, dXdx, normal);  in RiemannOutflow_Jacobian()
123     StoredValuesUnpack(Q, i, 0, 5, jac_data_sur, qi);  in RiemannOutflow_Jacobian()
124     StoredValuesUnpack(Q, i, 5, 6, jac_data_sur, kmstress);  in RiemannOutflow_Jacobian()
125     for (int j = 0; j < 5; j++) dqi[j] = dq[j][i];  in RiemannOutflow_Jacobian()
147 …  StatePhysicalGradientFromReference_Boundary(Q, i, gas, s_int, state_var, Grad_dq, dXdx, grad_ds);  in RiemannOutflow_Jacobian()
161     for (int j = 0; j < 5; j++) v[j][i] = -wdetJb * dFlux[j];  in RiemannOutflow_Jacobian()
196   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in PressureOutflow()  local
197     const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]};  in PressureOutflow()
202     QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, dXdx, normal);  in PressureOutflow()
206     StatePhysicalGradientFromReference_Boundary(Q, i, gas, s, state_var, Grad_q, dXdx, grad_s);  in PressureOutflow()
220     for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j];  in PressureOutflow()
224       StoredValuesPack(Q, i, 0, 5, qi, jac_data_sur);  in PressureOutflow()
225       StoredValuesPack(Q, i, 5, 6, kmstress, jac_data_sur);  in PressureOutflow()
258   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in PressureOutflow_Jacobian()  local
260     QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, dXdx, normal);  in PressureOutflow_Jacobian()
264     StoredValuesUnpack(Q, i, 0, 5, jac_data_sur, qi);  in PressureOutflow_Jacobian()
265     StoredValuesUnpack(Q, i, 5, 6, jac_data_sur, kmstress);  in PressureOutflow_Jacobian()
266     for (int j = 0; j < 5; j++) dqi[j] = dq[j][i];  in PressureOutflow_Jacobian()
274     StatePhysicalGradientFromReference_Boundary(Q, i, gas, s, state_var, Grad_dq, dXdx, grad_ds);  in PressureOutflow_Jacobian()
289     for (int j = 0; j < 5; j++) v[j][i] = -wdetJb * dFlux[j];  in PressureOutflow_Jacobian()