The idea is that we were missing some self interaction because we were setting
cj_nodeID = -1 and then the condition is done in between
cj. I have checked the number of neighbors between
swift_mpi with 1 rank and they agree now.
@matthieu This is my try to fix the bug identified by Alexei. I have only one problem (not sure if it is related). When I try to write a snapshot, it happens that I receive the following error
 [00005.1] cell.c:cell_drift_part():3759: Attempt to drift to the past.
@aborissov can you test this please? On my test, I have the same mass evolution for both serial and parallel.