MomentClosure.jl
MomentClosure.jl copied to clipboard
Support Catalyst named parameter/initial state vectors
Parameters in ModelingToolkit
can be passed as (name, value)
pairs or dictionaries. Example:
rn = @reaction_network begin
s, 0 -> A
d, A -> 0
end
dprob = DiscreteProblem(rn, [ :A => 5 ], (0, 10.), [ :s => 10, :d => 1])
This seems to work via ModelingToolkit.varmap_to_vars
and should probably be added to MomentClosure
Note that with MTK9, and hence Catalyst 14, using named mappings is a requirement as under the hood parameter and species orders can now change.
And a different parameter object is now generated for use in solvers: https://docs.sciml.ai/ModelingToolkit/stable/basics/FAQ/#Why-are-my-parameters-some-obscure-object?
(wrong link initially)