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