code: def log_likelihood(self):
          
          ll = GP.log_likelihood(self)
          jacobian = self.warping_function.fgrad_y(self.Y_untransformed)
          return ll + np.log(jacobian).sum()
docs: Notice we add the jacobian of the warping function here.