Hello, We're dealing with a 5000 x 5000 site system, and the time to interpolate the densities and currents using kwant.plotter.interpolate_density and kwant.plotter.interpolate_current are by far the rate limiting steps in the calculation:
System assembly : 15 mins Computing scattering states: 3 hrs Interpolating density and current (for use in plotting): 6 hrs! Is there any way to speed up the interpolation routines? I see that the present implementation in plotter.py involves a for loop (in python) over each site/hopping. Perhaps a vectorized implementation may help? Thanks, Mani