  An efficient method for calculating the electronic structure in large systems
with a fully converged BZ sampling is presented. The method is based on a
k.p-like approximation developed in the framework of the density functional
perturbation theory. The reliability and efficiency of the method are
demostrated in test calculations on Ar and Si supercells.
