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/hamiltonian/hamiltonian_elec_base_inc.F90

@@ -860,6 +860,10 @@ subroutine X(hamiltonian_elec_base_nlocal_finish)(this, mesh, bnd, std, projecti
 
       end if
 
+      ! MFT: DEBUG
+      write(*,*)'src/hamiltonian/hamiltonian_elec_base_inc.F90 L861:projection%X(projection)=', projection%X(projection)
+      ! MFT: DEBUG-end
+
       if (npoints /=  0) then
 
         SAFE_ALLOCATE(psi(1:nst_linear, 1:npoints))
@@ -1040,6 +1044,22 @@ contains
         call accel_finish()
 
         call profiling_out(cl_prof)
+        ! MFT: DEBUG
+#ifdef R_TCOMPLEX
+        SAFE_ALLOCATE(projection%X(projection)(1:vpsib%pack_size(1), 1:this%full_projection_size))
+        projection%X(projection) = M_ZERO
+        call accel_read_buffer(projection%buff_projection, &
+          this%full_projection_size*vpsib%pack_size(1), projection%X(projection))
+        write(*,*)'src/hamiltonian/hamiltonian_elec_base_inc.F90 L1048: projection%zprojection=', projection%X(projection)
+#else
+        SAFE_ALLOCATE(projection%X(projection)(1:vpsib%pack_size_real(1), 1:this%full_projection_size))
+        projection%X(projection) = M_ZERO
+        call accel_read_buffer(projection%buff_projection, &
+          this%full_projection_size*vpsib%pack_size_real(1), projection%X(projection))
+        write(*,*)'src/hamiltonian/hamiltonian_elec_base_inc.F90 L1048: projection%dprojection=', projection%X(projection)
+#endif
+        SAFE_DEALLOCATE_A(projection%X(projection))
+        ! MFT: DEBUG-end
       else
 
         buff_proj => projection%buff_projection