mirror of
https://github.com/LCPQ/quantum_package
synced 2024-11-08 15:13:52 +01:00
95 lines
2.2 KiB
FortranFixed
95 lines
2.2 KiB
FortranFixed
|
BEGIN_PROVIDER [ integer, mo_tot_num ]
|
||
|
|
implicit none
|
||
|
|
BEGIN_DOC
|
||
|
|
! Total number of molecular orbitals and the size of the keys corresponding
|
||
|
|
END_DOC
|
||
|
|
PROVIDE ezfio_filename
|
||
|
|
call ezfio_get_mo_basis_mo_tot_num(mo_tot_num)
|
||
|
|
ASSERT (mo_tot_num > 0)
|
||
|
|
END_PROVIDER
|
||
|
|
|
||
|
|
|
||
|
|
BEGIN_PROVIDER [ integer, mo_tot_num_align ]
|
||
|
|
implicit none
|
||
|
|
BEGIN_DOC
|
||
|
|
! Aligned variable for dimensioning of arrays
|
||
|
|
END_DOC
|
||
|
|
integer :: align_double
|
||
|
|
mo_tot_num_align = align_double(mo_tot_num)
|
||
|
|
END_PROVIDER
|
||
|
|
|
||
|
|
|
||
|
|
BEGIN_PROVIDER [ character*(64), mo_label ]
|
||
|
|
implicit none
|
||
|
|
BEGIN_DOC
|
||
|
|
! Label characterizing the MOS (local, canonical, natural, etc)
|
||
|
|
END_DOC
|
||
|
|
logical :: exists
|
||
|
|
PROVIDE ezfio_filename
|
||
|
|
call ezfio_has_mo_basis_mo_label(exists)
|
||
|
|
if (exists) then
|
||
|
|
call ezfio_get_mo_basis_mo_label(mo_label)
|
||
|
|
else
|
||
|
|
mo_label = 'no_label'
|
||
|
|
call ezfio_set_mo_basis_mo_label(mo_label)
|
||
|
|
endif
|
||
|
|
|
||
|
|
END_PROVIDER
|
||
|
|
|
||
|
|
|
||
|
|
BEGIN_PROVIDER [ double precision, mo_energy, (mo_tot_num) ]
|
||
|
|
implicit none
|
||
|
|
BEGIN_DOC
|
||
|
|
! Fock diagonal elements
|
||
|
|
END_DOC
|
||
|
|
PROVIDE ezfio_filename
|
||
|
|
mo_energy = 0.d0
|
||
|
|
call ezfio_get_mo_basis_mo_energy(mo_energy)
|
||
|
|
END_PROVIDER
|
||
|
|
|
||
|
|
|
||
|
|
BEGIN_PROVIDER [ double precision, mo_coef, (ao_num_align,mo_tot_num) ]
|
||
|
|
implicit none
|
||
|
|
BEGIN_DOC
|
||
|
|
! Molecular orbital coefficients on AO basis set
|
||
|
|
! mo_coef(i,j) = coefficient of the ith ao on the jth mo
|
||
|
|
END_DOC
|
||
|
|
integer :: i, j
|
||
|
|
double precision, allocatable :: buffer(:,:)
|
||
|
|
!DIR$ ATTRIBUTES ALIGN : $IRP_ALIGN :: buffer
|
||
|
|
|
||
|
|
PROVIDE ezfio_filename
|
||
|
|
allocate(buffer(ao_num,mo_tot_num))
|
||
|
|
buffer = 0.d0
|
||
|
|
call ezfio_get_mo_basis_mo_coef(buffer)
|
||
|
|
do i=1,mo_tot_num
|
||
|
|
do j=1,ao_num
|
||
|
|
mo_coef(j,i) = buffer(j,i)
|
||
|
|
enddo
|
||
|
|
do j=ao_num+1,ao_num_align
|
||
|
|
mo_coef(j,i) = 0.d0
|
||
|
|
enddo
|
||
|
|
enddo
|
||
|
|
deallocate(buffer)
|
||
|
|
|
||
|
|
END_PROVIDER
|
||
|
|
|
||
|
|
BEGIN_PROVIDER [ double precision, mo_coef_transp, (mo_tot_num_align,ao_num) ]
|
||
|
|
implicit none
|
||
|
|
BEGIN_DOC
|
||
|
|
! Molecular orbital coefficients on AO basis set
|
||
|
|
END_DOC
|
||
|
|
integer :: i, j
|
||
|
|
|
||
|
|
do j=1,ao_num
|
||
|
|
do i=1,mo_tot_num
|
||
|
|
mo_coef_transp(i,j) = mo_coef(j,i)
|
||
|
|
enddo
|
||
|
|
do i=mo_tot_num+1,mo_tot_num_align
|
||
|
|
mo_coef_transp(i,j) = 0.d0
|
||
|
|
enddo
|
||
|
|
enddo
|
||
|
|
|
||
|
|
END_PROVIDER
|
||
|
|
|