2019-07-04 16:16:57 +02:00
BEGIN_PROVIDER [ double precision , state_av_act_two_rdm_openmp_alpha_alpha_mo , ( n_act_orb , n_act_orb , n_act_orb , n_act_orb ) ]
implicit none
double precision , allocatable :: state_weights ( : )
BEGIN_DOC
! state_av_act_two_rdm_openmp_alpha_alpha_mo(i,j,k,l) = state average physicist two-body rdm restricted to the ACTIVE indices for alpha-alpha electron pairs
! = <Psi| a^{\dagger}_i a^{\dagger}_j a_l a_k |Psi>
END_DOC
allocate ( state_weights ( N_states ) )
2019-10-25 01:30:30 +02:00
state_weights = state_average_weight
2019-07-04 16:16:57 +02:00
integer :: ispin
! condition for alpha/beta spin
ispin = 1
state_av_act_two_rdm_openmp_alpha_alpha_mo = 0.D0
call orb_range_two_rdm_state_av_openmp ( state_av_act_two_rdm_openmp_alpha_alpha_mo , n_act_orb , n_act_orb , list_act , state_weights , ispin , psi_coef , size ( psi_coef , 2 ) , size ( psi_coef , 1 ) )
END_PROVIDER
BEGIN_PROVIDER [ double precision , state_av_act_two_rdm_openmp_beta_beta_mo , ( n_act_orb , n_act_orb , n_act_orb , n_act_orb ) ]
implicit none
double precision , allocatable :: state_weights ( : )
BEGIN_DOC
! state_av_act_two_rdm_openmp_beta_beta_mo(i,j,k,l) = state average physicist two-body rdm restricted to the ACTIVE indices for beta-beta electron pairs
! = <Psi| a^{\dagger}_i a^{\dagger}_j a_l a_k |Psi>
END_DOC
allocate ( state_weights ( N_states ) )
2019-10-25 01:30:30 +02:00
state_weights = state_average_weight
2019-07-04 16:16:57 +02:00
integer :: ispin
! condition for alpha/beta spin
ispin = 2
state_av_act_two_rdm_openmp_beta_beta_mo = 0.d0
call orb_range_two_rdm_state_av_openmp ( state_av_act_two_rdm_openmp_beta_beta_mo , n_act_orb , n_act_orb , list_act , state_weights , ispin , psi_coef , size ( psi_coef , 2 ) , size ( psi_coef , 1 ) )
END_PROVIDER
BEGIN_PROVIDER [ double precision , state_av_act_two_rdm_openmp_alpha_beta_mo , ( n_act_orb , n_act_orb , n_act_orb , n_act_orb ) ]
implicit none
double precision , allocatable :: state_weights ( : )
BEGIN_DOC
! state_av_act_two_rdm_openmp_alpha_beta_mo(i,j,k,l) = state average physicist two-body rdm restricted to the ACTIVE indices for alpha-beta electron pairs
! = <Psi| a^{\dagger}_{i,alpha} a^{\dagger}_{j,beta} a_{l,beta} a_{k,alpha} |Psi>
END_DOC
allocate ( state_weights ( N_states ) )
2019-10-25 01:30:30 +02:00
state_weights = state_average_weight
2019-07-04 16:16:57 +02:00
integer :: ispin
! condition for alpha/beta spin
print * , ''
print * , ''
print * , ''
print * , 'providint state_av_act_two_rdm_openmp_alpha_beta_mo '
ispin = 3
print * , 'ispin = ' , ispin
state_av_act_two_rdm_openmp_alpha_beta_mo = 0.d0
call orb_range_two_rdm_state_av_openmp ( state_av_act_two_rdm_openmp_alpha_beta_mo , n_act_orb , n_act_orb , list_act , state_weights , ispin , psi_coef , size ( psi_coef , 2 ) , size ( psi_coef , 1 ) )
END_PROVIDER
BEGIN_PROVIDER [ double precision , state_av_act_two_rdm_openmp_spin_trace_mo , ( n_act_orb , n_act_orb , n_act_orb , n_act_orb ) ]
implicit none
BEGIN_DOC
! state_av_act_two_rdm_openmp_spin_trace_mo(i,j,k,l) = state average physicist spin trace two-body rdm restricted to the ACTIVE indices
! The active part of the two-electron energy can be computed as:
!
! \sum_{i,j,k,l = 1, n_act_orb} state_av_act_two_rdm_openmp_spin_trace_mo(i,j,k,l) * < ii jj | kk ll >
!
! with ii = list_act(i), jj = list_act(j), kk = list_act(k), ll = list_act(l)
END_DOC
double precision , allocatable :: state_weights ( : )
allocate ( state_weights ( N_states ) )
2019-10-25 01:30:30 +02:00
state_weights = state_average_weight
2019-07-04 16:16:57 +02:00
integer :: ispin
! condition for alpha/beta spin
ispin = 4
state_av_act_two_rdm_openmp_spin_trace_mo = 0.d0
integer :: i
double precision :: wall_0 , wall_1
call wall_time ( wall_0 )
print * , 'providing the state average TWO-RDM ...'
call orb_range_two_rdm_state_av_openmp ( state_av_act_two_rdm_openmp_spin_trace_mo , n_act_orb , n_act_orb , list_act , state_weights , ispin , psi_coef , size ( psi_coef , 2 ) , size ( psi_coef , 1 ) )
call wall_time ( wall_1 )
print * , 'Time to provide the state average TWO-RDM' , wall_1 - wall_0
END_PROVIDER