MWE:
using Enzyme
function testODE(du,u,p)
du[1:2] = -p[1]*ones(2,2)*u[1:2]
return nothing
end
tmp1 = [1.0, 1.0]
tmp2 = copy(tmp1)
tmp3 = copy(tmp1)
tmp4 = copy(tmp1)
p = [1.0]
Enzyme.autodiff(testODE,Enzyme.Duplicated(tmp1, tmp2),
Enzyme.Duplicated(tmp3, tmp4),
p)