preliminary. in particular solvers cannot be a simple function, they
should handle checkpointing --- especially if Adams-Bashforth particles
are to be reimplemented with the new classes.