Commensurate matrices for the phonon calculations

We have to dynamically determine the most compact matrix (i.e. smallest determinant) to commensurate with the maximum special BZ points in order to reproduce vibrational thermodymamics as precise as possible with the minimum costs.