Some q points produces wrong results, and the CG algorithm for the hessian does not work.
In this case, a perturbation along a specific q point seems to produce wrong results. For the SnTe, it is the q=1 point (the L point).
This must be fixed before merging the pull request
#18