If I have periodic boundary condition in all direction for my simulation. I know particle wise PBC is done in domain.cpp pbc() function. where does the pair wise pbc is incorporate in Lammps. Is the neighbour list calculation done already with that ? Which is that function?.
When atoms are communicated, PBC box lengths are added
as needed (atom_vec_atomic.cpp for example). Thus the neigh
lists and pair computations don't have to worry about PBC directly.