I’m trying to understand why there are separate implementations of the force_clear() method in the Verlet and Respa classes. The method in Verlet seems to be more powerful (using memset if group ‘all’ is being integrated).
Wouldn’t it make sense to have force_clear() be a method of the Integrate base class?
P.S.: I’m asking as I’m working on a proxy pair style that needs to clear forces. And as it currently is I have to duplicate code.