4
1
mirror of https://github.com/pfloos/quack synced 2024-07-07 20:05:58 +02:00
quack/src/eDFT/fock_exchange_individual_energy.f90

32 lines
734 B
Fortran
Raw Normal View History

subroutine fock_exchange_individual_energy(nBas,Pw,P,ERI,Ex)
2020-03-31 22:15:45 +02:00
! Compute the Fock exchange potential
implicit none
! Input variables
integer,intent(in) :: nBas
double precision,intent(in) :: Pw(nBas,nBas)
2020-03-31 22:15:45 +02:00
double precision,intent(in) :: P(nBas,nBas)
double precision,intent(in) :: ERI(nBas,nBas,nBas,nBas)
! Local variables
double precision,allocatable :: Fx(:,:)
double precision,external :: trace_matrix
! Output variables
double precision,intent(out) :: Ex
! Compute HF exchange matrix
allocate(Fx(nBas,nBas))
2020-04-02 09:53:36 +02:00
call fock_exchange_potential(nBas,P(:,:),ERI(:,:,:,:),Fx(:,:))
2020-03-31 22:15:45 +02:00
2020-04-02 09:53:36 +02:00
Ex = 0.5d0*trace_matrix(nBas,matmul(P(:,:),Fx(:,:))) !&
2020-03-31 22:15:45 +02:00
end subroutine fock_exchange_individual_energy