From 1e3e94a33cf943594a848e2f1c04176d220a517d Mon Sep 17 00:00:00 2001 From: Anthony Scemama Date: Tue, 9 Feb 2021 12:02:14 +0100 Subject: [PATCH] Removed dependency of AO1 on AO2 --- src/ao_one_e_ints/EZFIO.cfg | 8 ++++++++ src/ao_one_e_ints/screening.irp.f | 2 +- src/mo_two_e_ints/integrals_3_index.irp.f | 4 ++-- 3 files changed, 11 insertions(+), 3 deletions(-) diff --git a/src/ao_one_e_ints/EZFIO.cfg b/src/ao_one_e_ints/EZFIO.cfg index ed9cdc35..8d4fff57 100644 --- a/src/ao_one_e_ints/EZFIO.cfg +++ b/src/ao_one_e_ints/EZFIO.cfg @@ -98,3 +98,11 @@ type: Threshold doc: Remove linear dependencies when the eigenvalues of the overlap matrix are below this value interface: ezfio,provider,ocaml default: 1.e-6 + +[ao_one_e_integrals_threshold] +type: Threshold +doc: If | (p|q) | < `ao_one_e_integrals_threshold` then (p|q) is zero +interface: ezfio,provider,ocaml +default: 1.e-15 +ezfio_name: threshold_ao + diff --git a/src/ao_one_e_ints/screening.irp.f b/src/ao_one_e_ints/screening.irp.f index 1bbe3c73..3f436743 100644 --- a/src/ao_one_e_ints/screening.irp.f +++ b/src/ao_one_e_ints/screening.irp.f @@ -4,7 +4,7 @@ logical function ao_one_e_integral_zero(i,k) ao_one_e_integral_zero = .False. if (.not.((io_ao_integrals_overlap/='None').or.is_periodic)) then - if (ao_overlap_abs(i,k) < ao_integrals_threshold) then + if (ao_overlap_abs(i,k) < ao_one_e_integrals_threshold) then ao_one_e_integral_zero = .True. return endif diff --git a/src/mo_two_e_ints/integrals_3_index.irp.f b/src/mo_two_e_ints/integrals_3_index.irp.f index 73e31182..4ffb0134 100644 --- a/src/mo_two_e_ints/integrals_3_index.irp.f +++ b/src/mo_two_e_ints/integrals_3_index.irp.f @@ -2,9 +2,9 @@ &BEGIN_PROVIDER [double precision, big_array_exchange_integrals,(mo_num,mo_num, mo_num)] implicit none BEGIN_DOC - ! big_array_coulomb_integrals(i,j) = = (ii|jj) + ! big_array_coulomb_integrals(j,i,k) = = (ik|jj) ! - ! big_array_exchange_integrals(i,j) = = (ij|ij) + ! big_array_exchange_integrals(i,j,k) = = (ij|kj) END_DOC integer :: i,j,k,l double precision :: get_two_e_integral