smooth_erf_locpol function

Compute smoothed erf with locpol approach