DiffEqBase.jl
DiffEqBase.jl copied to clipboard
Incorrect default Jacobian for DDE problems
The default implementation for DDE problems in lines https://github.com/JuliaDiffEq/DiffEqBase.jl/blob/master/src/diffeqfunction.jl#L536-L542 (and similar below) does not make sense since for DDE problems the Jacobian does also depend on h. However, I guess the DiffEqOperator interface currently assumes that the arguments of update_coefficients! and update_coefficients are (J, u, p, t) and (u, p, t), respectively, so I'm not sure how to fix this issue.