Closed ranocha closed 1 month ago
Note that vector version always worked fine:
using OrdinaryDiffEq
f1(v, q, params, t) = -sin.(q)
f2(v, q, params, t) = v
v0 = [1.0]
q0 = [-1.2]
tspan = (0.0, 55.0)
ode = DynamicalODEProblem(f1, f2, v0, q0, tspan)
sol = solve(ode, Tsit5())
sol(0.5, idxs = [2, 1])
it's just arraypartition of scalars. This has a fix coming in today.
Describe the bug 🐞
Expected behavior
The interpolation works.
Minimal Reproducible Example 👇
See above
Error & Stacktrace ⚠️
See above
Environment (please complete the following information):
Additional context
Found in https://github.com/ranocha/BSeries.jl/actions/runs/8965125375/job/24618039980#step:9:44