diff --git a/src/equations/compressible_euler_quasi_1d.jl b/src/equations/compressible_euler_quasi_1d.jl index 6656b4cbece..7e4b0228cff 100644 --- a/src/equations/compressible_euler_quasi_1d.jl +++ b/src/equations/compressible_euler_quasi_1d.jl @@ -199,7 +199,7 @@ end # Calculate 1D flux for a single point @inline function flux(u, orientation::Integer, - equations::CompressibleEulerEquationsQuasi1D) + equations::CompressibleEulerEquationsQuasi1D) a_rho, a_rho_v1, a_e, a = u rho, v1, p, a = cons2prim(u, equations) e = a_e / a @@ -235,7 +235,7 @@ Further details are available in the paper: # For flux differencing using non-conservative terms, we return the # non-conservative flux scaled by 2. This cancels with a factor of 0.5 # in the arithmetic average of {p}. - p_avg = p_ll + p_rr + p_avg = p_ll + p_rr z = zero(eltype(u_ll)) diff --git a/test/test_unit.jl b/test/test_unit.jl index 8276ee87b53..b3ed29d38e3 100644 --- a/test/test_unit.jl +++ b/test/test_unit.jl @@ -662,8 +662,10 @@ end orientations = [1] for orientation in orientations - @test flux_chan_etal(u, u, orientation, equations) ≈ flux(u, orientation, equations) + @test flux_chan_etal(u, u, orientation, equations) ≈ + flux(u, orientation, equations) end +end @timed_testset "Consistency check for HLL flux (naive): LEE" begin flux_hll = FluxHLL(min_max_speed_naive)