From 4e174c0c5c7ff7debcbc1e2b2783acad5ceef8d7 Mon Sep 17 00:00:00 2001
From: Berenger Bramas <bbramas@mpcdf.mpg.de>
Date: Wed, 26 Oct 2016 11:38:52 +0200
Subject: [PATCH] forgot to remove a bug ffrom mpi periodic

---
 Src/Core/FFmmAlgorithmThreadProcPeriodic.hpp | 21 ++++++++++----------
 1 file changed, 10 insertions(+), 11 deletions(-)

diff --git a/Src/Core/FFmmAlgorithmThreadProcPeriodic.hpp b/Src/Core/FFmmAlgorithmThreadProcPeriodic.hpp
index 80571ba1e..f280c213e 100644
--- a/Src/Core/FFmmAlgorithmThreadProcPeriodic.hpp
+++ b/Src/Core/FFmmAlgorithmThreadProcPeriodic.hpp
@@ -1701,17 +1701,16 @@ protected:
                                             myThreadkernels->P2PRemote(currentIter.coord,currentIter.targets,
                                                                        currentIter.sources, periodicNeighbors, periodicNeighborPositions, periodicNeighborsCounter);
 
-                                            for(int idxNeig = 0 ; idxNeig < counter ; ++idxNeig){
-                                                if( periodicNeighbors[idxNeig] ){
-                                                    FReal*const positionsX = periodicNeighbors[idxNeig]->getWPositions()[0];
-                                                    FReal*const positionsY = periodicNeighbors[idxNeig]->getWPositions()[1];
-                                                    FReal*const positionsZ = periodicNeighbors[idxNeig]->getWPositions()[2];
-
-                                                    for(FSize idxPart = 0; idxPart < periodicNeighbors[idxNeig]->getNbParticles() ; ++idxPart){
-                                                        positionsX[idxPart] -= boxWidth * FReal(offsets[idxNeig].getX());
-                                                        positionsY[idxPart] -= boxWidth * FReal(offsets[idxNeig].getY());
-                                                        positionsZ[idxPart] -= boxWidth * FReal(offsets[idxNeig].getZ());
-                                                    }
+                                            for(int idxNeig = 0 ; idxNeig < periodicNeighborsCounter ; ++idxNeig){
+                                                FAssertLF( periodicNeighbors[idxNeig] );
+                                                FReal*const positionsX = periodicNeighbors[idxNeig]->getWPositions()[0];
+                                                FReal*const positionsY = periodicNeighbors[idxNeig]->getWPositions()[1];
+                                                FReal*const positionsZ = periodicNeighbors[idxNeig]->getWPositions()[2];
+
+                                                for(FSize idxPart = 0; idxPart < periodicNeighbors[idxNeig]->getNbParticles() ; ++idxPart){
+                                                    positionsX[idxPart] -= boxWidth * FReal(offsets[idxNeig].getX());
+                                                    positionsY[idxPart] -= boxWidth * FReal(offsets[idxNeig].getY());
+                                                    positionsZ[idxPart] -= boxWidth * FReal(offsets[idxNeig].getZ());
                                                 }
                                             }
                                         }
-- 
GitLab