Sam Isaacson

We could add something like this, but it would require propagating it throughout JumpProcesses as a new jump type, and updating ModelingToolkit JumpSystems too. Also, in your example I assume...

As an alternative, perhaps structural simplification could be moved into `symbolic_discretize`, since that is still a symbolic component, and then users could be pointed to using that function that to...

@bjarthur what you are asking for was recently added by @YingboMa. I don't think it is documented yet, but you can see an example in https://github.com/SciML/ModelingToolkit.jl/blob/2e53e8947e5761daad95be4aef4f9d1591ee1bb9/test/sdesystem.jl#L588-L612 of constructing an `SDESystem`...

@ChrisRackauckas maybe it would make sense to allow passing non-regular jump types to regular jump methods too so they can specialize their code? This would allow @c-louhichi to specialize on...

@c-louhichi can you also update https://github.com/SciML/JumpProcesses.jl/blob/master/test/regular_jumps.jl with an example using / testing the output of your new solver?

We probably don't want to have an expression-based analysis here, but let's get a test in and then I can leave more detailed feedback.

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

We should add `maxiters` keyword support for `SSAStepper` at least.

I don't think this can be done currently unfortunately, but @Vilin97 can hopefully better comment on this. Allowing such reactions is definitely on our TODO list, but I'm not sure...

How is the information in the change in the rate due to the change in the number of `A` going to be propagated across all spatial locations?