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