Add numerical support of other real types (Float32)

src/Trixi.jl

@@ -216,10 +216,11 @@ export initial_condition_eoc_test_coupled_euler_gravity,
        source_terms_eoc_test_coupled_euler_gravity, source_terms_eoc_test_euler
 
 export cons2cons, cons2prim, prim2cons, cons2macroscopic, cons2state, cons2mean,
-       cons2entropy, entropy2cons
+       cons2entropy, entropy2cons, cons2entropy_guermond_etal
 export density, pressure, density_pressure, velocity, global_mean_vars,
        equilibrium_distribution, waterheight_pressure
-export entropy, energy_total, energy_kinetic, energy_internal, energy_magnetic,
+export entropy, entropy_math, entropy_thermodynamic, entropy_guermond_etal,
+       energy_total, energy_kinetic, energy_internal, energy_magnetic,
        cross_helicity,
        enstrophy
 export lake_at_rest_error

src/equations/compressible_euler_3d.jl

@@ -418,17 +418,18 @@ end
     return SVector(f1, f2, f3, f4, f5)
 end
 
-@inline function flux(u, normal::AbstractVector,
+@inline function flux(u, normal_direction::AbstractVector,
                       equations::CompressibleEulerEquations3D)
     rho_e = last(u)
     rho, v1, v2, v3, p = cons2prim(u, equations)
 
-    v_normal = v1 * normal[1] + v2 * normal[2] + v3 * normal[3]
+    v_normal = v1 * normal_direction[1] + v2 * normal_direction[2] +
+               v3 * normal_direction[3]
     rho_v_normal = rho * v_normal
     f1 = rho_v_normal
-    f2 = rho_v_normal * v1 + p * normal[1]
-    f3 = rho_v_normal * v2 + p * normal[2]
-    f4 = rho_v_normal * v3 + p * normal[3]
+    f2 = rho_v_normal * v1 + p * normal_direction[1]
+    f3 = rho_v_normal * v2 + p * normal_direction[2]
+    f4 = rho_v_normal * v3 + p * normal_direction[3]
     f5 = (rho_e + p) * v_normal
     return SVector(f1, f2, f3, f4, f5)
 end