Hi, I am so confused of the wave function we get from :
wf = kwant.solvers.default.wave_function(some_sys, some_energy) and the eigen vector we can get form theļ¼ ham=sys.hamiltonian_submatrix(args=[Zeeman]) eval,evec=la.eigh(ham) because we can reshape the eigen vector and get its norm, hence give us the wave mode distribution along the wire, which is the real space. So I thought that is the wave function. But when I tried the default.wave_function above, I got something different. I want to make these two consistent with each other. So I wonder what the relation between these two ways. Many Thanks, Chuanchang