HumpyBlumpy

My solution currently is to redefine the Hamiltonian every time step, i.e calling `non_uniform_ham = repeat(ham,length(ψ));`. Is it correct that I should be recalculating the environments like so `ψ_envs =...