hamiltonian that doesn't produce onsite terms fails to go through discretizer
import discretizer
kx = discretizer.momentum_operators[0]
discretizer.Discretizer(kx)
import discretizer
kx = discretizer.momentum_operators[0]
discretizer.Discretizer(kx)