Do not set to zero inside the P2M

Src/Kernels/Chebyshev/FChebInterpolator.hpp

@@ -638,7 +638,7 @@ inline void FChebInterpolator<ORDER,MatrixKernelClass>::applyP2M(const FPoint& c
                                                                  const ContainerClass *const inParticles) const
 {
     // set all multipole expansions to zero
-    FBlas::setzero(2*nRhs*nnodes, multipoleExpansion);
+    // Should not be set to zero here (TODO remove this line) FBlas::setzero(2*nRhs*nnodes, multipoleExpansion);
 
     // allocate stuff
     const map_glob_loc map(center, width);

Src/Kernels/Uniform/FUnifInterpolator.hpp

@@ -416,7 +416,7 @@ inline void FUnifInterpolator<ORDER,MatrixKernelClass>::applyP2M(const FPoint& c
                                                                  const ContainerClass *const inParticles) const
 {
   // set all multipole expansions to zero
-  FBlas::setzero(nRhs*nnodes, multipoleExpansion);
+  // Should not be set to zero here (TODO remove this line) FBlas::setzero(nRhs*nnodes, multipoleExpansion);
 
   // allocate stuff
   const map_glob_loc map(center, width);