17 static PetscErrorCode RiemannSolverUnitTests(NewtonianIGProperties gas, CeedScalar rtol);
141   State S_infty = StateFromPrimitive(newtonian_ig_ctx->gas, Y_inf);  in FreestreamBCSetup()
172     if (run_unit_tests) PetscCall(RiemannSolverUnitTests(newtonian_ig_ctx->gas, 5e-7));  in FreestreamBCSetup()
217 static PetscErrorCode TestRiemannHLL_fwd(NewtonianIGProperties gas, CeedScalar rtol_0, CeedScalar r…  in TestRiemannHLL_fwd()  argument
222   const CeedScalar p           = (HeatCapacityRatio(gas) - 1) * rho * gas.cv * T;  in TestRiemannHLL_fwd()
230   State left0  = StateFromY(gas, Y0_left);  in TestRiemannHLL_fwd()
231   State right0 = StateFromY(gas, Y0_right);  in TestRiemannHLL_fwd()
245       State dleft0  = StateFromY_fwd(gas, left0, dY_left);  in TestRiemannHLL_fwd()
246       State dright0 = StateFromY_fwd(gas, right0, dY_right);  in TestRiemannHLL_fwd()
248 …  StateConservative dFlux_state = RiemannFlux_HLL_fwd(gas, left0, dleft0, right0, dright0, normal);  in TestRiemannHLL_fwd()
262       State left1  = StateFromY(gas, Y1_left);  in TestRiemannHLL_fwd()
263       State right1 = StateFromY(gas, Y1_right);  in TestRiemannHLL_fwd()
265       StateConservative Flux0_state = RiemannFlux_HLL(gas, left0, right0, normal);  in TestRiemannHLL_fwd()
266       StateConservative Flux1_state = RiemannFlux_HLL(gas, left1, right1, normal);  in TestRiemannHLL_fwd()
279 static PetscErrorCode TestRiemannHLLC_fwd(NewtonianIGProperties gas, CeedScalar rtol_0, CeedScalar …  in TestRiemannHLLC_fwd()  argument
284   const CeedScalar p           = (HeatCapacityRatio(gas) - 1) * rho * gas.cv * T;  in TestRiemannHLLC_fwd()
292   State left0  = StateFromY(gas, Y0_left);  in TestRiemannHLLC_fwd()
293   State right0 = StateFromY(gas, Y0_right);  in TestRiemannHLLC_fwd()
307       State dleft0  = StateFromY_fwd(gas, left0, dY_left);  in TestRiemannHLLC_fwd()
308       State dright0 = StateFromY_fwd(gas, right0, dY_right);  in TestRiemannHLLC_fwd()
310 … StateConservative dFlux_state = RiemannFlux_HLLC_fwd(gas, left0, dleft0, right0, dright0, normal);  in TestRiemannHLLC_fwd()
324       State left1  = StateFromY(gas, Y1_left);  in TestRiemannHLLC_fwd()
325       State right1 = StateFromY(gas, Y1_right);  in TestRiemannHLLC_fwd()
327       StateConservative Flux0_state = RiemannFlux_HLLC(gas, left0, right0, normal);  in TestRiemannHLLC_fwd()
328       StateConservative Flux1_state = RiemannFlux_HLLC(gas, left1, right1, normal);  in TestRiemannHLLC_fwd()
341 static PetscErrorCode TestComputeHLLSpeeds_Roe_fwd(NewtonianIGProperties gas, CeedScalar rtol) {  in TestComputeHLLSpeeds_Roe_fwd()  argument
346   const CeedScalar p           = (HeatCapacityRatio(gas) - 1) * rho * gas.cv * T;  in TestComputeHLLSpeeds_Roe_fwd()
354   State left0  = StateFromY(gas, Y0_left);  in TestComputeHLLSpeeds_Roe_fwd()
355   State right0 = StateFromY(gas, Y0_right);  in TestComputeHLLSpeeds_Roe_fwd()
371       State      dleft0   = StateFromY_fwd(gas, left0, dY_left);  in TestComputeHLLSpeeds_Roe_fwd()
372       State      dright0  = StateFromY_fwd(gas, right0, dY_right);  in TestComputeHLLSpeeds_Roe_fwd()
377 …ComputeHLLSpeeds_Roe_fwd(gas, left0, dleft0, u_left0, du_left, right0, dright0, u_right0, du_right…  in TestComputeHLLSpeeds_Roe_fwd()
389       State      left1    = StateFromY(gas, Y1_left);  in TestComputeHLLSpeeds_Roe_fwd()
390       State      right1   = StateFromY(gas, Y1_right);  in TestComputeHLLSpeeds_Roe_fwd()
395       ComputeHLLSpeeds_Roe(gas, left0, u_left0, right0, u_right0, &s_left0, &s_right0);  in TestComputeHLLSpeeds_Roe_fwd()
396       ComputeHLLSpeeds_Roe(gas, left1, u_left1, right1, u_right1, &s_left1, &s_right1);  in TestComputeHLLSpeeds_Roe_fwd()
416 static PetscErrorCode TestTotalSpecificEnthalpy_fwd(NewtonianIGProperties gas, CeedScalar rtol) {  in TestTotalSpecificEnthalpy_fwd()  argument
421   const CeedScalar p      = (HeatCapacityRatio(gas) - 1) * rho * gas.cv * T;  in TestTotalSpecificEnthalpy_fwd()
427   State state0 = StateFromY(gas, Y0);  in TestTotalSpecificEnthalpy_fwd()
435       State dstate0 = StateFromY_fwd(gas, state0, dY);  in TestTotalSpecificEnthalpy_fwd()
436       dH            = TotalSpecificEnthalpy_fwd(gas, state0, dstate0);  in TestTotalSpecificEnthalpy_fwd()
443       State state1 = StateFromY(gas, Y1);  in TestTotalSpecificEnthalpy_fwd()
445       H0    = TotalSpecificEnthalpy(gas, state0);  in TestTotalSpecificEnthalpy_fwd()
446       H1    = TotalSpecificEnthalpy(gas, state1);  in TestTotalSpecificEnthalpy_fwd()
463 static PetscErrorCode TestRowSetup_fwd(NewtonianIGProperties gas, CeedScalar rtol) {  in TestRowSetup_fwd()  argument
504 static PetscErrorCode RiemannSolverUnitTests(NewtonianIGProperties gas, CeedScalar rtol) {  in RiemannSolverUnitTests()  argument
506   PetscCall(TestRiemannHLL_fwd(gas, rtol, rtol, rtol));  in RiemannSolverUnitTests()
507   PetscCall(TestRiemannHLLC_fwd(gas, rtol, rtol, rtol));  in RiemannSolverUnitTests()
508   PetscCall(TestComputeHLLSpeeds_Roe_fwd(gas, rtol));  in RiemannSolverUnitTests()
509   PetscCall(TestTotalSpecificEnthalpy_fwd(gas, rtol));  in RiemannSolverUnitTests()
510   PetscCall(TestRowSetup_fwd(gas, rtol));  in RiemannSolverUnitTests()