On simulating strongly-interacting, stochastic population models

Simulating strongly-interacting biological populations over many orders-of-magnitude can be challenging even on fast modern computers. Here we describe a regime-switching technique that can speed up simulations from stochastic models by large factors, while retaining third-order accuracy at each time-step. We apply it to an elementary model of the immune response to an infectious organism.