Lines Matching refs:s
24 CEED_QFUNCTION_HELPER void UnpackState_U(StateConservative s, CeedScalar U[5]) { in UnpackState_U() argument
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()
30 CEED_QFUNCTION_HELPER void UnpackState_Y(StatePrimitive s, CeedScalar Y[5]) { in UnpackState_Y() argument
31 Y[0] = s.pressure; in UnpackState_Y()
32 for (int i = 0; i < 3; i++) Y[i + 1] = s.velocity[i]; in UnpackState_Y()
33 Y[4] = s.temperature; in UnpackState_Y()
36 CEED_QFUNCTION_HELPER void UnpackState_V(StateEntropy s, CeedScalar V[5]) { in UnpackState_V() argument
37 V[0] = s.S_density; in UnpackState_V()
38 for (int i = 0; i < 3; i++) V[i + 1] = s.S_momentum[i]; in UnpackState_V()
39 V[4] = s.S_energy; in UnpackState_V()
52 CEED_QFUNCTION_HELPER CeedScalar TotalSpecificEnthalpy(NewtonianIGProperties gas, const State s) { in TotalSpecificEnthalpy() argument
53 CeedScalar e_kinetic = 0.5 * Dot3(s.Y.velocity, s.Y.velocity); in TotalSpecificEnthalpy()
54 CeedScalar e_internal = gas.cv * s.Y.temperature; in TotalSpecificEnthalpy()
55 return e_internal + e_kinetic + s.Y.pressure / s.U.density; in TotalSpecificEnthalpy()
58 CEED_QFUNCTION_HELPER CeedScalar TotalSpecificEnthalpy_fwd(NewtonianIGProperties gas, const State s… in TotalSpecificEnthalpy_fwd() argument
59 CeedScalar de_kinetic = Dot3(ds.Y.velocity, s.Y.velocity); in TotalSpecificEnthalpy_fwd()
61 …return de_internal + de_kinetic + ds.Y.pressure / s.U.density - s.Y.pressure / Square(s.U.density)… in TotalSpecificEnthalpy_fwd()
75 …imitive StatePrimitiveFromConservative_fwd(NewtonianIGProperties gas, State s, StateConservative d… in StatePrimitiveFromConservative_fwd() argument
78 dY.velocity[i] = (dU.momentum[i] - s.Y.velocity[i] * dU.density) / s.U.density; in StatePrimitiveFromConservative_fwd()
80 CeedScalar e_kinetic = .5 * Dot3(s.Y.velocity, s.Y.velocity); in StatePrimitiveFromConservative_fwd()
81 CeedScalar de_kinetic = Dot3(dY.velocity, s.Y.velocity); 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()
105 CEED_QFUNCTION_HELPER StateEntropy StateEntropyFromPrimitive_fwd(NewtonianIGProperties gas, State s… in StateEntropyFromPrimitive_fwd() argument
108 …r drho = (dY.pressure * s.Y.temperature - s.Y.pressure * dY.temperature) / (GasConstant(gas… in StateEntropyFromPrimitive_fwd()
110 const CeedScalar e_kinetic = .5 * Dot3(s.Y.velocity, s.Y.velocity); in StateEntropyFromPrimitive_fwd()
111 const CeedScalar de_kinetic = Dot3(dY.velocity, s.Y.velocity); in StateEntropyFromPrimitive_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()
118 …for (CeedInt i = 0; i < 3; i++) dV.S_momentum[i] = rho_div_p * dY.velocity[i] + drho_div_p * s.Y.v… in StateEntropyFromPrimitive_fwd()
135 …StatePrimitive StatePrimitiveFromEntropy_fwd(NewtonianIGProperties gas, State s, StateEntropy dV) { in StatePrimitiveFromEntropy_fwd() argument
137 StateEntropy V = StateEntropyFromPrimitive(gas, s.Y); in StatePrimitiveFromEntropy_fwd()
141 const CeedScalar e_kinetic = 0.5 * Dot3(s.Y.velocity, s.Y.velocity); in StatePrimitiveFromEntropy_fwd()
142 const CeedScalar de_kinetic = Dot3(dY.velocity, s.Y.velocity); in StatePrimitiveFromEntropy_fwd()
144 …dY.pressure = s.Y.pressure * (-dentropy - gamma * dV.S_energy / V.S_energy) / (gam… in StatePrimitiveFromEntropy_fwd()
159 …rvative StateConservativeFromPrimitive_fwd(NewtonianIGProperties gas, State s, StatePrimitive dY) { in StateConservativeFromPrimitive_fwd() argument
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()
165 CeedScalar e_kinetic = .5 * Dot3(s.Y.velocity, s.Y.velocity); in StateConservativeFromPrimitive_fwd()
166 CeedScalar de_kinetic = Dot3(dY.velocity, s.Y.velocity); in StateConservativeFromPrimitive_fwd()
167 CeedScalar e_internal = gas.cv * s.Y.temperature; in StateConservativeFromPrimitive_fwd()
171 dU.E_total = dU.density * e_total + s.U.density * de_total; in StateConservativeFromPrimitive_fwd()
189 …teEntropy StateEntropyFromConservative_fwd(NewtonianIGProperties gas, State s, StateConservative d… in StateEntropyFromConservative_fwd() argument
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()
195 const CeedScalar p = s.Y.pressure; 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()
224 …onservative StateConservativeFromEntropy_fwd(NewtonianIGProperties gas, State s, StateEntropy dV) { in StateConservativeFromEntropy_fwd() argument
227 StateEntropy V = StateEntropyFromPrimitive(gas, s.Y); in StateConservativeFromEntropy_fwd()
228 …for (int i = 0; i < 3; i++) dvelocity[i] = (-dV.S_momentum[i] - s.Y.velocity[i] * dV.S_energy) / V… in StateConservativeFromEntropy_fwd()
230 const CeedScalar e_kinetic = 0.5 * Dot3(s.Y.velocity, s.Y.velocity); in StateConservativeFromEntropy_fwd()
231 const CeedScalar de_kinetic = Dot3(dvelocity, s.Y.velocity); 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()
248 State s; in StateFromPrimitive() local
249 s.U = U; in StateFromPrimitive()
250 s.Y = Y; in StateFromPrimitive()
251 return s; in StateFromPrimitive()
254 CEED_QFUNCTION_HELPER State StateFromPrimitive_fwd(NewtonianIGProperties gas, State s, StatePrimiti… in StateFromPrimitive_fwd() argument
255 StateConservative dU = StateConservativeFromPrimitive_fwd(gas, s, dY); in StateFromPrimitive_fwd()
324 State s; in StateFromU() local
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()
331 return s; in StateFromU()
334 CEED_QFUNCTION_HELPER State StateFromU_fwd(NewtonianIGProperties gas, State s, const CeedScalar dU[… in StateFromU_fwd() argument
341 ds.Y = StatePrimitiveFromConservative_fwd(gas, s, ds.U); in StateFromU_fwd()
346 State s; in StateFromY() local
347 s.Y.pressure = Y[0]; in StateFromY()
348 s.Y.velocity[0] = Y[1]; in StateFromY()
349 s.Y.velocity[1] = Y[2]; in StateFromY()
350 s.Y.velocity[2] = Y[3]; in StateFromY()
351 s.Y.temperature = Y[4]; in StateFromY()
352 s.U = StateConservativeFromPrimitive(gas, s.Y); in StateFromY()
353 return s; in StateFromY()
356 CEED_QFUNCTION_HELPER State StateFromY_fwd(NewtonianIGProperties gas, State s, const CeedScalar dY[… in StateFromY_fwd() argument
363 ds.U = StateConservativeFromPrimitive_fwd(gas, s, ds.Y); in StateFromY_fwd()
368 State s; in StateFromV() local
375 s.U = StateConservativeFromEntropy(gas, state_V); in StateFromV()
376 s.Y = StatePrimitiveFromEntropy(gas, state_V); in StateFromV()
377 return s; in StateFromV()
380 CEED_QFUNCTION_HELPER State StateFromV_fwd(NewtonianIGProperties gas, State s, const CeedScalar dV[… in StateFromV_fwd() argument
388 ds.U = StateConservativeFromEntropy_fwd(gas, s, state_dV); in StateFromV_fwd()
389 ds.Y = StatePrimitiveFromEntropy_fwd(gas, s, state_dV); in StateFromV_fwd()
394 State s; in StateFromQ() local
397 s = StateFromU(gas, Q); in StateFromQ()
400 s = StateFromY(gas, Q); in StateFromQ()
403 s = StateFromV(gas, Q); in StateFromQ()
406 return s; in StateFromQ()
409 CEED_QFUNCTION_HELPER State StateFromQ_fwd(NewtonianIGProperties gas, State s, const CeedScalar dQ[… in StateFromQ_fwd() argument
413 ds = StateFromU_fwd(gas, s, dQ); in StateFromQ_fwd()
416 ds = StateFromY_fwd(gas, s, dQ); in StateFromQ_fwd()
419 ds = StateFromV_fwd(gas, s, dQ); in StateFromQ_fwd()
425 CEED_QFUNCTION_HELPER void FluxInviscid(const NewtonianIGProperties gas, const State s, StateConser… in FluxInviscid() argument
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()
433 CEED_QFUNCTION_HELPER void FluxInviscid_fwd(const NewtonianIGProperties gas, const State s, const S… in FluxInviscid_fwd() argument
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()
443 …ervative FluxInviscidDotNormal(const NewtonianIGProperties gas, const State s, const CeedScalar no… in FluxInviscidDotNormal() argument
445 FluxInviscid(gas, s, Flux); in FluxInviscidDotNormal()
454 …servative FluxInviscidDotNormal_fwd(const NewtonianIGProperties gas, const State s, const State ds, in FluxInviscidDotNormal_fwd() argument
457 FluxInviscid_fwd(gas, s, ds, dFlux); in FluxInviscidDotNormal_fwd()
466 CEED_QFUNCTION_HELPER void FluxInviscidStrong(const NewtonianIGProperties gas, const State s, const… in FluxInviscidStrong() argument
470 FluxInviscid_fwd(gas, s, ds[i], dF); in FluxInviscidStrong()
581 …radientFromReference(CeedInt Q, CeedInt i, NewtonianIGProperties gas, State s, StateVariable state… in StatePhysicalGradientFromReference() argument
591 grad_s[j] = StateFromQ_fwd(gas, s, dqi, state_var); in StatePhysicalGradientFromReference()
595 …atePhysicalGradientFromReference_Boundary(CeedInt Q, CeedInt i, NewtonianIGProperties gas, State s, in StatePhysicalGradientFromReference_Boundary() argument
603 grad_s[k] = StateFromQ_fwd(gas, s, dqi, state_var); in StatePhysicalGradientFromReference_Boundary()