Automated selection of optimal Gaussian fits to arbitrary functions in electronic structure theory