LCPQ/quantum_package
10
0
Fork 0
mirror of https://github.com/LCPQ/quantum_package synced 2024-12-22 20:35:19 +01:00
Code Issues Releases Wiki Activity

Merge branch 'master' of github.com:scemama/quantum_package

Browse Source 
Conflicts:
	plugins/read_integral/read_integrals_mo.irp.f
This commit is contained in:
Anthony Scemama 2017-06-08 22:16:42 +02:00
commit 7ccd2c962c
3 changed files with 20 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
README.md
View File

@ -5,7 +5,7 @@
Set of quantum chemistry programs and libraries. Set of quantum chemistry programs and libraries.
(under GNU GENERAL PUBLIC LICENSE v2) (under GNU GENERAL PUBLIC LICENSE v2)
For more information, you can visit the [wiki of the project](http://github.com/LCPQ/quantum_package/wiki>), or below for the installation instructions. For more information, you can visit the [wiki of the project](http://github.com/LCPQ/quantum_package/wiki), or below for the installation instructions.

2
plugins/read_integral/read_integrals_mo.irp.f
View File

@ -1,4 +1,4 @@
program print_integrals program read_integrals
PROVIDE ezfio_filename PROVIDE ezfio_filename
call ezfio_set_integrals_monoelec_disk_access_mo_one_integrals("None") call ezfio_set_integrals_monoelec_disk_access_mo_one_integrals("None")
call run call run

25
src/Integrals_Monoelec/mo_mono_ints.irp.f
View File

@ -3,13 +3,24 @@ BEGIN_PROVIDER [ double precision, mo_mono_elec_integral,(mo_tot_num_align,mo_to
integer :: i,j,n,l integer :: i,j,n,l
BEGIN_DOC BEGIN_DOC
! array of the mono electronic hamiltonian on the MOs basis : ! array of the mono electronic hamiltonian on the MOs basis :
! sum of the kinetic and nuclear electronic potential ! sum of the kinetic and nuclear electronic potential (and pseudo potential if needed)
END_DOC END_DOC
print*,'Providing the mono electronic integrals' print*,'Providing the mono electronic integrals'
do j = 1, mo_tot_num
do i = 1, mo_tot_num IF (do_pseudo) THEN
mo_mono_elec_integral(i,j) = mo_nucl_elec_integral(i,j) + & do j = 1, mo_tot_num
mo_kinetic_integral(i,j) + mo_pseudo_integral(i,j) do i = 1, mo_tot_num
enddo mo_mono_elec_integral(i,j) = mo_nucl_elec_integral(i,j) + mo_kinetic_integral(i,j) &
enddo + mo_pseudo_integral(i,j)
enddo
enddo
ELSE
do j = 1, mo_tot_num
do i = 1, mo_tot_num
mo_mono_elec_integral(i,j) = mo_nucl_elec_integral(i,j) + mo_kinetic_integral(i,j)
enddo
enddo
END IF
END_PROVIDER END_PROVIDER