"impose model constraints" should only be called *after* all Runge Kutta
substeps have been computed.
Otherwise it messes up the numerical method.

Since rhs objects can contain collision counters, it means we need to be
able to update them, they cannot be const.

problem found by gcc 10