Consistent separation of left-, right-propagating and evanescent modes
Right now if a new propagating mode opens, Kwant may count one mode out of the pair as propagating, another one as evanescent due to numerical errors. It should be relatively easy to enforce consistency by using that evanescent modes always come in pairs λ
1/λ^*
.
The relevant code is in the make_proper_modes
function over here.
Edited by Anton Akhmerov