def variational_expectation(self, y, m, v, hyp=None):
        """
        If no custom variational expectation method is provided, we use Gauss-Hermite quadrature.
        """
        return self.variational_expectation_quadrature(y, m, v, hyp)