Sam Isaacson

Yeah, this is an even more general problem I guess. i.e. probably in DiffEqJump user `rate` and `affect` functions should take the `rng` we're using as a parameter in some...

FWIW I was able to get sampling from distributions to work in https://github.com/SciML/Catalyst.jl/pull/493. I also checked and registering a two-argument `rand` that takes a `rng` works fine too. So the...

You mean generate a function from `build_function` that takes the rng as the last parameter and then make a closure?

Is this what gets called? https://github.com/SciML/DiffEqNoiseProcess.jl/blob/3130f51e6acdad1502574dd16b1220b362f8bfe1/src/wiener.jl#L74 I guess a dispatch is needed on `REAL_WHITE_NOISE_DIST`? (Though I can imagine this might require changes here too that make sure `dW` and `W`...

Nice find! I looked for this paper previously without luck...

Thanks for the info! @YingboMa would you consider changing that to preserve the first form?

Or alternatively to give a useful error about double including `tstops`.

Is there a use case for having non-unique initial `tstops` or is that always a user setup error?

This is exactly what I've been thinking about looking at now that most of the non-spatial work is finished, so that sounds great! (And spatial tooling / solvers should also...

I think PDESystems -> MethodOfLines.jl are in good enough shape that we could tackle the reaction-diffusion PDE case today. It is really just making sure everything works when the species...