fixup in quick exit enthalpy

sxs-collaboration · Oct 10, 2023 · c709337 · c709337
1 parent ba2a32a
commit c709337
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/Evolution/Systems/GrMhd/ValenciaDivClean/PrimitiveFromConservative.cpp b/src/Evolution/Systems/GrMhd/ValenciaDivClean/PrimitiveFromConservative.cpp
@@ -133,7 +133,7 @@ bool PrimitiveFromConservative<OrderedListOfPrimitiveRecoverySchemes,
                 Scalar<double>{specific_energy_at_point},
                 Scalar<double>{get(*electron_fraction)[s]}));
         const double enthalpy_density_at_point =
-            floorD + specific_energy_at_point + pressure_at_point;
+            floorD + specific_energy_at_point * floorD + pressure_at_point;
         primitive_data = PrimitiveRecoverySchemes::PrimitiveRecoveryData{
             floorD,
             1.0,
@@ -151,7 +151,7 @@ bool PrimitiveFromConservative<OrderedListOfPrimitiveRecoverySchemes,
                 Scalar<double>{floorD},
                 Scalar<double>{specific_energy_at_point}));
         const double enthalpy_density_at_point =
-            floorD + specific_energy_at_point + pressure_at_point;
+            floorD + specific_energy_at_point * floorD + pressure_at_point;
         primitive_data = PrimitiveRecoverySchemes::PrimitiveRecoveryData{
             floorD,
             1.0,
@@ -166,7 +166,7 @@ bool PrimitiveFromConservative<OrderedListOfPrimitiveRecoverySchemes,
         const double pressure_at_point = get(
             equation_of_state.pressure_from_density(Scalar<double>{floorD}));
         const double enthalpy_density_at_point =
-            floorD + specific_energy_at_point + pressure_at_point;
+            floorD + specific_energy_at_point * floorD + pressure_at_point;
         primitive_data = PrimitiveRecoverySchemes::PrimitiveRecoveryData{
             floorD,
             1.0,