def conditional_moments(self, f, hyp=None):
        """
        The first two conditional moments of a Gaussian are the mean and variance:
            E[y|f] = f
            Var[y|f] = σ²
        """
        hyp = softplus(self.hyp) if hyp is None else hyp
        return f, hyp