425 …LPER void FluxInviscid(const NewtonianIGProperties gas, const State s, StateConservative Flux[3]) {  in FluxInviscid()
427     Flux[i].density = s.U.momentum[i];  in FluxInviscid()
428 …for (CeedInt j = 0; j < 3; j++) Flux[i].momentum[j] = s.U.momentum[i] * s.Y.velocity[j] + s.Y.pres…  in FluxInviscid()
429     Flux[i].E_total = (s.U.E_total + s.Y.pressure) * s.Y.velocity[i];  in FluxInviscid()
444   StateConservative Flux[3], Flux_dot_n = {0};  in FluxInviscidDotNormal()  local
445   FluxInviscid(gas, s, Flux);  in FluxInviscidDotNormal()
447     Flux_dot_n.density += Flux[i].density * normal[i];  in FluxInviscidDotNormal()
448     for (CeedInt j = 0; j < 3; j++) Flux_dot_n.momentum[j] += Flux[i].momentum[j] * normal[i];  in FluxInviscidDotNormal()
449     Flux_dot_n.E_total += Flux[i].E_total * normal[i];  in FluxInviscidDotNormal()
478                                      CeedScalar Flux[5][3]) {  in FluxTotal()
480     Flux[0][j] = F_inviscid[j].density;  in FluxTotal()
481     for (CeedInt k = 0; k < 3; k++) Flux[k + 1][j] = F_inviscid[j].momentum[k] - stress[k][j];  in FluxTotal()
482     Flux[4][j] = F_inviscid[j].E_total + Fe[j];  in FluxTotal()
487                                               const CeedScalar normal[3], CeedScalar Flux[5]) {  in FluxTotal_Boundary()
488   for (CeedInt j = 0; j < 5; j++) Flux[j] = 0.;  in FluxTotal_Boundary()
490     Flux[0] += F_inviscid[j].density * normal[j];  in FluxTotal_Boundary()
492       Flux[k + 1] += (F_inviscid[j].momentum[k] - stress[k][j]) * normal[j];  in FluxTotal_Boundary()
494     Flux[4] += (F_inviscid[j].E_total + Fe[j]) * normal[j];  in FluxTotal_Boundary()
499 …                                                  const CeedScalar normal[3], CeedScalar Flux[5]) {  in FluxTotal_RiemannBoundary()
500   Flux[0] = F_inviscid_normal.density;  in FluxTotal_RiemannBoundary()
501   for (CeedInt k = 0; k < 3; k++) Flux[k + 1] = F_inviscid_normal.momentum[k];  in FluxTotal_RiemannBoundary()
502   Flux[4] = F_inviscid_normal.E_total;  in FluxTotal_RiemannBoundary()
505       Flux[k + 1] -= stress[k][j] * normal[j];  in FluxTotal_RiemannBoundary()
507     Flux[4] += Fe[j] * normal[j];  in FluxTotal_RiemannBoundary()