Sam Isaacson

Maybe something like this: ```julia rn1 = @reaction_network rn_observed begin @variables X1(t) @observables X2 ~ X1 k, A --> 0 end t = Catalyst.default_t() @variables X2(t) # this is false...

@TorkelE @vyudu this is fixed now right?

Yes, it would be nice if all of these are supported (assuming they work in the symbolic interface).

This needs to be updated to master.

@vyudu when you have a chance can you update to master? (No rush though as I’m traveling this week as we discussed.)

bump

(This is also consistent with the performance difference between `MassActionJump`s and a hand-coded `Direct` method that evaluates propensities via a single function.)

Yeah, I'm not yet convinced of this result, but the conversion would be easy to add and allow us to test this systematically from Catalyst-generated models in SciMLBenchmarks.

I think we really want to avoid defining new versions of existing types like systems and problems. They have enormous amounts of functionality and dispatches defined on them. Having (separate)...

We can try to have an accessor to indicate / tag that a property is dependent on parameters in initial conditions.