Many physical quantities in solid-state physics are calculated from $k$-space summation. For spectral functions, the frequency-dependent factor can be decomposed into energy-conserving $\ensuremath{\delta}$-function part and nondissipative principal value part. A very useful scheme for this summation is tetrahedron method. Tetrahedron methods have been widely used to calculate of part, such as ...