code: def dKdiag_dtheta(self, dL_dKdiag, X, target):
          
          target += 2.*self.mapping.df_dtheta(dL_dKdiag[:, None]*self.mapping.f(X), X)
docs: Gradient of diagonal of covariance with respect to parameters.