Efficient computation of quantum mechanical matrix elements over Gaussian basis sets
