47 int locIC = this->_sLattice.getLoadBalancer().loc(globic);
50 T physPosP[] { p->getPos()[0], p->getPos()[1], p->getPos()[2] };
53 int latticeRoundedPosP[] {0, 0, 0};
54 this->_sLattice.getCuboidGeometry().get(globic).getLatticeR( latticeRoundedPosP, physPosP );
56 auto eul2LagrRho = _sLattice.getBlock(locIC).get (
57 latticeRoundedPosP[0], latticeRoundedPosP[1], latticeRoundedPosP[2] ).template getFieldPointer<descriptors::EUL2LAGR>();
58 eul2LagrRho[0] += 1. / _nParticles;