J, I = np.meshgrid(ix, ix)