I am creating functionality which makes the Hamiltonian block diagonal. Currently this is implemented by taking the total z-magnetization as a good quantum number.