do not allocate/free M dynamically in demag convolution
send various messages to DEBUG instead of LOG
reduce number of integration points in phonon heat capacity estimation
add PI^4 constant
improve calculation of phonon's heat capacity
improve phonon's heat capacity unit test