renamed _gaussian_error_function to _gaussian_sf which correctly reflects that it is the survival function