Matrix elements (integrals) evaluation over Cartesian Gaussian functions
https://github.com/evaleev/libint