1
0
mirror of https://gitlab.com/scemama/qp_plugins_scemama.git synced 2024-09-01 08:43:46 +02:00

Compare commits

..

6 Commits

3 changed files with 9 additions and 6 deletions

View File

@ -72,14 +72,17 @@ subroutine run
close(iunit)
iunit = getunitandopen('W.qp','w')
do l=1,ao_num
do k=1,ao_num
do j=1,ao_num
do i=1,ao_num
do j=l,ao_num
do i=k,ao_num
if (i>=j) then
integral = get_ao_two_e_integral(i,j,k,l,ao_integrals_map)
if (integral /= 0.d0) then
write (iunit,'(4(I5,2X),E22.15)') i,j,k,l, integral
endif
endif
enddo
enddo
enddo

View File

@ -137,8 +137,8 @@ subroutine run
enddo
12 continue
close(iunit)
call ezfio_set_ao_one_e_ints_ao_integrals_e_n(A(1:ao_num, 1:ao_num))
call ezfio_set_ao_one_e_ints_io_ao_integrals_e_n('Read')
call ezfio_set_ao_one_e_ints_ao_integrals_n_e(A(1:ao_num, 1:ao_num))
call ezfio_set_ao_one_e_ints_io_ao_integrals_n_e('Read')
allocate(buffer_i(ao_num**3), buffer_values(ao_num**3))
iunit = getunitandopen('W.qp','r')

View File

@ -113,8 +113,8 @@ subroutine run
enddo
12 continue
close(iunit)
call ezfio_set_mo_one_e_ints_mo_integrals_e_n(A)
call ezfio_set_mo_one_e_ints_io_mo_integrals_e_n('Read')
call ezfio_set_mo_one_e_ints_mo_integrals_n_e(A)
call ezfio_set_mo_one_e_ints_io_mo_integrals_n_e('Read')
iunit = getunitandopen('Wmo.qp','r')
n_integrals=0