How to compute partial derivatives efficiently with Zygote

Thank you for the pointers, and let’s hope so!