Abstract
A new flexible working synthesis and extension of the methods for linear integration over the Brillouin zone, used for obtaining densities of states and generalised susceptibilities is given. The well known Gilat-Raubenheimer and tetrahedron techniques are analysed within this model. Its power and simplicity are demonstrated in a direct calculation of the imaginary part of the generalised susceptibility. The model is based on convex polytope geometry.