title: Partial derivatives of Gaussian Process wrt features
body: <p>Given a Gaussian Process Model with multidimensional features and scalar observations, how do I compute derivatives of the output wrt to each input, in GPyTorch or GPFlow (or scikit-learn)?</p>