Torkel Loman

Just also mentioning that this would be useful in Catalyst.jl as well. Every once in a while we get an issue where someone has tried to run something in the...

Seems that (yet again) when I start to believe that I'm getting a vague grip on how to solve things numerically, something will pop up to rudely prove me the...

On a general level, it is kinda cool that there's this kind of strange behaviours in numerical systems, makes the field seem quite cool. On a local level, it kind...

Actually, I think I have found similar errors, but for other parameter sets, for the Stratonovich methods as well.

For discrete spatial models, would boundary conditions need to be defined? I guess one could e.g. define a series of containers with diffusion, with the two ends connected to the...

Yes, that's it. Since what the `@add_reactions` macro essentially does is creating a new `ReactionSystem` of your network, and then merge it with the previous one, every set of input...

I vaguely think this was the case in the beginning, but it is not that any longer. The disadvantage of the current system is that it is less flexible. Also,...

The second one seems good for now. If this is something that will actually happen occasionally, we probably should implement a less hacky solution. ```julia @reaction_network begin k1*X1, ∅ →...

I'd say that if ```julia rn = @reaction_network begin r1s... end p1... @add_reactions rn begin r2s... end p2... ``` is a valid notation, then ```julia rn = @reaction_network begin r1s......

Would it make sense to allow an arbitrary number of graphs to be in the input, and then reactions could use "neighbours" through any select graph? In the case of...