def METHOD_NAME(self, Xc, Xe):
    Xall = self.trans(Xc, Xe)
    py, ps2 = self.gp.METHOD_NAME(Xall)
    mu = self.yscaler.inverse_transform(py.reshape((-1, 1)))
    var = self.yscaler.scale_**2 * ps2.reshape((-1, 1))
    return mu, np.clip(var, a_min=np.finfo(var.dtype).eps, a_max=None)