diff --git a/cpp/particles/interpolation/particles_field_computer.hpp b/cpp/particles/interpolation/particles_field_computer.hpp
index 1a7b3d1816189aeefa018d378e5362bf9f3d1240..5cdd9eee01dd2ac8c53e95af57c3920f515cbe0b 100644
--- a/cpp/particles/interpolation/particles_field_computer.hpp
+++ b/cpp/particles/interpolation/particles_field_computer.hpp
@@ -93,6 +93,7 @@ public:
     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);
diff --git a/cpp/particles/p2p/p2p_computer.hpp b/cpp/particles/p2p/p2p_computer.hpp
index 49115bd2f5cc626fbae7758e45cf2ec6ef9d8344..52e34734fa6fd2ecd89c896be9c04881e71861d8 100644
--- a/cpp/particles/p2p/p2p_computer.hpp
+++ b/cpp/particles/p2p/p2p_computer.hpp
@@ -55,6 +55,7 @@ public:
 
     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);
     }