Fast Rcpp implementation of Gauss-Hermite quadrature