Sam Isaacson

Also note, `rx.rate` is not the full state-dependent transition rate for a multiparticle reaction like `A --> B`, which should be `rx.rate*A(t)`. This is given by the associated `jumpratelaw`.

@rveltz once you've got a symbolic representation for the generator you should be able to create a ModelingToolkit `ODESystem` with the corresponding ODEs, from which you can instantiate it by...

In this case it seems like what we want is a differential-difference equation system type to encode the infinite system (or can the PDE system type handle that?).

Got it. I'll try to put an example together for you.

Is this what you are looking for? ```julia using Catalyst, ModelingToolkit, Symbolics, LinearAlgebra rn = @reaction_network begin k1, S1 --> S2 k2, S2 --> S1 k3, S3 --> S1 end...

Some modalities to consider: - Reactions on graphs or meshes; given an implicit or explicit connectivity (i.e. a graph / adj. matrix) allow users to specify several types of reactions...

I noticed the `maxiters` warning is also showing up in DiffEqBiological tests now, from "func_test.jl", when `solve` is called on each of the three `SDEProblem`s. I'm not sure what it...

On how to handle solutions going negative; my understanding is that the underlying CLE models are not really well-defined once a solution component becomes zero. So I'm not sure one...

Another fix, which would also be a fair amount of work, would be to provide a hybrid algorithm that dynamically partitions reactions involving species with small populations to be represented...

Sure, I'm in complete agreement that is something we should build towards. I'd like to be able to use that functionality myself (actually, to use it in spatial systems...). At...