Mentions légales du service

Skip to content
Snippets Groups Projects
Commit 0862ae97 authored by BRAMAS Berenger's avatar BRAMAS Berenger
Browse files

update kernel from new leaf

parent 9df34042
No related branches found
No related tags found
No related merge requests found
...@@ -455,12 +455,12 @@ public: ...@@ -455,12 +455,12 @@ public:
FReal legendre[SizeArray]; FReal legendre[SizeArray];
// For all particles in the leaf box // For all particles in the leaf box
const FReal*const physicalValues = inParticles.getPhysicalValues(); const FReal*const physicalValues = inParticles->getPhysicalValues();
const FReal*const positionsX = inParticles.getPositions()[0]; const FReal*const positionsX = inParticles->getPositions()[0];
const FReal*const positionsY = inParticles.getPositions()[1]; const FReal*const positionsY = inParticles->getPositions()[1];
const FReal*const positionsZ = inParticles.getPositions()[2]; const FReal*const positionsZ = inParticles->getPositions()[2];
for(int idxPart = 0 ; idxPart < inParticles.getNbParticles() ; ++ idxPart){ for(int idxPart = 0 ; idxPart < inParticles->getNbParticles() ; ++ idxPart){
// P2M // P2M
const FPoint position(positionsX[idxPart],positionsY[idxPart],positionsZ[idxPart]); const FPoint position(positionsX[idxPart],positionsY[idxPart],positionsZ[idxPart]);
const FSpherical sph(position - cellPosition); const FSpherical sph(position - cellPosition);
...@@ -666,16 +666,16 @@ public: ...@@ -666,16 +666,16 @@ public:
const FPoint cellPosition = getLeafCenter(inLocal->getCoordinate()); const FPoint cellPosition = getLeafCenter(inLocal->getCoordinate());
// For all particles in the leaf box // For all particles in the leaf box
const FReal*const physicalValues = inParticles.getPhysicalValues(); const FReal*const physicalValues = inParticles->getPhysicalValues();
const FReal*const positionsX = inParticles.getPositions()[0]; const FReal*const positionsX = inParticles->getPositions()[0];
const FReal*const positionsY = inParticles.getPositions()[1]; const FReal*const positionsY = inParticles->getPositions()[1];
const FReal*const positionsZ = inParticles.getPositions()[2]; const FReal*const positionsZ = inParticles->getPositions()[2];
const FReal*const forcesX = inParticles.getForcesX(); FReal*const forcesX = inParticles->getForcesX();
const FReal*const forcesY = inParticles.getForcesY(); FReal*const forcesY = inParticles->getForcesY();
const FReal*const forcesZ = inParticles.getForcesZ(); FReal*const forcesZ = inParticles->getForcesZ();
const FReal*const potentials = inParticles.getPotentials(); FReal*const potentials = inParticles->getPotentials();
for(int idxPart = 0 ; idxPart < inParticles.getNbParticles() ; ++ idxPart){ for(int idxPart = 0 ; idxPart < inParticles->getNbParticles() ; ++ idxPart){
// L2P // L2P
const FPoint position(positionsX[idxPart],positionsY[idxPart],positionsZ[idxPart]); const FPoint position(positionsX[idxPart],positionsY[idxPart],positionsZ[idxPart]);
const FSpherical sph(position - cellPosition); const FSpherical sph(position - cellPosition);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment