Benchmarking of tau leap

[j-c-gibson]

I've noticed that stochastic simulations run without end if all state populations hit zero. I think this is because the adaptive tau leap algorithm is allocating very small or possibly zero step sizes.

It would be good to fix this and also in general to run some checks to verify that the step sizes are sensible.