I'm trying to add capability for self gravitation for N-Body dust simulations. As I understand it, LIGGGHTS includes coulombic long range solvers, backed up by the "getting started" section of the documentation, however, I've been having a difficult time finding where in the SRC the long range solvers are implemented, such that I can modify them as needed. It should be relatively straight forward, as they are both solving poisson equations, gravity should be as simple as changing a constant, adding a minus sign, and swapping m for q.
Does anybody have any experience with this side of LIGGGHTS?
Thanks,
Thomas