354
350
! done with nspin=4 (grid_nspin)?
355
351
V(1:nsp,1:nspin) = Vscf(1:nsp,ip,1:nspin)
358
353
C Factor two for nondiagonal elements for non-collinear spin
360
355
V(1:nsp,7:8) = V(1:nsp,3:4)
362
357
V(1:nsp,3:4) = 2.0_dp * V(1:nsp,3:4)
366
360
C Loop on first orbital of mesh point