Fix and improvement of the Slater potential

Description

Fix some problems with the Slater potential and improve the numerical stability for the case of singular matrices. This implies a local rotation to the frame for which the spin-density matrix is singular, and to solve the Sylvester equation in the local frame, and to set the last, undefined, component of the exchange potential to zero.

News snippet

Fix some problems with the Slater potential and improve the numerical stability for the case of singular matrices.

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/system/v_ks.F90

@@ -396,7 +396,11 @@ contains
 
       if(bitand(ks%xc_family, XC_FAMILY_OEP) /= 0) then
         if (gr%have_fine_mesh) call messages_not_implemented("OEP functionals with UseFineMesh", namespace=namespace)
-        call xc_oep_init(ks%oep, namespace, ks%xc_family, gr, st, geo, mc)
+        if (ks%xc%functional(FUNC_X,1)%id /= XC_OEP_X_SLATER) then 
+          call xc_oep_init(ks%oep, namespace, ks%xc_family, gr, st, geo, mc)
+        else
+          ks%oep%level = XC_OEP_NONE
+        end if
       else
         ks%oep%level = XC_OEP_NONE
       end if