speed optimization of Vacbd (symmetry) + Gauss-Hermite integration (64 points)