From a6dcee2381679def5471838191df2a70d4b431ce Mon Sep 17 00:00:00 2001 From: JackHsu Date: Tue, 30 Jul 2019 11:50:20 +0800 Subject: [PATCH] Forgot to use lengh() instead of dlength_vec(). --- acc/mhd_solver/stencil_process.sps | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/acc/mhd_solver/stencil_process.sps b/acc/mhd_solver/stencil_process.sps index 2cd6633..4b8bae9 100644 --- a/acc/mhd_solver/stencil_process.sps +++ b/acc/mhd_solver/stencil_process.sps @@ -72,7 +72,7 @@ gravity() { const Vector sink_pos = (Vector){DCONST_REAL(AC_sink_pos_x), DCONST_REAL(AC_sink_pos_y), DCONST_REAL(AC_sink_pos_z)}; - const Scalar distance = length_vec(grid_pos - sink_pos); + const Scalar distance = length(grid_pos - sink_pos); force_gravity = (6.67e-8 * sink_mass) / (distance * distance); return force_gravity; }//first attempt to do a self-containing LSINK module