From 547fb37f74f871a20b6a69b937986a9487a5e471 Mon Sep 17 00:00:00 2001
From: James Willis <james.s.willis@durham.ac.uk>
Date: Thu, 24 Aug 2017 15:58:59 +0100
Subject: [PATCH] Account for separations of zero.
---
src/runner_doiact_vec.c | 3 +++
1 file changed, 3 insertions(+)
diff --git a/src/runner_doiact_vec.c b/src/runner_doiact_vec.c
index 3be71c7816..e00409c157 100644
--- a/src/runner_doiact_vec.c
+++ b/src/runner_doiact_vec.c
@@ -746,6 +746,9 @@ for (int pid = 0; pid < count; pid++) {
v_hj.v = vec_load(&cell_cache->h[pjd]);
v_hj_inv = vec_reciprocal(v_hj);
+ /* To stop floating point exceptions for when particle separations are 0. */
+ v_r2.v = vec_add(v_r2.v, vec_set1(FLT_MIN));
+
runner_iact_nonsym_1_vec_force(
&v_r2, &v_dx, &v_dy, &v_dz, v_vix, v_viy, v_viz,
v_rhoi, v_grad_hi, v_pOrhoi2, v_balsara_i, v_ci,
--
GitLab