i, j = np.mgrid[-radius : radius + 1, -radius : radius + 1].astype(float) / radius