i, j = np.mgrid[0 : len(radii_a), 0 : len(radii_b)]