Hello, I am plotting the wavefunction of my system as such: import scipy.linalg as la
ham = sys.hamiltonian_submatrix(params=dict(pot=0)) evecs = la.eigh(ham)[1] wf = abs(evecs[:, 1])**2 kwant.plotter.density(sys, wf) I am wondering how/if I could give plotter.density some condition to only plot part of my system? Some region's wavefunction values don't appear as they are too small, but those are the ones I am trying to observe.
