vf = np.vectorize(f, **kwargs)