Motivation. Among the molecular integrals, the two-center overlap integrals play a major role in any accurate molecular structure calculation. They are central to the calculation of multicenter overlap integrals when using the series expansion formulae for Slater type functions about a new center. Consequently, these integrals require an accurate and fast numerical evaluation. Recently, we show...