states_elec_total_density should be rewritten using BLAS operations

This is a minor optimization but this also makes the code easier to read.

Similarly, one could use BLAS axpy for hamiltonian_elec_update_pot too.