Skip to content

remake failing for several cases with vector variables #2660

Closed
@TorkelE

Description

@TorkelE

When you have vector variables, it seems like remake is not working for several cases:

using ModelingToolkit, OrdinaryDiffEq

@variables t (X(t))[1:2]
@parameters p[1:2]
D = Differential(t)
eqs = [
    D(X[1]) ~ p[1] - X[1],
    D(X[2]) ~ p[2] - X[2]
]
@mtkbuild osys = ODESystem(eqs, t)

u0 = [X[1] => 1.0, X[2] => 2.0]
tspan = (0.0, 1.0)
ps = [p[1] => 0.1, p[2] => 0.2]
oprob = ODEProblem(osys, u0, tspan, ps)

oprob_rmk = remake(oprob; u0 = [X => [10.0, 20.0]]) # Errors
oprob_rmk = remake(oprob; u0 = [X[1] => 10.0, X[2] => 20.0]) # Errors
oprob_rmk = remake(oprob; u0 = [osys.X => [10.0, 20.0]]) # Errors
oprob_rmk = remake(oprob; u0 = [osys.X[1] => 10.0, osys.X[2] => 20.0]) # Errors
oprob_rmk = remake(oprob; u0 = [:X => [10.0, 20.0]]) # Errors
oprob_rmk = remake(oprob; ps = [p => [10.0, 20.0]]) # Works
oprob_rmk = remake(oprob; ps = [p[1] => 10.0, p[2] => 20.0]) # Works
oprob_rmk = remake(oprob; ps = [osys.p => [10.0, 20.0]]) # Works
oprob_rmk = remake(oprob; ps = [osys.p[1] => 10.0, osys.p[2] => 20.0]) # Works
oprob_rmk = remake(oprob; ps = [:p => [10.0, 20.0]]) # Works

For parameters it seems fine though.

Metadata

Metadata

Labels

bugSomething isn't working

Type

No type

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions