Lines Matching refs:U
27 StateConservative U; member
31 CEED_QFUNCTION_HELPER void UnpackState_U(StateConservative s, CeedScalar U[5]) { in UnpackState_U()
32 U[0] = s.density; in UnpackState_U()
33 for (int i = 0; i < 3; i++) U[i + 1] = s.momentum[i]; in UnpackState_U()
34 U[4] = s.E_total; in UnpackState_U()
62 return e_internal + e_kinetic + s.Y.pressure / s.U.density; in TotalSpecificEnthalpy()
68 …turn de_internal + de_kinetic + ds.Y.pressure / s.U.density - s.Y.pressure / Square(s.U.density) *… in TotalSpecificEnthalpy_fwd()
71 … StatePrimitive StatePrimitiveFromConservative(NewtonianIdealGasContext gas, StateConservative U) { in StatePrimitiveFromConservative() argument
73 for (CeedInt i = 0; i < 3; i++) Y.velocity[i] = U.momentum[i] / U.density; in StatePrimitiveFromConservative()
75 CeedScalar e_total = U.E_total / U.density; in StatePrimitiveFromConservative()
78 Y.pressure = (HeatCapacityRatio(gas) - 1) * U.density * e_internal; in StatePrimitiveFromConservative()
85 dY.velocity[i] = (dU.momentum[i] - s.Y.velocity[i] * dU.density) / s.U.density; in StatePrimitiveFromConservative_fwd()
89 CeedScalar e_total = s.U.E_total / s.U.density; in StatePrimitiveFromConservative_fwd()
90 CeedScalar de_total = (dU.E_total - e_total * dU.density) / s.U.density; in StatePrimitiveFromConservative_fwd()
94 …dY.pressure = (HeatCapacityRatio(gas) - 1) * (dU.density * e_internal + s.U.density * d… in StatePrimitiveFromConservative_fwd()
119 const CeedScalar rho_div_p = s.U.density / s.Y.pressure; in StateEntropyFromPrimitive_fwd()
120 …const CeedScalar drho_div_p = (drho * s.Y.pressure - s.U.density * dY.pressure) / Square(s.Y.press… in StateEntropyFromPrimitive_fwd()
122 CeedScalar dentropy = dY.pressure / s.Y.pressure - gamma * drho / s.U.density; in StateEntropyFromPrimitive_fwd()
156 StateConservative U; in StateConservativeFromPrimitive() local
157 U.density = Y.pressure / (GasConstant(gas) * Y.temperature); in StateConservativeFromPrimitive()
158 for (int i = 0; i < 3; i++) U.momentum[i] = U.density * Y.velocity[i]; in StateConservativeFromPrimitive()
162 U.E_total = U.density * e_total; in StateConservativeFromPrimitive()
163 return U; in StateConservativeFromPrimitive()
170 dU.momentum[i] = dU.density * s.Y.velocity[i] + s.U.density * dY.velocity[i]; in StateConservativeFromPrimitive_fwd()
178 dU.E_total = dU.density * e_total + s.U.density * de_total; in StateConservativeFromPrimitive_fwd()
182 …LPER StateEntropy StateEntropyFromConservative(NewtonianIdealGasContext gas, StateConservative U) { in StateEntropyFromConservative() argument
185 const CeedScalar e_kinetic = .5 * Dot3(U.momentum, U.momentum) / U.density; in StateEntropyFromConservative()
186 const CeedScalar e_internal = U.E_total - e_kinetic; in StateEntropyFromConservative()
188 const CeedScalar entropy = log(p) - gamma * log(U.density); 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()
199 const CeedScalar e_kinetic = .5 * Dot3(s.U.momentum, s.U.momentum) / s.U.density; in StateEntropyFromConservative_fwd()
200 …const CeedScalar de_kinetic = (Dot3(s.U.momentum, dU.momentum) - e_kinetic * dU.density) / s.U.de… in StateEntropyFromConservative_fwd()
205 CeedScalar dentropy = dp / p - gamma * dU.density / s.U.density; in StateEntropyFromConservative_fwd()
209 dV.S_momentum[i] = (dU.momentum[i] - s.U.momentum[i] * dp / p) / p; in StateEntropyFromConservative_fwd()
211 dV.S_energy = -(dU.density - s.U.density * dp / p) / p; in StateEntropyFromConservative_fwd()
216 StateConservative U; in StateConservativeFromEntropy() local
223 U.density = exp(log_rho); in StateConservativeFromEntropy()
224 for (int i = 0; i < 3; i++) U.momentum[i] = U.density * velocity[i]; in StateConservativeFromEntropy()
227 U.E_total = U.density * (e_internal + e_kinetic); in StateConservativeFromEntropy()
228 return U; in StateConservativeFromEntropy()
244 …for (int i = 0; i < 3; i++) dU.momentum[i] = dU.density * s.Y.velocity[i] + s.U.density * dvelocit… in StateConservativeFromEntropy_fwd()
249 dU.E_total = dU.density * e_total + s.U.density * (de_internal + de_kinetic); in StateConservativeFromEntropy_fwd()
254 StateConservative U = StateConservativeFromPrimitive(gas, Y); in StateFromPrimitive() local
256 s.U = U; in StateFromPrimitive()
264 ds.U = dU; in StateFromPrimitive_fwd()
289 …oU(NewtonianIdealGasContext gas, const State input, CeedScalar U[5]) { UnpackState_U(input.U, U); } in StateToU()
328 CEED_QFUNCTION_HELPER State StateFromU(NewtonianIdealGasContext gas, const CeedScalar U[5]) { in StateFromU()
330 s.U.density = U[0]; in StateFromU()
331 s.U.momentum[0] = U[1]; in StateFromU()
332 s.U.momentum[1] = U[2]; in StateFromU()
333 s.U.momentum[2] = U[3]; in StateFromU()
334 s.U.E_total = U[4]; in StateFromU()
335 s.Y = StatePrimitiveFromConservative(gas, s.U); in StateFromU()
341 ds.U.density = dU[0]; in StateFromU_fwd()
342 ds.U.momentum[0] = dU[1]; in StateFromU_fwd()
343 ds.U.momentum[1] = dU[2]; in StateFromU_fwd()
344 ds.U.momentum[2] = dU[3]; in StateFromU_fwd()
345 ds.U.E_total = dU[4]; in StateFromU_fwd()
346 ds.Y = StatePrimitiveFromConservative_fwd(gas, s, ds.U); in StateFromU_fwd()
357 s.U = StateConservativeFromPrimitive(gas, s.Y); in StateFromY()
368 ds.U = StateConservativeFromPrimitive_fwd(gas, s, ds.Y); in StateFromY_fwd()
380 s.U = StateConservativeFromEntropy(gas, state_V); in StateFromV()
393 ds.U = StateConservativeFromEntropy_fwd(gas, s, state_dV); in StateFromV_fwd()
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()
440 dFlux[i].density = ds.U.momentum[i]; in FluxInviscid_fwd()
442 …dFlux[i].momentum[j] = ds.U.momentum[i] * s.Y.velocity[j] + s.U.momentum[i] * ds.Y.velocity[j] + d… in FluxInviscid_fwd()
444 …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()