  The contour integration technique applied to calculate the optical
conductivity tensor at finite temperatures in the case of layered systems
within the framework of the spin-polarized relativistic screened
Korringa-Kohn-Rostoker band structure method is improved from the computational
point of view by applying the Gauss-Konrod quadrature for the integrals along
the different parts of the contour and by designing a cumulative special points
scheme for two-dimensional Brillouin zone integrals corresponding to cubic
systems.
