Cannot interpolate currents if there are coinciding sites
I recently faced this problem when trying to plot currents of systems where two sites have coinciding positions. As an application example, if one creates a tight-binding model for bilayer graphene that is purely two-dimensional, two sites at different layers have coinciding positions in 2D.
The problem occurs because interpolate_current
generates an array with several nan
. The origin of this problem are the 0/0
operations when computing the current direction.