From d7bc6088205ee678b34cbc946b79b6f6f62b5964 Mon Sep 17 00:00:00 2001 From: Kevin Gasperich Date: Thu, 30 Jan 2020 17:00:44 -0600 Subject: [PATCH] minor change to complex integral maps --- src/ao_two_e_ints/map_integrals.irp.f | 19 +++++++++++-------- 1 file changed, 11 insertions(+), 8 deletions(-) diff --git a/src/ao_two_e_ints/map_integrals.irp.f b/src/ao_two_e_ints/map_integrals.irp.f index 5ad63c6a..ac32382d 100644 --- a/src/ao_two_e_ints/map_integrals.irp.f +++ b/src/ao_two_e_ints/map_integrals.irp.f @@ -364,18 +364,17 @@ subroutine ao_two_e_integral_periodic_map_idx_sign(i,j,k,l,use_map1,idx,sign) ! ! ! for , conditionals are [a | | | | | | | | + ! | NEW | | | | | | | | | ! +---------+---------+---------+---------+---------+---------+---------+---------+---------+ ! | | m1 | m1* | m2 | m2* | ! +---------+---------+---------+---------+---------+---------+---------+---------+---------+ ! | | TTT | TTF | FFT | FFF | FTT | TFF | TFT | FTF | ! | | 0TT | T0F | 0FT | F0F | | | | | ! | | T0T | 0TF | F0T | 0FF | | | | | - ! | | | | | | TT0 | | FF0 | | + ! | | TT0 | | FF0 | | FT0(r) | TF0(r) | | | ! +---------+---------+---------+---------+---------+---------+---------+---------+---------+ - ! | | FT0 | TF0 | | | | | | | ! | | 00T | 00F | | | | | | | ! | | 000 | | | | | | | | ! +---------+---------+---------+---------+---------+---------+---------+---------+---------+ @@ -396,17 +395,21 @@ subroutine ao_two_e_integral_periodic_map_idx_sign(i,j,k,l,use_map1,idx,sign) idx = 2*idx-1 - if (ij==kl) then !real, map1 + if (ij==kl) then !real, J -> map1, K -> map2 sign=0.d0 - use_map1=.True. + if (i==k) then + use_map1=.True. + else + use_map1=.False. + endif else if (ik.eq.jl) then if (i.lt.k) then !TT0 sign=1.d0 - use_map1=.False. + use_map1=.True. else !FF0 sign=-1.d0 - use_map1=.False. + use_map1=.True. endif else if (i.eq.k) then if (j.lt.l) then !0T*