14   CeedScalar density;  member
25   U[0] = s.density;  in UnpackState_U()
55   return e_internal + e_kinetic + s.Y.pressure / s.U.density;  in TotalSpecificEnthalpy()
61 …de_internal + de_kinetic + ds.Y.pressure / s.U.density - s.Y.pressure / Square(s.U.density) * ds.U…  in TotalSpecificEnthalpy_fwd()
66   for (CeedInt i = 0; i < 3; i++) Y.velocity[i] = U.momentum[i] / U.density;  in StatePrimitiveFromConservative()
68   CeedScalar e_total    = U.E_total / U.density;  in StatePrimitiveFromConservative()
71   Y.pressure            = (HeatCapacityRatio(gas) - 1) * U.density * e_internal;  in StatePrimitiveFromConservative()
78     dY.velocity[i] = (dU.momentum[i] - s.Y.velocity[i] * dU.density) / s.U.density;  in StatePrimitiveFromConservative_fwd()
82   CeedScalar e_total     = s.U.E_total / s.U.density;  in StatePrimitiveFromConservative_fwd()
83   CeedScalar de_total    = (dU.E_total - e_total * dU.density) / s.U.density;  in StatePrimitiveFromConservative_fwd()
87 …dY.pressure            = (HeatCapacityRatio(gas) - 1) * (dU.density * e_internal + s.U.density * d…  in StatePrimitiveFromConservative_fwd()
112   const CeedScalar rho_div_p  = s.U.density / s.Y.pressure;  in StateEntropyFromPrimitive_fwd()
113 …const CeedScalar drho_div_p = (drho * s.Y.pressure - s.U.density * dY.pressure) / Square(s.Y.press…  in StateEntropyFromPrimitive_fwd()
115   CeedScalar dentropy = dY.pressure / s.Y.pressure - gamma * drho / s.U.density;  in StateEntropyFromPrimitive_fwd()
150   U.density = Y.pressure / (GasConstant(gas) * Y.temperature);  in StateConservativeFromPrimitive()
151   for (int i = 0; i < 3; i++) U.momentum[i] = U.density * Y.velocity[i];  in StateConservativeFromPrimitive()
155   U.E_total             = U.density * e_total;  in StateConservativeFromPrimitive()
161 …dU.density = (dY.pressure * s.Y.temperature - s.Y.pressure * dY.temperature) / (GasConstant(gas) *…  in StateConservativeFromPrimitive_fwd()
163     dU.momentum[i] = dU.density * s.Y.velocity[i] + s.U.density * dY.velocity[i];  in StateConservativeFromPrimitive_fwd()
171   dU.E_total             = dU.density * e_total + s.U.density * de_total;  in StateConservativeFromPrimitive_fwd()
178   const CeedScalar e_kinetic  = .5 * Dot3(U.momentum, U.momentum) / U.density;  in StateEntropyFromConservative()
181   const CeedScalar entropy    = log(p) - gamma * log(U.density);  in StateEntropyFromConservative()
185   V.S_energy = -U.density / p;  in StateEntropyFromConservative()
192   const CeedScalar e_kinetic   = .5 * Dot3(s.U.momentum, s.U.momentum) / s.U.density;  in StateEntropyFromConservative_fwd()
193 … CeedScalar de_kinetic  = (Dot3(s.U.momentum, dU.momentum) - e_kinetic * dU.density) / s.U.density;  in StateEntropyFromConservative_fwd()
198   CeedScalar dentropy = dp / p - gamma * dU.density / s.U.density;  in StateEntropyFromConservative_fwd()
204   dV.S_energy = -(dU.density - s.U.density * dp / p) / p;  in StateEntropyFromConservative_fwd()
216   U.density                  = exp(log_rho);  in StateConservativeFromEntropy()
217   for (int i = 0; i < 3; i++) U.momentum[i] = U.density * velocity[i];  in StateConservativeFromEntropy()
220   U.E_total                   = U.density * (e_internal + e_kinetic);  in StateConservativeFromEntropy()
236   dU.density                  = -rho / (gamma - 1) * (dentropy + dV.S_energy / V.S_energy);  in StateConservativeFromEntropy_fwd()
237 …for (int i = 0; i < 3; i++) dU.momentum[i] = dU.density * s.Y.velocity[i] + s.U.density * dvelocit…  in StateConservativeFromEntropy_fwd()
242   dU.E_total                   = dU.density * e_total + s.U.density * (de_internal + de_kinetic);  in StateConservativeFromEntropy_fwd()
266     R.density += a[i] * X[i].density;  in StateConservativeMult()
276   R.density = a * X.density + b * Y.density + c * Z.density;  in StateConservativeAXPBYPCZ()
325   s.U.density     = U[0];  in StateFromU()
336   ds.U.density     = dU[0];  in StateFromU_fwd()
427     Flux[i].density = s.U.momentum[i];  in FluxInviscid()
435     dFlux[i].density = ds.U.momentum[i];  in FluxInviscid_fwd()
447     Flux_dot_n.density += Flux[i].density * normal[i];  in FluxInviscidDotNormal()
459     Flux_dot_n.density += dFlux[i].density * normal[i];  in FluxInviscidDotNormal_fwd()
480     Flux[0][j] = F_inviscid[j].density;  in FluxTotal()
490     Flux[0] += F_inviscid[j].density * normal[j];  in FluxTotal_Boundary()
500   Flux[0] = F_inviscid_normal.density;  in FluxTotal_RiemannBoundary()