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