JumpProcesses.jl icon indicating copy to clipboard operation
JumpProcesses.jl copied to clipboard

Set max number of jumps

Open gaurav-arya opened this issue 1 year ago • 3 comments

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

gaurav-arya avatar Jul 07 '23 14:07 gaurav-arya

Couldn't you track this via a DiscreteCallback and just call terminate! when you hit your limit?

isaacsas avatar Jul 07 '23 15:07 isaacsas

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.

gaurav-arya avatar Jul 07 '23 16:07 gaurav-arya

We should add maxiters keyword support for SSAStepper at least.

isaacsas avatar Feb 29 '24 17:02 isaacsas