JumpProcesses.jl
JumpProcesses.jl copied to clipboard
Set max number of jumps
Is there a way of setting a max number of jumps, at which point to abort with an incomplete solution? This would be very useful when dealing with time-dependent jumps with Coevolve
, where I would like my simulation to naturally abort if my jump rates are larger than I expected
Couldn't you track this via a DiscreteCallback
and just call terminate!
when you hit your limit?
It seems like it's getting stuck in between DiscreteCallback
's, so that I can't catch the infinite loop in the DiscreteCallback
. It'll take me a bit to make an MWE, but from my logs it seems like Coevolve
is getting stuck in a loop at a point where all reaction rates are 0.
We should add maxiters
keyword support for SSAStepper
at least.