Dear all,

  Imagine that I have an external potential in a form of a list of values, corresponding to the site numbers in finalized system (this has appeared from mean-field account of Hubbard interactions).  Now I want to feed this back to Kwant  to compute the modes.

 The standard function to specify the potential has Site as an argument, with no access to its number in the finalized system.  Is there a fast way to get the site number?

  I found an old reply of Anton:

> site_numbers = {tuple(site.pos(i)): i for i in xrange(sys.graph.num_nodes)}
> site_numbers[(1, 1)]  # Find the number of the site with coordinates (1, 1).
but it does not look efficient.

Any suggestions?

Thanks,
Sergey