This branch should divide / 2 the simulation time step occurring during the Rho update, due to the absence of the dyagonalization of the starting dynamical matrix. Moreover, if the current dynamical matrix is the same as the original one, no diagonalization is performed at all in the step update.
This branch should divide / 2 the simulation time step occurring during the Rho update, due to the absence of the dyagonalization of the starting dynamical matrix. Moreover, if the current dynamical matrix is the same as the original one, no diagonalization is performed at all in the step update.