Skip to content
Snippets Groups Projects

Fix OEP full - spin polarized case

Merged Fix OEP full - spin polarized case
fix_oep_polarized into hotfix-12.2
Merged Nicolas Tancogne-Dejean requested to merge fix_oep_polarized into hotfix-12.2
Viewing commit 471a7e6c
Show latest version
5 files
+ 17
− 17
Compare changes
  • Side-by-side
  • Inline
Files
5
src/hamiltonian/xc_vxc_inc.F90
+ 3
− 3
@@ -741,12 +741,12 @@ contains
dpol = sqrt((d(1) - d(2))**2 + &
M_FOUR*(quantities%rho(ip, 3)**2 + quantities%rho(ip, 4)**2))
vpol = (dedd(ip, 1) - dedd(ip, 2))*(d(1) - d(2))/(dpol + tiny)
vpol = (dedd(ip, 1) - dedd(ip, 2))*(d(1) - d(2))/(SAFE_TOL(dpol, tiny))
vxc(ip, 1) = vxc(ip, 1) + M_HALF*(dedd(ip, 1) + dedd(ip, 2) + vpol)
vxc(ip, 2) = vxc(ip, 2) + M_HALF*(dedd(ip, 1) + dedd(ip, 2) - vpol)
vxc(ip, 3) = vxc(ip, 3) + (dedd(ip, 1) - dedd(ip, 2))*quantities%rho(ip, 3)/(dpol + tiny)
vxc(ip, 4) = vxc(ip, 4) + (dedd(ip, 1) - dedd(ip, 2))*quantities%rho(ip, 4)/(dpol + tiny)
vxc(ip, 3) = vxc(ip, 3) + (dedd(ip, 1) - dedd(ip, 2))*quantities%rho(ip, 3)/(SAFE_TOL(dpol, tiny))
vxc(ip, 4) = vxc(ip, 4) + (dedd(ip, 1) - dedd(ip, 2))*quantities%rho(ip, 4)/(SAFE_TOL(dpol, tiny))
end do
elseif (ispin == SPIN_POLARIZED) then
call lalg_axpy(der%mesh%np, M_ONE, dedd(:, 1), vxc(:, 1))