From 8025d9809389a6fe678ddc8b3614aa616c044470 Mon Sep 17 00:00:00 2001
From: Berenger Bramas <bbramas@mpcdf.mpg.de>
Date: Thu, 27 Apr 2017 16:05:32 +0200
Subject: [PATCH] Use floor to find the correct field cell

---
 bfps/cpp/particles/particles_field_computer.hpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/bfps/cpp/particles/particles_field_computer.hpp b/bfps/cpp/particles/particles_field_computer.hpp
index 96d58a26..8d91ac82 100644
--- a/bfps/cpp/particles/particles_field_computer.hpp
+++ b/bfps/cpp/particles/particles_field_computer.hpp
@@ -192,7 +192,7 @@ public:
 
     virtual int pbc_field_layer(const real_number& a_z_pos, const int idx_dim) const final {
         const real_number shifted_pos = a_z_pos - spatial_box_offset[idx_dim];
-        const int nb_level_to_pos = int(shifted_pos/box_step_width[idx_dim]);
+        const int nb_level_to_pos = int(floor(shifted_pos/box_step_width[idx_dim]));
         const int int_field_grid_dim = int(field_grid_dim[idx_dim]);
         const int pbc_level = ((nb_level_to_pos%int_field_grid_dim)+int_field_grid_dim)%int_field_grid_dim;
         assert(0 <= pbc_level && pbc_level < int(field_grid_dim[idx_dim]));
-- 
GitLab