I wonder how to generate tetrahedra meshes, necessary for Brillouin zone integration. I prefer very simple one based on Lehmann and Taut (1972). Is there any algorithm to generate meshes by giving crystal symmetry, lattice constants, and number of desired tetrahedra etc? I think the algorithm for integration will depend on the quantity that one wants to evaluate. So, only mesh generation part would be sufficient for me.