code: def _get_warped_mean(self, mean, std, pred_init=None, deg_gauss_hermite=20):
          
          gh_samples, gh_weights = np.polynomial.hermite.hermgauss(deg_gauss_hermite)
          gh_samples = gh_samples[:, None]
          gh_weights = gh_weights[None, :]
          return gh_weights.dot(self._get_warped_term(mean, std, gh_samples)) / np.sqrt(np.pi)
docs: Calculate the warped mean by using Gauss-Hermite quadrature.