From 9d02f75287b62a7597e90409b341870f6e1d0bb2 Mon Sep 17 00:00:00 2001 From: Berenger Bramas Date: Thu, 5 Nov 2015 21:42:58 +0100 Subject: [PATCH] debug periodic --- Src/Core/FFmmAlgorithmPeriodic.hpp | 16 ++++++++-------- 1 file changed, 8 insertions(+), 8 deletions(-) diff --git a/Src/Core/FFmmAlgorithmPeriodic.hpp b/Src/Core/FFmmAlgorithmPeriodic.hpp index 53082658..c2189793 100644 --- a/Src/Core/FFmmAlgorithmPeriodic.hpp +++ b/Src/Core/FFmmAlgorithmPeriodic.hpp @@ -352,19 +352,19 @@ protected: for(int idxNeig = 0 ; idxNeig < counter ; ++idxNeig){ if( !offsets[idxNeig].equals(0,0,0) ){ // Put periodic neighbors into other array - periodicNeighbors[periodicNeighborsCounter] = neighbors[idxNeig]; - periodicNeighborPositions[periodicNeighborsCounter] = neighborPositions[idxNeig]; - ++periodicNeighborsCounter; - - FReal*const positionsX = periodicNeighbors[idxNeig]->getWPositions()[0]; - FReal*const positionsY = periodicNeighbors[idxNeig]->getWPositions()[1]; - FReal*const positionsZ = periodicNeighbors[idxNeig]->getWPositions()[2]; + FReal*const positionsX = neighbors[idxNeig]->getWPositions()[0]; + FReal*const positionsY = neighbors[idxNeig]->getWPositions()[1]; + FReal*const positionsZ = neighbors[idxNeig]->getWPositions()[2]; - for(FSize idxPart = 0; idxPart < periodicNeighbors[idxNeig]->getNbParticles() ; ++idxPart){ + for(FSize idxPart = 0; idxPart < neighbors[idxNeig]->getNbParticles() ; ++idxPart){ positionsX[idxPart] += boxWidth * FReal(offsets[idxNeig].getX()); positionsY[idxPart] += boxWidth * FReal(offsets[idxNeig].getY()); positionsZ[idxPart] += boxWidth * FReal(offsets[idxNeig].getZ()); } + + periodicNeighbors[periodicNeighborsCounter] = neighbors[idxNeig]; + periodicNeighborPositions[periodicNeighborsCounter] = neighborPositions[idxNeig]; + ++periodicNeighborsCounter; } else{ neighbors[idxNeig-periodicNeighborsCounter] = neighbors[idxNeig]; -- GitLab