This method is especially suited for getting band edge states in the framework of linear scaling self-consistent density functional or Hartree-Fock theory.