From bdbb08207956051bb427051f184449985628fd23 Mon Sep 17 00:00:00 2001
From: Emmanuel Giner <giner.emmanuel@gmail.com>
Date: Sun, 8 Nov 2020 17:53:37 +0100
Subject: [PATCH] fixed a little bun in src/utils/one_e_integration.irp.f

---
 src/utils/one_e_integration.irp.f | 14 +++++++-------
 1 file changed, 7 insertions(+), 7 deletions(-)

diff --git a/src/utils/one_e_integration.irp.f b/src/utils/one_e_integration.irp.f
index 97eef89d..2c597104 100644
--- a/src/utils/one_e_integration.irp.f
+++ b/src/utils/one_e_integration.irp.f
@@ -53,13 +53,13 @@ subroutine overlap_gaussian_xyz(A_center,B_center,alpha,beta,power_A,&
   integer                        :: iorder_p(3)
 
   call give_explicit_poly_and_gaussian(P_new,P_center,p,fact_p,iorder_p,alpha,beta,power_A,power_B,A_center,B_center,dim)
-!  if(fact_p.lt.1d-20)then
-!    overlap_x = 0.d0
-!    overlap_y = 0.d0
-!    overlap_z = 0.d0
-!    overlap = 0.d0
-!    return
-!  endif
+   if(fact_p.lt.1d-20)then
+     overlap_x = 1.d-10
+     overlap_y = 1.d-10
+     overlap_z = 1.d-10
+     overlap = 1.d-10
+     return
+   endif
   integer                        :: nmax
   double precision               :: F_integral
   nmax = maxval(iorder_p)