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.

share/opencl/projector.cl

@@ -436,6 +436,7 @@ __kernel void projector_ket_phase(const int nmat,
   const int scal_offset   = offsets[OFFSET_SIZE*imat + 4];
 
   if(ip < ip_start || ip >= ip_end) return;
+  // if(ip >= npoints) return;
 
   double2 aa = 0.0;
   for(int ipj = 0; ipj < nprojs; ipj++){
@@ -471,6 +472,7 @@ __kernel void projector_ket_phase_spiral(const int nmat,
   const int scal_offset   = offsets[OFFSET_SIZE*imat + 4];
 
   if(ip < ip_start || ip >= ip_end) return;
+  // if(ip >= npoints) return;
 
   double2 aa = 0.0;
   for(int ipj = 0; ipj < nprojs; ipj++){
@@ -493,9 +495,9 @@ __kernel void projector_commutator_ket(const int nmat,
                __global double const * restrict projection, const int ldprojection,
                __global double * restrict cpsi1,
                __global double * restrict cpsi2,
-               __global double * restrict cpsi3, const int ldpsi,
-               const int ip_start, 
-               const int ip_end
+               __global double * restrict cpsi3, const int ldpsi
+               // const int ip_start, 
+               // const int ip_end
                ){
   
   const int ist = get_global_id(0);
@@ -555,6 +557,7 @@ __kernel void projector_commutator_ket_phase(const int nmat,
   const int scal_offset   = offsets[OFFSET_SIZE*imat + 4];
 
   if(ip < ip_start || ip >= ip_end) return;
+  // if(ip >= npoints) return;
 
   double2 aa0 = 0.0;
   double2 aa1 = 0.0;
@@ -608,6 +611,7 @@ __kernel void projector_commutator_ket_phase_spiral(const int nmat,
   const int scal_offset   = offsets[OFFSET_SIZE*imat + 4];
 
   if(ip < ip_start || ip >= ip_end) return;
+  // if(ip >= npoints) return;
 
   double2 aa0 = 0.0;
   double2 aa1 = 0.0;