diff --git a/src/exx_erigto.f90 b/src/exx_erigto.f90
index 0386420b..72ffecb9 100644
--- a/src/exx_erigto.f90
+++ b/src/exx_erigto.f90
@@ -666,6 +666,8 @@ subroutine compute_eri_hoh(i,j,k,l,i_nsp,j_nsp,k_nsp,l_nsp,i_xyz,j_xyz,k_xyz,l_x
     
     real(double), dimension(5) :: F
 
+    eri_gto = zero
+
     i_nt = trim(gto( i_nsp )%sf( i )%nt)
     j_nt = trim(gto( j_nsp )%sf( j )%nt)
     k_nt = trim(gto( k_nsp )%sf( k )%nt)