Likely, we have not wrapped some particle, as running with 1 MPI process is not going through comm_plasma_y.
(gdb) p x[1][iv]
$2 = -0.00022396883585063811
The boundary conditions are applied in comm_plasma_x, which is a bit weird, but allows us to only use 2 if's per particle to determine if they must move outside the chunk.
We could move the boundary condition outside the comm_plasma, but then we couldn't check again we have no lost particles left.
Likely, we have not wrapped some particle, as running with 1 MPI process is not going through comm_plasma_y.
The boundary conditions are applied in comm_plasma_x, which is a bit weird, but allows us to only use 2 if's per particle to determine if they must move outside the chunk.
We could move the boundary condition outside the comm_plasma, but then we couldn't check again we have no lost particles left.