How is periodic boundary condition is obtained in pair_style dpd ?
I mean where it is calculating periodic boundary in pair_dpd.cpp (in compute function?) ?
or where it checks when the distance between two particle is calculated?
PBC are enforced for all potentials in domain::pbc().
Within dpd, the distances are calculated between owned
and ghost particles as necessary, so pbc enforcement is