Open
Description
this
using Catalyst
rn = @reaction_network begin
(k1,k2), 2X1 <--> X2
(k3,k4), X1 + X2 <--> X3
end
@unpack X1, X2, X3 = rn
u0 = [X1 => 1.0, X2 => 1.0, X3 => 1.0]
ps = [:k1 => 0.1, :k2 => 0.2, :k3 => 0.3, :k4 => 0.4]
prob_old = ODEProblem(rn, u0, 1.0, ps; remove_conserved = true)
prob_new = remake(prob_old; p = [:Γ => [10.0]], warn_initialize_determined = false)
gives a
┌ Warning: Initialization system is overdetermined. 1 equations for 0 unknowns. Initialization will default to using least squares. `SCCNonlinearProblem` can only be used for initialization of fully determined systems and hence will not be used here. To suppress this warning pass warn_initialize_determined = false. To make this warning into an error, pass fully_determined = true
└ @ ModelingToolkit C:\Users\Torkel\.julia\packages\ModelingToolkit\r9J1w\src\systems\diffeqs\abstractodesystem.jl:1446