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