91 lines
3.3 KiB
Fortran
91 lines
3.3 KiB
Fortran
BEGIN_PROVIDER [double precision, f_psi_cas_ab, (n_points_final_grid,N_states)]
|
|
&BEGIN_PROVIDER [double precision, on_top_cas_mu_r, (n_points_final_grid,N_states)]
|
|
implicit none
|
|
BEGIN_DOC
|
|
!
|
|
! Function f_{\Psi^B}(r,r) of Eq. (22) of J. Chem. Phys. 149, 194301 (2018) on each point of the grid and for all states and for a CAS wave function
|
|
!
|
|
! Assumes that the wave function in psi_det is developped within an active space defined
|
|
!
|
|
END_DOC
|
|
integer :: ipoint,istate
|
|
double precision :: wall0,wall1,r(3)
|
|
double precision :: f_ii_val_ab,two_bod_dens_ii,f_ia_val_ab,two_bod_dens_ia,f_aa_val_ab,two_bod_dens_aa
|
|
double precision :: accu
|
|
accu = 0.d0
|
|
r = 0.d0
|
|
istate = 1
|
|
! To initialize parallelization
|
|
call give_f_ii_val_ab(r,r,f_ii_val_ab,two_bod_dens_ii)
|
|
call give_f_ia_val_ab(r,r,f_ia_val_ab,two_bod_dens_ia,istate)
|
|
call give_f_aa_val_ab(r,r,f_aa_val_ab,two_bod_dens_aa,istate)
|
|
provide final_grid_points act_2_rdm_ab_mo
|
|
|
|
print*,'Providing f_psi_cas_ab..... '
|
|
|
|
|
|
call wall_time(wall0)
|
|
do istate = 1, N_states
|
|
!$OMP PARALLEL &
|
|
!$OMP DEFAULT (NONE) &
|
|
!$OMP PRIVATE (ipoint,r,f_ii_val_ab,two_bod_dens_ii,f_ia_val_ab,two_bod_dens_ia,f_aa_val_ab,two_bod_dens_aa) &
|
|
!$OMP SHARED (n_points_final_grid,f_psi_cas_ab,on_top_cas_mu_r,final_grid_points,istate)
|
|
!$OMP DO
|
|
do ipoint = 1, n_points_final_grid
|
|
r(1) = final_grid_points(1,ipoint)
|
|
r(2) = final_grid_points(2,ipoint)
|
|
r(3) = final_grid_points(3,ipoint)
|
|
! inactive-inactive part of f_psi(r1,r2)
|
|
call give_f_ii_val_ab(r,r,f_ii_val_ab,two_bod_dens_ii)
|
|
! inactive-active part of f_psi(r1,r2)
|
|
call give_f_ia_val_ab(r,r,f_ia_val_ab,two_bod_dens_ia,istate)
|
|
! active-active part of f_psi(r1,r2)
|
|
call give_f_aa_val_ab(r,r,f_aa_val_ab,two_bod_dens_aa,istate)
|
|
f_psi_cas_ab(ipoint,istate) = f_ii_val_ab + f_ia_val_ab + f_aa_val_ab
|
|
on_top_cas_mu_r(ipoint,istate) = two_bod_dens_ii + two_bod_dens_ia + two_bod_dens_aa
|
|
enddo
|
|
!$OMP END DO
|
|
!$OMP END PARALLEL
|
|
enddo
|
|
call wall_time(wall1)
|
|
print*,'Time to provide f_psi_cas_ab = ',wall1 - wall0
|
|
print*,'accu = ',accu
|
|
|
|
END_PROVIDER
|
|
|
|
BEGIN_PROVIDER [double precision, f_psi_hf_ab, (n_points_final_grid)]
|
|
&BEGIN_PROVIDER [double precision, on_top_hf_mu_r, (n_points_final_grid)]
|
|
implicit none
|
|
BEGIN_DOC
|
|
!
|
|
! Function f_{\Psi^B}(r,r) of Eq. (22) of J. Chem. Phys. 149, 194301 (2018) on each point of the grid for a HF wave function
|
|
!
|
|
END_DOC
|
|
integer :: ipoint
|
|
double precision :: wall0,wall1,r(3),f_HF_val_ab,two_bod_dens
|
|
f_psi_hf_ab = 0.d0
|
|
r = 0.d0
|
|
! To initialize parallelization
|
|
call f_HF_valence_ab(r,r,f_HF_val_ab,two_bod_dens)
|
|
|
|
call wall_time(wall0)
|
|
!$OMP PARALLEL &
|
|
!$OMP DEFAULT (NONE) &
|
|
!$OMP PRIVATE (ipoint,r,f_HF_val_ab,two_bod_dens) &
|
|
!$OMP SHARED (n_points_final_grid,f_psi_hf_ab,on_top_hf_mu_r,final_grid_points)
|
|
!$OMP DO
|
|
do ipoint = 1, n_points_final_grid
|
|
r(1) = final_grid_points(1,ipoint)
|
|
r(2) = final_grid_points(2,ipoint)
|
|
r(3) = final_grid_points(3,ipoint)
|
|
call f_HF_valence_ab(r,r,f_HF_val_ab,two_bod_dens)
|
|
f_psi_hf_ab(ipoint) = f_HF_val_ab
|
|
on_top_hf_mu_r(ipoint) = two_bod_dens
|
|
enddo
|
|
!$OMP END DO
|
|
!$OMP END PARALLEL
|
|
call wall_time(wall1)
|
|
print*,'Time to provide f_psi_hf_ab = ',wall1 - wall0
|
|
|
|
END_PROVIDER
|