This method has been prototyped by Adrien Sorgniard with contributions
by Michal Nowak.  The prototype has been cleaned-up, debugged, optimized,
and integrated into Kwant by Christoph Groth.