Bug in integrals

2024-12-23 04:43:50 +01:00 · 2016-10-10 17:08:53 +02:00 · 2016-10-10 17:08:53 +02:00 · 6e897584e9
commit 6e897584e9
parent 764a46230f
1 changed files with 3 additions and 3 deletions
--- a/src/Integrals_Bielec/mo_bi_integrals.irp.f
+++ b/src/Integrals_Bielec/mo_bi_integrals.irp.f
@ -240,9 +240,9 @@ subroutine add_integrals_to_map(mask_ijkl)
              exit
            endif
            bielec_tmp_1(i) = c*bielec_tmp_3(i,j0,k0)
-!          enddo
-!          
-!          do i = 1, min(k,j1-i1)
+          enddo
+          
+          do i = 1, min(k,j1-i1)
            if (abs(bielec_tmp_1(i)) < mo_integrals_threshold) then
              cycle
            endif