diff --git a/cpp/particles/interpolation/particles_field_computer.hpp b/cpp/particles/interpolation/particles_field_computer.hpp
index 5cdd9eee01dd2ac8c53e95af57c3920f515cbe0b..172ca37a0b3a303d9a25a966b682920f40c5822b 100644
--- a/cpp/particles/interpolation/particles_field_computer.hpp
+++ b/cpp/particles/interpolation/particles_field_computer.hpp
@@ -89,16 +89,6 @@ public:
     /// Computation related
     ////////////////////////////////////////////////////////////////////////
 
-    template <int size_particle_rhs>
-    void init_result_array(real_number particles_current_rhs[],
-                                   const partsize_t nb_particles) const {
-        // Set values to zero initialy
-        // TODO: openmp
-        std::fill(particles_current_rhs,
-                  particles_current_rhs+nb_particles*size_particle_rhs,
-                  0);
-    }
-
     const std::array<real_number, 3> getSpatialBoxWidth() const
     {
         return this->spatial_box_width;
diff --git a/cpp/particles/p2p/p2p_computer.hpp b/cpp/particles/p2p/p2p_computer.hpp
index 52e34734fa6fd2ecd89c896be9c04881e71861d8..5a399e515da2bf59bdfcb3874c7591daf0af1b4c 100644
--- a/cpp/particles/p2p/p2p_computer.hpp
+++ b/cpp/particles/p2p/p2p_computer.hpp
@@ -53,12 +53,6 @@ class p2p_computer{
 public:
     p2p_computer() : isActive(true){}
 
-    template <int size_particle_rhs>
-    void init_result_array(real_number rhs[], const partsize_t nbParticles) const{
-        // TODO: openmp
-        memset(rhs, 0, sizeof(real_number)*nbParticles*size_particle_rhs);
-    }
-
     template <int size_particle_rhs>
     void reduce_particles_rhs(real_number rhs_dst[], const real_number rhs_src[], const partsize_t nbParticles) const{
         static_assert(size_particle_rhs == 6, "This kernel works only with 6 values per particle's rhs");
diff --git a/cpp/particles/p2p/p2p_computer_empty.hpp b/cpp/particles/p2p/p2p_computer_empty.hpp
index 5c43ef3c8a13d3095e3b172787bca8a8bdbe8c3b..b13a95e1f2bb8fd4f51f51c4cb3b6f74e00b079b 100644
--- a/cpp/particles/p2p/p2p_computer_empty.hpp
+++ b/cpp/particles/p2p/p2p_computer_empty.hpp
@@ -31,9 +31,6 @@
 template <class real_number, class partsize_t>
 class p2p_computer_empty{
 public:
-    template <int size_particle_rhs>
-    void init_result_array(real_number /*rhs*/[], const partsize_t /*nbParticles*/) const{
-    }
 
     template <int size_particle_rhs>
     void reduce_particles_rhs(real_number /*rhs_dst*/[], const real_number /*rhs_src*/[], const partsize_t /*nbParticles*/) const{
diff --git a/cpp/particles/p2p/p2p_ghost_collisions.hpp b/cpp/particles/p2p/p2p_ghost_collisions.hpp
index 4341fd39c866561388f704a634cd50749c54fd4b..5dee282ac94c3f8b4a9427477e8a5a14936139b0 100644
--- a/cpp/particles/p2p/p2p_ghost_collisions.hpp
+++ b/cpp/particles/p2p/p2p_ghost_collisions.hpp
@@ -118,10 +118,6 @@ public:
         return min_distance;
     }
 
-    template <int size_particle_rhs>
-    void init_result_array(real_number /*rhs*/[], const partsize_t /*nbParticles*/) const{
-    }
-
     template <int size_particle_rhs>
     void reduce_particles_rhs(real_number /*rhs_dst*/[], const real_number /*rhs_src*/[], const partsize_t /*nbParticles*/) const{
     }
diff --git a/cpp/particles/p2p/p2p_merge_collisions.hpp b/cpp/particles/p2p/p2p_merge_collisions.hpp
index cf93f3e8001699b3cb8de2db4c043409bae4861b..4b728c8f3c3db1c26e41947fe1a23325fbe688e8 100644
--- a/cpp/particles/p2p/p2p_merge_collisions.hpp
+++ b/cpp/particles/p2p/p2p_merge_collisions.hpp
@@ -37,10 +37,6 @@ public:
     // Copy constructor use a counter set to zero
     p2p_merge_collisions(const p2p_merge_collisions&){}
 
-    template <int size_particle_rhs>
-    void init_result_array(real_number /*rhs*/[], const partsize_t /*nbParticles*/) const{
-    }
-
     template <int size_particle_rhs>
     void reduce_particles_rhs(real_number /*rhs_dst*/[], const real_number /*rhs_src*/[], const partsize_t /*nbParticles*/) const{
     }