Fix OEP full - spin polarized case

Description

The OEP full was not taking the right occupations for the spin-polarized case.

News snippet

The OEP full was not taking the right occupations for the spin-polarized case.

Checklist

• I have checked that my code follows the Octopus coding standards
• I have added tests for all the new features added in this request.

src/td/propagator_elec.F90

@@ -566,7 +566,6 @@ contains
     if (ion_dynamics_ions_move(ions_dyn) .and. tr%method .ne. PROP_EXPLICIT_RUNGE_KUTTA4) then
       call forces_calculate(gr, namespace, ions, hm, ext_partners, st, ks, t = abs(nt*dt), dt = dt)
       call ion_dynamics_propagate_vel(ions_dyn, ions, atoms_moved = generate)
-      if (generate) call hamiltonian_elec_epot_generate(hm, namespace, space, gr, ions, ext_partners, st, time = abs(nt*dt))
       call ions%update_kinetic_energy()
     else
       if (outp%what(OPTION__OUTPUT__FORCES) .or. write_handler%out(OUT_SEPARATE_FORCES)%write) then