I have looked at lammps developer.pdf and followed the guidence of wirting a new Fix which stores atoms position from previous time step but it didn't work! I think the code is not complete. For example, the int setmask() is a pure virtual function in the fix.h, it must be defined in child classes. But is it enough? How to complete this fix?