calculate_1contracted_teints_direct Subroutine

public subroutine calculate_1contracted_teints_direct(env, atm, bas, first_bf_of_shell, n_basis, lc, c_teint, accuracy)

returns two-electron integrals, where the first index has been contracted with LCAO orbital coefficients LC

Todo

:parallelization doesn't work her with large basis sets

Arguments

Type IntentOptional Attributes Name
real(kind=long), intent(in) :: env(:)
integer, intent(in) :: atm(:,:)
integer, intent(in) :: bas(:,:)
integer, intent(in) :: first_bf_of_shell(:)
integer, intent(in) :: n_basis
real(kind=long), intent(in) :: lc(:,:)
real(kind=long), intent(out) :: c_teint(size(lc,2),n_basis,n_basis*(n_basis+1)/2)
real(kind=long), optional :: accuracy