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