Add first optimization of phonon initialization routines.

Refactor few phonon routines of the initial part of the calculation that represent a bottleneck for very large systems

• d2ionq and rotate_patterns ( from !732 (closed) of @bonfus which can be closed now )
• compute_nldyn

more merge requests to follow for the other parts.