From 0e0ace397056a60af507de13ac3a468c7a823a27 Mon Sep 17 00:00:00 2001 From: jpekkila Date: Tue, 6 Aug 2019 18:07:29 +0300 Subject: [PATCH] Pure hydro now works with autotests --- src/standalone/model/model_rk3.cc | 6 ++---- 1 file changed, 2 insertions(+), 4 deletions(-) diff --git a/src/standalone/model/model_rk3.cc b/src/standalone/model/model_rk3.cc index 3907f93..e033aa1 100644 --- a/src/standalone/model/model_rk3.cc +++ b/src/standalone/model/model_rk3.cc @@ -560,11 +560,9 @@ momentum(const ModelVectorData& uu, const ModelScalarData& lnrho #else // !!!!!!!!!!!!!!!!%JP: NOTE TODO IMPORTANT!!!!!!!!!!!!!!!!!!!!!!!! // NOT CHECKED FOR CORRECTNESS: USE AT YOUR OWN RISK - const ModelMatrix S = stress_tensor(uu); - const ModelScalar cs2 = get(AC_cs2_sound) * - expl((get(AC_gamma) - 1) * (value(lnrho) - get(AC_lnrho0))); + const ModelMatrix S = stress_tensor(uu); - const ModelVector mom = -mul(gradients(uu), value(uu)) - cs2 * gradient(lnrho) + + const ModelVector mom = -mul(gradients(uu), value(uu)) - get(AC_cs2_sound) * gradient(lnrho) + get(AC_nu_visc) * (laplace_vec(uu) + ModelScalar(1. / 3.) * gradient_of_divergence(uu) + ModelScalar(2.) * mul(S, gradient(lnrho))) +