From 77e32d0e1052bfae617bfced6cce37b738b6a304 Mon Sep 17 00:00:00 2001 From: huiyuxie Date: Fri, 3 May 2024 12:11:46 -0700 Subject: [PATCH] remove TODO --- src/equations/compressible_euler_1d.jl | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/equations/compressible_euler_1d.jl b/src/equations/compressible_euler_1d.jl index 92e6b93ce1..677a661b32 100644 --- a/src/equations/compressible_euler_1d.jl +++ b/src/equations/compressible_euler_1d.jl @@ -354,8 +354,8 @@ Entropy conserving two-point flux by # Compute the necessary mean values rho_avg = 0.5f0 * (rho_ll + rho_rr) - rho_mean = ln_mean(rho_ll, rho_rr) # TODO: Modify ln_mean() from math.jl - beta_mean = ln_mean(beta_ll, beta_rr) # TODO: Modify ln_mean() from math.jl + rho_mean = ln_mean(rho_ll, rho_rr) + beta_mean = ln_mean(beta_ll, beta_rr) beta_avg = 0.5f0 * (beta_ll + beta_rr) v1_avg = 0.5f0 * (v1_ll + v1_rr) p_mean = 0.5f0 * rho_avg / beta_avg @@ -392,12 +392,12 @@ See also rho_rr, v1_rr, p_rr = cons2prim(u_rr, equations) # Compute the necessary mean values - rho_mean = ln_mean(rho_ll, rho_rr) # TODO: Modify ln_mean() from math.jl + rho_mean = ln_mean(rho_ll, rho_rr) # Algebraically equivalent to `inv_ln_mean(rho_ll / p_ll, rho_rr / p_rr)` # in exact arithmetic since # log((ϱₗ/pₗ) / (ϱᵣ/pᵣ)) / (ϱₗ/pₗ - ϱᵣ/pᵣ) # = pₗ pᵣ log((ϱₗ pᵣ) / (ϱᵣ pₗ)) / (ϱₗ pᵣ - ϱᵣ pₗ) - inv_rho_p_mean = p_ll * p_rr * inv_ln_mean(rho_ll * p_rr, rho_rr * p_ll) # TODO: Modify inv_ln_mean() from math.jl + inv_rho_p_mean = p_ll * p_rr * inv_ln_mean(rho_ll * p_rr, rho_rr * p_ll) v1_avg = 0.5f0 * (v1_ll + v1_rr) p_avg = 0.5f0 * (p_ll + p_rr) velocity_square_avg = 0.5f0 * (v1_ll * v1_rr)