430 …NCTION_HELPER void FluxInviscid(NewtonianIdealGasContext gas, State s, StateConservative Flux[3]) {  in FluxInviscid()
432     Flux[i].density = s.U.momentum[i];  in FluxInviscid()
433 …for (CeedInt j = 0; j < 3; j++) Flux[i].momentum[j] = s.U.momentum[i] * s.Y.velocity[j] + s.Y.pres…  in FluxInviscid()
434     Flux[i].E_total = (s.U.E_total + s.Y.pressure) * s.Y.velocity[i];  in FluxInviscid()
449   StateConservative Flux[3], Flux_dot_n = {0};  in FluxInviscidDotNormal()  local
450   FluxInviscid(gas, s, Flux);  in FluxInviscidDotNormal()
452     Flux_dot_n.density += Flux[i].density * normal[i];  in FluxInviscidDotNormal()
453     for (CeedInt j = 0; j < 3; j++) Flux_dot_n.momentum[j] += Flux[i].momentum[j] * normal[i];  in FluxInviscidDotNormal()
454     Flux_dot_n.E_total += Flux[i].E_total * normal[i];  in FluxInviscidDotNormal()
481 …tateConservative F_inviscid[3], CeedScalar stress[3][3], CeedScalar Fe[3], CeedScalar Flux[5][3]) {  in FluxTotal()
483     Flux[0][j] = F_inviscid[j].density;  in FluxTotal()
484     for (CeedInt k = 0; k < 3; k++) Flux[k + 1][j] = F_inviscid[j].momentum[k] - stress[k][j];  in FluxTotal()
485     Flux[4][j] = F_inviscid[j].E_total + Fe[j];  in FluxTotal()
490                                               const CeedScalar normal[3], CeedScalar Flux[5]) {  in FluxTotal_Boundary()
491   for (CeedInt j = 0; j < 5; j++) Flux[j] = 0.;  in FluxTotal_Boundary()
493     Flux[0] += F_inviscid[j].density * normal[j];  in FluxTotal_Boundary()
495       Flux[k + 1] += (F_inviscid[j].momentum[k] - stress[k][j]) * normal[j];  in FluxTotal_Boundary()
497     Flux[4] += (F_inviscid[j].E_total + Fe[j]) * normal[j];  in FluxTotal_Boundary()
502 …                                                  const CeedScalar normal[3], CeedScalar Flux[5]) {  in FluxTotal_RiemannBoundary()
503   Flux[0] = F_inviscid_normal.density;  in FluxTotal_RiemannBoundary()
504   for (CeedInt k = 0; k < 3; k++) Flux[k + 1] = F_inviscid_normal.momentum[k];  in FluxTotal_RiemannBoundary()
505   Flux[4] = F_inviscid_normal.E_total;  in FluxTotal_RiemannBoundary()
508       Flux[k + 1] -= stress[k][j] * normal[j];  in FluxTotal_RiemannBoundary()
510     Flux[4] += Fe[j] * normal[j];  in FluxTotal_RiemannBoundary()