Switch (back) to Einfeldt wave-speed estimate from paper (#1751)

* Switch (back) to wave-speed estimate from paper * upadte test vals
trixi-framework · Nov 24, 2023 · 5c08cf4 · 5c08cf4
1 parent 004bc84
commit 5c08cf4
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 8 deletions.
diff --git 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
@@ -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)