Sam Isaacson

That was an internal hack to allow unpacking all the Symbolics in a model, and was never exported or advertised/made api. It really isn’t needed anymore now that `Symbol`-based indexing...

This doesn't work, but at least it errors: ```julia rn = @reaction_network begin @parameters M @species A(t) = M (1.0,1.0), A B end osys = complete(convert(ODESystem, rn; remove_conserved = true))...

~~Continuing with respect to `Γ[1]`, the conservation constant, seems to (silently) give wrong answers:~~ ```julia using Catalyst, BifurcationKit, Plots rn = @reaction_network begin (1.0,1.0), A B end osys = structural_simplify(convert(ODESystem,...

Seems like we need to disallow changing initial conditions and/or the conservation constant in the current extension setup.

@rveltz is there any way you could use `remake` for updating parameters and initial conditions in the bifurcation solvers (at least when using problems coming via MTK)? I guess right...

Does that hook into the symbolic `SciMLBase.remake`?

Yeah, so I suspect that if one has parameters that are functions of initial conditions, initial conditions that are functions of parameters, or parameters that are functions of other parameters,...

Maybe this would be easier by having a `reeval` type function that a user can pass in (in this case Catalyst/MTK)? That could handle updating all other parameters when one...

I guess trying to worry about changing a parameter in an initial condition is a bit much, but I'll try to work up a simple example where there are dependent...

Yeah, this is the kind of issue I was worried about. Thanks for the simple example.