Forgot to use lengh() instead of dlength_vec().

2019-07-30 11:50:20 +08:00
parent 0acf52d1b4
commit a6dcee2381
1 changed files with 1 additions and 1 deletions
--- 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