Sparse matrices for both J and Jperp. Now generates 2D Ising models on periodic torii (which are easier computationally).