From 0850fa6f723755a0d5a725cd9a476940e14e2e3c Mon Sep 17 00:00:00 2001 From: Emmanuel Giner Date: Mon, 8 Jun 2020 15:17:53 +0200 Subject: [PATCH] fixed bug in g0_UEG_mu_inf, src/dft_utils_func/on_top_from_ueg.irp.f --- src/dft_utils_func/on_top_from_ueg.irp.f | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/dft_utils_func/on_top_from_ueg.irp.f b/src/dft_utils_func/on_top_from_ueg.irp.f index 24e40353..717081a7 100644 --- a/src/dft_utils_func/on_top_from_ueg.irp.f +++ b/src/dft_utils_func/on_top_from_ueg.irp.f @@ -33,7 +33,7 @@ double precision function g0_UEG_mu_inf(rho_a,rho_b) D = -0.01277d0 E = 0.001859d0 x = -d2*rs - if (dabs(rho) > 1.d-20.and.dabs(x).lt.20.d0) then + if (dabs(rho) > 1.d-20) then rs = (3d0 / (4d0*pi*rho))**(1d0/3d0) ! JT: serious bug fixed 20/03/19 x = -d2*rs if(dabs(x).lt.50.d0)then