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