Autodiff jacobians sometimes don't work

[jonniedie]

using ComponentArrays
using DifferentialEquations
using Parameters: @unpack

tspan = (0.0, 20.0);

function lorenz!(D, u, p, t; f=0.0)
    @unpack σ, ρ, β = p
    @unpack x, y, z = u
    
    D.x = σ*(y - x)
    D.y = x*(ρ - z) - y - f
    D.z = x*y - β*z
    return nothing
end

lorenz_p = (σ=10.0, ρ=28.0, β=8/3)
lorenz_ic = CArray(x=1.0, y=0.0, z=0.0)

lorenz_prob = ODEProblem(lorenz!, lorenz_ic, tspan, lorenz_p)
lorenz_sol = solve(lorenz_prob, Rodas5())

[jonniedie]

It’s failing in another way now for a different problem. It’s really hard to pin down with the debugger due to all the generated functions in Base.Broadcast. I can’t seem to get it to fail outside of a diff eq solving context. I’ll need to add some tests with full implicit solves after this.