Fast Rcpp implementation of Gauss–Hermite quadrature
