Lines Matching refs:U
20 StateConservative U; member
24 CEED_QFUNCTION_HELPER void UnpackState_U(StateConservative s, CeedScalar U[5]) { in UnpackState_U()
25 U[0] = s.density; in UnpackState_U()
26 for (int i = 0; i < 3; i++) U[i + 1] = s.momentum[i]; in UnpackState_U()
27 U[4] = s.E_total; in UnpackState_U()
55 return e_internal + e_kinetic + s.Y.pressure / s.U.density; in TotalSpecificEnthalpy()
61 …turn de_internal + de_kinetic + ds.Y.pressure / s.U.density - s.Y.pressure / Square(s.U.density) *… in TotalSpecificEnthalpy_fwd()
64 …PER StatePrimitive StatePrimitiveFromConservative(NewtonianIGProperties gas, StateConservative U) { in StatePrimitiveFromConservative() argument
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()
149 StateConservative U; in StateConservativeFromPrimitive() local
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()
156 return U; in StateConservativeFromPrimitive()
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()
175 …_HELPER StateEntropy StateEntropyFromConservative(NewtonianIGProperties gas, StateConservative U) { in StateEntropyFromConservative() argument
178 const CeedScalar e_kinetic = .5 * Dot3(U.momentum, U.momentum) / U.density; in StateEntropyFromConservative()
179 const CeedScalar e_internal = U.E_total - e_kinetic; in StateEntropyFromConservative()
181 const CeedScalar entropy = log(p) - gamma * log(U.density); in StateEntropyFromConservative()
184 for (int i = 0; i < 3; i++) V.S_momentum[i] = U.momentum[i] / p; 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 …const CeedScalar de_kinetic = (Dot3(s.U.momentum, dU.momentum) - e_kinetic * dU.density) / s.U.de… in StateEntropyFromConservative_fwd()
198 CeedScalar dentropy = dp / p - gamma * dU.density / s.U.density; in StateEntropyFromConservative_fwd()
202 dV.S_momentum[i] = (dU.momentum[i] - s.U.momentum[i] * dp / p) / p; in StateEntropyFromConservative_fwd()
204 dV.S_energy = -(dU.density - s.U.density * dp / p) / p; in StateEntropyFromConservative_fwd()
209 StateConservative U; in StateConservativeFromEntropy() local
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()
221 return U; in StateConservativeFromEntropy()
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()
247 StateConservative U = StateConservativeFromPrimitive(gas, Y); in StateFromPrimitive() local
249 s.U = U; in StateFromPrimitive()
257 ds.U = dU; in StateFromPrimitive_fwd()
282 …teToU(NewtonianIGProperties gas, const State input, CeedScalar U[5]) { UnpackState_U(input.U, U); } in StateToU()
323 CEED_QFUNCTION_HELPER State StateFromU(NewtonianIGProperties gas, const CeedScalar U[5]) { in StateFromU()
325 s.U.density = U[0]; in StateFromU()
326 s.U.momentum[0] = U[1]; in StateFromU()
327 s.U.momentum[1] = U[2]; in StateFromU()
328 s.U.momentum[2] = U[3]; in StateFromU()
329 s.U.E_total = U[4]; in StateFromU()
330 s.Y = StatePrimitiveFromConservative(gas, s.U); in StateFromU()
336 ds.U.density = dU[0]; in StateFromU_fwd()
337 ds.U.momentum[0] = dU[1]; in StateFromU_fwd()
338 ds.U.momentum[1] = dU[2]; in StateFromU_fwd()
339 ds.U.momentum[2] = dU[3]; in StateFromU_fwd()
340 ds.U.E_total = dU[4]; in StateFromU_fwd()
341 ds.Y = StatePrimitiveFromConservative_fwd(gas, s, ds.U); in StateFromU_fwd()
352 s.U = StateConservativeFromPrimitive(gas, s.Y); in StateFromY()
363 ds.U = StateConservativeFromPrimitive_fwd(gas, s, ds.Y); in StateFromY_fwd()
375 s.U = StateConservativeFromEntropy(gas, state_V); in StateFromV()
388 ds.U = StateConservativeFromEntropy_fwd(gas, s, state_dV); in StateFromV_fwd()
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()
435 dFlux[i].density = ds.U.momentum[i]; in FluxInviscid_fwd()
437 …dFlux[i].momentum[j] = ds.U.momentum[i] * s.Y.velocity[j] + s.U.momentum[i] * ds.Y.velocity[j] + d… in FluxInviscid_fwd()
439 …dFlux[i].E_total = (ds.U.E_total + ds.Y.pressure) * s.Y.velocity[i] + (s.U.E_total + s.Y.pressure)… in FluxInviscid_fwd()