We describe an efficient numerical approach to calculate the longitudinal and transverse Kubo conductivities of large systems using Bastin's formulation [1]. We expand the Green's functions in terms of Chebyshev polynomials and compute the conductivity tensor for any temperature and chemical potential in a single step. To illustrate the power and generality of the approach, we calculate the conductivity tensor for the quantum Hall effect in disordered graphene and analyze the effect of the disorder in a Chern insulator in Haldane's model on a honeycomb lattice.PACS numbers: 71.23. An,72.15.Rn,71.30.+h One of the most important experimental probes in condensed matter physics is the electrical response to an external electrical field. In addition to the longitudinal conductivity, in specific circumstances, a system can present a transverse conductivity under an electrical perturbation. The Hall effect [2] and the anomalous Hall effect in magnetic materials [3] are two examples of this type of response. Paramagnetic materials with spin-orbit interaction can also present transverse spin currents [4]. There are also the quantized versions of the three phenomena: while the quantum Hall effect (QHE) was observed more than 30 years ago [5], the quantum spin Hall effect (QSHE) and the quantum anomalous Hall effect (QAHE) could only be observed [6, 7] with the recent discovery of topological insulators, a new class of quantum matter [8].In the linear response regime, the conductivity tensor can be calculated using the Kubo formalism [9]. The Hall conductivity can be easily obtained in momentum space in terms of the Berry curvature associated with the bands [10]. The downside of working in momentum space, however, is that the robustness of a topological state in the presence of disorder can only be calculated perturbatively [11]. Real-space implementations of the Kubo formalism for the Hall conductivity, on the other hand, allow the incorporation of different types of disorder in varying degrees, while providing flexibility to treat different geometries. Real-space techniques, however, normally require a large computational effort. This has generally restricted their use to either small systems at any temperature [12,13], or large systems at zero temperature [14].In this Letter, we propose a new efficient numerical approach to calculate the conductivity tensor in solids. We use a real space implementation of the Kubo formalism where both diagonal and off-diagonal conductivities are treated in the same footing. We adopt a formulation of the Kubo theory that is known as Bastin formula [1] and expand the Green's functions involved in terms of Chebyshev polynomials using the kernel polynomial method [16]. There are few numerical methods that use Chebyshev expansions to calculate the longitudinal DC conductivity [17][18][19][20] and transverse conductivity [14,21] at zero temperature. An advantage of our approach is the possibility of obtaining both conductivities for large systems in a single calculation step, independe...