Hugely speed up smearing using np.apply_along_axis; remove verbose option