From 98713ff9d2ade911472e8a24315eb9334312f1c8 Mon Sep 17 00:00:00 2001 From: Miikka Vaisala Date: Wed, 3 Jul 2019 14:49:10 +0800 Subject: [PATCH] A possible bug note added. Will look into late. --- acc/mhd_solver/stencil_process.sps | 1 + 1 file changed, 1 insertion(+) diff --git a/acc/mhd_solver/stencil_process.sps b/acc/mhd_solver/stencil_process.sps index 4960554..8fbd8ae 100644 --- a/acc/mhd_solver/stencil_process.sps +++ b/acc/mhd_solver/stencil_process.sps @@ -69,6 +69,7 @@ momentum(in Vector uu, in Scalar lnrho, in Scalar ss, in Vector aa) { const Scalar cs2 = cs2_sound * exp(gamma * value(ss) / cp_sound + (gamma - 1) * (value(lnrho) - LNRHO0)); const Vector j = (Scalar(1.) / mu0) * (gradient_of_divergence(aa) - laplace_vec(aa)); // Current density const Vector B = curl(aa); + //TODO: DOES INTHERMAL VERSTION INCLUDE THE MAGNETIC FIELD? const Scalar inv_rho = Scalar(1.) / exp(value(lnrho)); // Regex replace CPU constants with get\(AC_([a-zA-Z_0-9]*)\)