Lines Matching refs:CeedScalar
29 CeedScalar left, right;
32 CEED_QFUNCTION_HELPER RoeWeights RoeSetup(CeedScalar rho_left, CeedScalar rho_right) { in RoeSetup()
33 CeedScalar sqrt_left = sqrt(rho_left), sqrt_right = sqrt(rho_right); in RoeSetup()
38 …UNCTION_HELPER RoeWeights RoeSetup_fwd(CeedScalar rho_left, CeedScalar rho_right, CeedScalar drho_… in RoeSetup_fwd()
39 CeedScalar sqrt_left = sqrt(rho_left), sqrt_right = sqrt(rho_right); in RoeSetup_fwd()
40 CeedScalar square_sum_root = Square(sqrt_left + sqrt_right); in RoeSetup_fwd()
41 …CeedScalar r_right = (sqrt_left / (2 * sqrt_right * square_sum_root)) * drho_right - (sqrt_right /… in RoeSetup_fwd()
42 …CeedScalar r_left = (sqrt_right / (2 * sqrt_left * square_sum_root)) * drho_left - (sqrt_left / (… in RoeSetup_fwd()
47 CEED_QFUNCTION_HELPER CeedScalar RoeAverage(RoeWeights r, CeedScalar q_left, CeedScalar q_right) { … in RoeAverage()
49 CEED_QFUNCTION_HELPER CeedScalar RoeAverage_fwd(RoeWeights r, RoeWeights dr, CeedScalar q_left, Cee… in RoeAverage_fwd()
50 CeedScalar dq_right) { in RoeAverage_fwd()
55 CeedScalar s_left, CeedScalar s_right) { in Flux_HLL()
56 CeedScalar U_left[5], U_right[5], F_right[5], F_left[5], F_hll[5]; in Flux_HLL()
74 … CeedScalar S_l, CeedScalar S_r, CeedScalar dS_l, CeedScalar dS_r) { in Flux_HLL_fwd()
75 CeedScalar U_l[5], U_r[5], F_r[5], F_l[5]; in Flux_HLL_fwd()
81 CeedScalar dU_l[5], dU_r[5], dF_r[5], dF_l[5], dF_hll[5] = {0.}; in Flux_HLL_fwd()
87 const CeedScalar S_diff = S_r - S_l; in Flux_HLL_fwd()
101 …LLSpeeds_Roe(NewtonianIdealGasContext gas, State left, CeedScalar u_left, State right, CeedScalar … in ComputeHLLSpeeds_Roe()
102 CeedScalar *s_left, CeedScalar *s_right) { in ComputeHLLSpeeds_Roe()
103 const CeedScalar gamma = HeatCapacityRatio(gas); in ComputeHLLSpeeds_Roe()
109 CeedScalar u_roe = RoeAverage(r, u_left, u_right); in ComputeHLLSpeeds_Roe()
111 CeedScalar H_left = TotalSpecificEnthalpy(gas, left); in ComputeHLLSpeeds_Roe()
112 CeedScalar H_right = TotalSpecificEnthalpy(gas, right); in ComputeHLLSpeeds_Roe()
113 CeedScalar H_roe = RoeAverage(r, H_left, H_right); in ComputeHLLSpeeds_Roe()
114 CeedScalar a_roe = sqrt((gamma - 1) * (H_roe - 0.5 * Square(u_roe))); in ComputeHLLSpeeds_Roe()
121 …oe_fwd(NewtonianIdealGasContext gas, State left, State dleft, CeedScalar u_left, CeedScalar du_lef… in ComputeHLLSpeeds_Roe_fwd()
122 … State right, State dright, CeedScalar u_right, CeedScalar du_right, CeedScalar *s_left, in ComputeHLLSpeeds_Roe_fwd()
123 … CeedScalar *ds_left, CeedScalar *s_right, CeedScalar *ds_right) { in ComputeHLLSpeeds_Roe_fwd()
124 const CeedScalar gamma = HeatCapacityRatio(gas); in ComputeHLLSpeeds_Roe_fwd()
131 CeedScalar u_roe = RoeAverage(r, u_left, u_right); in ComputeHLLSpeeds_Roe_fwd()
132 CeedScalar du_roe = RoeAverage_fwd(r, dr, u_left, u_right, du_left, du_right); in ComputeHLLSpeeds_Roe_fwd()
134 CeedScalar H_left = TotalSpecificEnthalpy(gas, left); in ComputeHLLSpeeds_Roe_fwd()
135 CeedScalar H_right = TotalSpecificEnthalpy(gas, right); in ComputeHLLSpeeds_Roe_fwd()
136 CeedScalar dH_left = TotalSpecificEnthalpy_fwd(gas, left, dleft); in ComputeHLLSpeeds_Roe_fwd()
137 CeedScalar dH_right = TotalSpecificEnthalpy_fwd(gas, right, dright); in ComputeHLLSpeeds_Roe_fwd()
139 CeedScalar H_roe = RoeAverage(r, H_left, H_right); in ComputeHLLSpeeds_Roe_fwd()
140 CeedScalar dH_roe = RoeAverage_fwd(r, dr, H_left, H_right, dH_left, dH_right); in ComputeHLLSpeeds_Roe_fwd()
141 CeedScalar a_roe = sqrt((gamma - 1) * (H_roe - 0.5 * Square(u_roe))); in ComputeHLLSpeeds_Roe_fwd()
142 …CeedScalar da_roe = 0.5 * sqrt((gamma - 1) / (H_roe - 0.5 * Square(u_roe))) * dH_roe; // (da/dH) … in ComputeHLLSpeeds_Roe_fwd()
163 …iemannFlux_HLL(NewtonianIdealGasContext gas, State left, State right, const CeedScalar normal[3]) { in RiemannFlux_HLL()
167 CeedScalar u_left = Dot3(left.Y.velocity, normal); in RiemannFlux_HLL()
168 CeedScalar u_right = Dot3(right.Y.velocity, normal); in RiemannFlux_HLL()
170 CeedScalar s_left, s_right; in RiemannFlux_HLL()
196 const CeedScalar normal[3]) { in RiemannFlux_HLL_fwd()
202 CeedScalar u_left = Dot3(left.Y.velocity, normal); in RiemannFlux_HLL_fwd()
203 CeedScalar u_right = Dot3(right.Y.velocity, normal); in RiemannFlux_HLL_fwd()
204 CeedScalar du_left = Dot3(dleft.Y.velocity, normal); in RiemannFlux_HLL_fwd()
205 CeedScalar du_right = Dot3(dright.Y.velocity, normal); in RiemannFlux_HLL_fwd()
207 CeedScalar s_left, ds_left, s_right, ds_right; in RiemannFlux_HLL_fwd()
220 … const CeedScalar normal[3], CeedScalar u_side, CeedScalar s_side, CeedScalar s_star) { in RiemannFlux_HLLC_Star()
221 CeedScalar fact = side.U.density * (s_side - u_side) / (s_side - s_star); in RiemannFlux_HLLC_Star()
222 CeedScalar denom = side.U.density * (s_side - u_side); in RiemannFlux_HLLC_Star()
238 … StateConservative dF_side, const CeedScalar normal[3], CeedScalar u_side, in RiemannFlux_HLLC_Star_fwd()
239 … CeedScalar du_side, CeedScalar s_side, CeedScalar ds_side, CeedScalar s_star, in RiemannFlux_HLLC_Star_fwd()
240 CeedScalar ds_star) { in RiemannFlux_HLLC_Star_fwd()
241 CeedScalar fact = side.U.density * (s_side - u_side) / (s_side - s_star); in RiemannFlux_HLLC_Star_fwd()
242 …CeedScalar dfact = (side.U.density * (ds_side - du_side) + dside.U.density * (s_side - u_side)) / … in RiemannFlux_HLLC_Star_fwd()
244 CeedScalar denom = side.U.density * (s_side - u_side); in RiemannFlux_HLLC_Star_fwd()
245 CeedScalar ddenom = side.U.density * (ds_side - du_side) + dside.U.density * (s_side - u_side); in RiemannFlux_HLLC_Star_fwd()
269 …const CeedScalar a[] = {1, ds_side * fact + s_side * dfact, s_side * fact, -ds_side, -s_sid… in RiemannFlux_HLLC_Star_fwd()
275 …emannFlux_HLLC(NewtonianIdealGasContext gas, State left, State right, const CeedScalar normal[3]) { in RiemannFlux_HLLC()
279 CeedScalar u_left = Dot3(left.Y.velocity, normal); in RiemannFlux_HLLC()
280 CeedScalar u_right = Dot3(right.Y.velocity, normal); in RiemannFlux_HLLC()
281 CeedScalar s_left, s_right; in RiemannFlux_HLLC()
285 CeedScalar rhou_left = left.U.density * u_left, rhou_right = right.U.density * u_right; in RiemannFlux_HLLC()
286 …CeedScalar numer = right.Y.pressure - left.Y.pressure + rhou_left * (s_left - u_left) - rhou_righ… in RiemannFlux_HLLC()
287 CeedScalar denom = left.U.density * (s_left - u_left) - right.U.density * (s_right - u_right); in RiemannFlux_HLLC()
288 CeedScalar s_star = numer / denom; in RiemannFlux_HLLC()
303 const CeedScalar normal[3]) { in RiemannFlux_HLLC_fwd()
309 CeedScalar u_left = Dot3(left.Y.velocity, normal); in RiemannFlux_HLLC_fwd()
310 CeedScalar u_right = Dot3(right.Y.velocity, normal); in RiemannFlux_HLLC_fwd()
311 CeedScalar du_left = Dot3(dleft.Y.velocity, normal); in RiemannFlux_HLLC_fwd()
312 CeedScalar du_right = Dot3(dright.Y.velocity, normal); in RiemannFlux_HLLC_fwd()
314 CeedScalar s_left, ds_left, s_right, ds_right; in RiemannFlux_HLLC_fwd()
318 …CeedScalar rhou_left = left.U.density * u_left, drhou_left = left.U.density * du_left + dleft.U.de… in RiemannFlux_HLLC_fwd()
319 …CeedScalar rhou_right = right.U.density * u_right, drhou_right = right.U.density * du_right + drig… in RiemannFlux_HLLC_fwd()
320 CeedScalar numer = right.Y.pressure - left.Y.pressure // in RiemannFlux_HLLC_fwd()
323 CeedScalar dnumer = dright.Y.pressure - dleft.Y.pressure // in RiemannFlux_HLLC_fwd()
326 CeedScalar denom = left.U.density * (s_left - u_left) - right.U.density * (s_right - u_right); in RiemannFlux_HLLC_fwd()
327 CeedScalar ddenom = left.U.density * (ds_left - du_left) + dleft.U.density * (s_left - u_left) // in RiemannFlux_HLLC_fwd()
329 CeedScalar s_star = numer / denom; in RiemannFlux_HLLC_fwd()
330 CeedScalar ds_star = dnumer / denom - numer / Square(denom) * ddenom; in RiemannFlux_HLLC_fwd()