Dear Pablo,

in the example here:

https://kwant-project.org/doc/dev/tutorial/kpm

you write:

..

fsyst_staggered = make_syst_staggered().finalized() # find 'A' and 'B' sites in the unit cell at the center of the disk center_tag = np.array([0, 0]) where = lambda s: s.tag == center_tag # make local vectors vector_factory = kwant.kpm.LocalVectors(fsyst_staggered, where)

and plot their respective local density of states.

# 'num_vectors' can be unspecified when using 'LocalVectors' local_dos = kwant.kpm.SpectralDensity(fsyst_staggered, num_vectors=None, vector_factory=vector_factory, mean=False, rng=0) energies, densities = local_dos()

.

My question is:  how do I know which sublattice LDOS comes out first in the result?

I plan to compute LDOS in a lot of point (say, 50 points) at a given set of energies.  So, what will be the ordering of vectors given by vector_factory.   Or, can I create those vectors manually without the vector factory?

Thank you,

Sergey