From 5c08cf4a4a7407d3ffa90182baafaf8c286b8c33 Mon Sep 17 00:00:00 2001
From: Daniel Doehring <daniel.doehring@rwth-aachen.de>
Date: Fri, 24 Nov 2023 10:45:07 +0100
Subject: [PATCH] Switch (back) to Einfeldt wave-speed estimate from paper
 (#1751)

* Switch (back) to wave-speed estimate from paper

* upadte test vals
---
 src/equations/ideal_glm_mhd_1d.jl |  2 +-
 test/test_tree_1d_mhd.jl          | 15 ++++++++-------
 2 files changed, 9 insertions(+), 8 deletions(-)

diff --git a/src/equations/ideal_glm_mhd_1d.jl b/src/equations/ideal_glm_mhd_1d.jl
index a465571989b..eca06c8c203 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 447572eee88..8895fe30e8b 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)