diff --git a/src/equations/ideal_glm_mhd_1d.jl b/src/equations/ideal_glm_mhd_1d.jl index a465571989..eca06c8c20 100644 --- a/src/equations/ideal_glm_mhd_1d.jl +++ b/src/equations/ideal_glm_mhd_1d.jl @@ -285,7 +285,7 @@ function flux_hllc(u_ll, u_rr, orientation::Integer, f_ll = flux(u_ll, orientation, equations) f_rr = flux(u_rr, orientation, equations) - SsL, SsR = min_max_speed_naive(u_ll, u_rr, orientation, equations) + SsL, SsR = min_max_speed_einfeldt(u_ll, u_rr, orientation, equations) sMu_L = SsL - v1_ll sMu_R = SsR - v1_rr if SsL >= 0 diff --git a/test/test_tree_1d_mhd.jl b/test/test_tree_1d_mhd.jl index 447572eee8..8895fe30e8 100644 --- a/test/test_tree_1d_mhd.jl +++ b/test/test_tree_1d_mhd.jl @@ -210,15 +210,16 @@ end @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_mhd_torrilhon_shock_tube.jl"), surface_flux=flux_hllc, l2=[ - 0.4574266553239646, 0.4794143154876439, 0.3407079689595056, - 0.44797768430829343, 0.9206916204424165, - 1.3216517820475193e-16, 0.2889748702415378, - 0.25529778018020927, + 0.45738965718253993, 0.479402222862685, 0.34069729746967664, + 0.44795514335568865, 0.9206813325913135, + 1.3216517820475193e-16, 0.2889672868491632, + 0.2552794220777942, ], linf=[ - 1.217943947570543, 0.8868438459815245, 0.878215340656725, - 0.9710882819266371, 1.6742759645320984, - 2.220446049250313e-16, 0.704710220504591, 0.6562122176458641, + 1.2181099854251536, 0.8869319941747589, 0.8763562906332134, + 0.9712221036087284, 1.6734231113527818, + 2.220446049250313e-16, 0.7035011427822779, + 0.6562884129650286, ]) # Ensure that we do not have excessive memory allocations # (e.g., from type instabilities)