Many thanks,
I have managed Zygote.jl to work, but unfortunately is much slower than ReverseDiff.jl (even for simple real-valued functions). I have open a new question about this, to see if I am missing anything (Zygote terribly slow. What am I doing wrong?).
Thanks!