plt.imshow(np.log(np.abs(J)))