in pair_dpd.cpp , in the “compute” function , when “j” which a neighbor of particle “i” is such that j>nlocal

how is the how is the force assigned to the correct neighbor, such that the pairwise cancellation of forces correctly

implemented ?

in pair_dpd.cpp , in the "compute" function , when "j" which a

neighbor of particle "i" is such that j>nlocal

how is the how is the force assigned to the correct neighbor, such

that the pairwise cancellation of forces correctly

implemented ?

what pairwise _cancellation_??? please explain.

please have a look at the documentation of newton_pair

and the papers on the lammps algorithm and implementation.

cheers,

axel.

For any pairwise style (not just dpd), the forces are assigned

to ghost atoms and then later communicated. The details

of whether or how this is done is controlled by the newton

command.

Steve