43 CEED_QFUNCTION_HELPER void UnpackState_V(StateEntropy s, CeedScalar V[5]) {  in UnpackState_V()
44   V[0] = s.S_density;  in UnpackState_V()
45   for (int i = 0; i < 3; i++) V[i + 1] = s.S_momentum[i];  in UnpackState_V()
46   V[4] = s.S_energy;  in UnpackState_V()
99   StateEntropy     V;  in StateEntropyFromPrimitive()  local
106   V.S_density = (gamma - entropy) / (gamma - 1) - rho_div_p * e_kinetic;  in StateEntropyFromPrimitive()
107   for (int i = 0; i < 3; i++) V.S_momentum[i] = rho_div_p * Y.velocity[i];  in StateEntropyFromPrimitive()
108   V.S_energy = -rho_div_p;  in StateEntropyFromPrimitive()
109   return V;  in StateEntropyFromPrimitive()
130 …ION_HELPER StatePrimitive StatePrimitiveFromEntropy(NewtonianIdealGasContext gas, StateEntropy V) {  in StatePrimitiveFromEntropy()  argument
132   for (int i = 0; i < 3; i++) Y.velocity[i] = -V.S_momentum[i] / V.S_energy;  in StatePrimitiveFromEntropy()
133   Y.temperature              = -1 / (GasConstant(gas) * V.S_energy);  in StatePrimitiveFromEntropy()
136   const CeedScalar entropy   = gamma - (gamma - 1) * (V.S_density - e_kinetic * V.S_energy);  in StatePrimitiveFromEntropy()
137   const CeedScalar log_P     = -(entropy + gamma * log(-V.S_energy)) / (gamma - 1);  in StatePrimitiveFromEntropy()
144   StateEntropy   V = StateEntropyFromPrimitive(gas, s.Y);  in StatePrimitiveFromEntropy_fwd()  local
145 …0; i < 3; i++) dY.velocity[i] = -(dV.S_momentum[i] - V.S_momentum[i] * dV.S_energy / V.S_energy) /…  in StatePrimitiveFromEntropy_fwd()
146   dY.temperature              = dV.S_energy / (GasConstant(gas) * V.S_energy * V.S_energy);  in StatePrimitiveFromEntropy_fwd()
150 …alar dentropy   = (1 - gamma) * (dV.S_density - e_kinetic * dV.S_energy - de_kinetic * V.S_energy);  in StatePrimitiveFromEntropy_fwd()
151 …dY.pressure                 = s.Y.pressure * (-dentropy - gamma * dV.S_energy / V.S_energy) / (gam…  in StatePrimitiveFromEntropy_fwd()
183   StateEntropy     V;  in StateEntropyFromConservative()  local
190   V.S_density = (gamma - entropy) / (gamma - 1) - e_kinetic / p;  in StateEntropyFromConservative()
191   for (int i = 0; i < 3; i++) V.S_momentum[i] = U.momentum[i] / p;  in StateEntropyFromConservative()
192   V.S_energy = -U.density / p;  in StateEntropyFromConservative()
193   return V;  in StateEntropyFromConservative()
215 …LPER StateConservative StateConservativeFromEntropy(NewtonianIdealGasContext gas, StateEntropy V) {  in StateConservativeFromEntropy()  argument
218   for (int i = 0; i < 3; i++) velocity[i] = -V.S_momentum[i] / V.S_energy;  in StateConservativeFromEntropy()
221   const CeedScalar entropy   = gamma - (gamma - 1) * (V.S_density - e_kinetic * V.S_energy);  in StateConservativeFromEntropy()
222   const CeedScalar log_rho   = -(entropy + log(-V.S_energy)) / (gamma - 1);  in StateConservativeFromEntropy()
226   const CeedScalar e_internal = -gas->cv / (GasConstant(gas) * V.S_energy);  in StateConservativeFromEntropy()
234   StateEntropy      V = StateEntropyFromPrimitive(gas, s.Y);  in StateConservativeFromEntropy_fwd()  local
235 …i = 0; i < 3; i++) dvelocity[i] = (-dV.S_momentum[i] - s.Y.velocity[i] * dV.S_energy) / V.S_energy;  in StateConservativeFromEntropy_fwd()
239   const CeedScalar entropy    = gamma - (gamma - 1) * (V.S_density - e_kinetic * V.S_energy);  in StateConservativeFromEntropy_fwd()
240 …const CeedScalar dentropy   = -(gamma - 1) * (dV.S_density - (de_kinetic * V.S_energy + e_kinetic …  in StateConservativeFromEntropy_fwd()
241   const CeedScalar log_rho    = -(entropy + log(-V.S_energy)) / (gamma - 1);  in StateConservativeFromEntropy_fwd()
243   dU.density                  = -rho / (gamma - 1) * (dentropy + dV.S_energy / V.S_energy);  in StateConservativeFromEntropy_fwd()
246   const CeedScalar e_internal  = -gas->cv / (GasConstant(gas) * V.S_energy);  in StateConservativeFromEntropy_fwd()
247 …const CeedScalar de_internal = gas->cv * dV.S_energy / (GasConstant(gas) * V.S_energy * V.S_energy…  in StateConservativeFromEntropy_fwd()
293 CEED_QFUNCTION_HELPER void StateToV(NewtonianIdealGasContext gas, const State input, CeedScalar V[5…  in StateToV()
295   UnpackState_V(state_V, V);  in StateToV()
372 CEED_QFUNCTION_HELPER State StateFromV(NewtonianIdealGasContext gas, const CeedScalar V[5]) {  in StateFromV()
375   state_V.S_density     = V[0];  in StateFromV()
376   state_V.S_momentum[0] = V[1];  in StateFromV()
377   state_V.S_momentum[1] = V[2];  in StateFromV()
378   state_V.S_momentum[2] = V[3];  in StateFromV()
379   state_V.S_energy      = V[4];  in StateFromV()