2019-01-25 11:39:31 +01:00
|
|
|
BEGIN_PROVIDER [double precision, big_array_coulomb_integrals, (mo_num,mo_num, mo_num)]
|
|
|
|
|
&BEGIN_PROVIDER [double precision, big_array_exchange_integrals,(mo_num,mo_num, mo_num)]
|
|
|
|
|
implicit none
|
|
|
|
|
BEGIN_DOC
|
|
|
|
|
! big_array_coulomb_integrals(i,j) = <ij|ij> = (ii|jj)
|
|
|
|
|
!
|
|
|
|
|
! big_array_exchange_integrals(i,j) = <ij|ji> = (ij|ij)
|
|
|
|
|
END_DOC
|
|
|
|
|
integer :: i,j,k,l
|
|
|
|
|
double precision :: get_two_e_integral
|
|
|
|
|
double precision :: integral
|
|
|
|
|
|
|
|
|
|
do k = 1, mo_num
|
|
|
|
|
do i = 1, mo_num
|
|
|
|
|
do j = 1, mo_num
|
|
|
|
|
l = j
|
|
|
|
|
integral = get_two_e_integral(i,j,k,l,mo_integrals_map)
|
|
|
|
|
big_array_coulomb_integrals(j,i,k) = integral
|
|
|
|
|
l = j
|
|
|
|
|
integral = get_two_e_integral(i,j,l,k,mo_integrals_map)
|
|
|
|
|
big_array_exchange_integrals(j,i,k) = integral
|
|
|
|
|
enddo
|
|
|
|
|
enddo
|
|
|
|
|
enddo
|
|
|
|
|
|
|
|
|
|
END_PROVIDER
|
|
|
|
|
|
2020-02-12 01:23:34 +01:00
|
|
|
BEGIN_PROVIDER [complex*16, big_array_coulomb_integrals_complex, (mo_num,mo_num, mo_num)]
|
|
|
|
|
&BEGIN_PROVIDER [complex*16, big_array_exchange_integrals_complex,(mo_num,mo_num, mo_num)]
|
2020-02-05 21:21:28 +01:00
|
|
|
implicit none
|
|
|
|
|
BEGIN_DOC
|
|
|
|
|
! big_array_coulomb_integrals(j,i,k) = <ij|kj> = (ik|jj)
|
|
|
|
|
! big_array_exchange_integrals(j,i,k) = <ij|jk> = (ij|jk)
|
|
|
|
|
! for both of these, i and k must be from same kpt for integral to be nonzero
|
|
|
|
|
! TODO: only loop over half, and assign two elements:
|
|
|
|
|
! b_a_coul_int(j,i,k) = b_a_coul_int(j,k,i)*
|
|
|
|
|
! b_a_exch_int(j,i,k) = b_a_exch_int(j,k,i)*
|
|
|
|
|
END_DOC
|
|
|
|
|
integer :: i,j,k,l
|
2020-02-12 01:23:34 +01:00
|
|
|
complex*16 :: get_two_e_integral_complex
|
2020-02-05 21:21:28 +01:00
|
|
|
complex*16 :: integral
|
|
|
|
|
|
|
|
|
|
do k = 1, mo_num
|
|
|
|
|
do i = 1, mo_num
|
|
|
|
|
do j = 1, mo_num
|
|
|
|
|
l = j
|
2020-02-12 01:23:34 +01:00
|
|
|
integral = get_two_e_integral_complex(i,j,k,l,mo_integrals_map,mo_integrals_map_2)
|
2020-03-06 18:09:29 +01:00
|
|
|
big_array_coulomb_integrals_complex(j,i,k) = integral
|
2020-02-05 21:21:28 +01:00
|
|
|
l = j
|
2020-02-12 01:23:34 +01:00
|
|
|
integral = get_two_e_integral_complex(i,j,l,k,mo_integrals_map,mo_integrals_map_2)
|
2020-03-06 18:09:29 +01:00
|
|
|
big_array_exchange_integrals_complex(j,i,k) = integral
|
2020-02-05 21:21:28 +01:00
|
|
|
enddo
|
|
|
|
|
enddo
|
|
|
|
|
enddo
|
|
|
|
|
|
|
|
|
|
END_PROVIDER
|
|
|
|
|
|
