Zygote.hessian of loss function when defining ODEs with modeling toolkit

[cefitzg]

Forwarddiff works:

hessian = ForwardDiff.hessian(z -> model_loss(z), optimized_parameters)