mirror of
https://github.com/LCPQ/quantum_package
synced 2024-11-14 01:53:55 +01:00
48 lines
1.2 KiB
Fortran
48 lines
1.2 KiB
Fortran
program MRPT_Utils
|
|
implicit none
|
|
read_wf = .True.
|
|
touch read_wf
|
|
! call routine
|
|
! call routine_2
|
|
call routine_3
|
|
end
|
|
|
|
|
|
subroutine routine_3
|
|
implicit none
|
|
integer :: i
|
|
!provide fock_virt_total_spin_trace
|
|
provide delta_ij
|
|
|
|
print *, 'N_det = ', N_det
|
|
print *, 'N_states = ', N_states
|
|
do i = 1, N_States
|
|
print*,'State',i
|
|
write(*,'(A12,X,I3,A3,XX,F16.10)') ' PT2 ', i,' = ', second_order_pt_new(i)
|
|
write(*,'(A12,X,I3,A3,XX,F16.09)') ' E ', i,' = ', CI_energy(i)
|
|
write(*,'(A12,X,I3,A3,XX,F16.09)') ' E+PT2 ', i,' = ', CI_energy(i)+second_order_pt_new(i)
|
|
write(*,'(A12,X,I3,A3,XX,F16.09)') ' E dressed ', i,' = ', CI_dressed_pt2_new_energy(i)
|
|
write(*,'(A12,X,I3,A3,XX,F16.09)') ' S^2 ', i,' = ', CI_dressed_pt2_new_eigenvectors_s2(i)
|
|
enddo
|
|
|
|
end
|
|
|
|
subroutine routine_2
|
|
implicit none
|
|
integer :: i
|
|
do i = 1, n_core_inact_orb
|
|
print*,fock_core_inactive_total(i,1,1),fock_core_inactive(i)
|
|
enddo
|
|
double precision :: accu
|
|
accu = 0.d0
|
|
do i = 1, n_act_orb
|
|
integer :: j_act_orb
|
|
j_act_orb = list_act(i)
|
|
accu += one_body_dm_mo_alpha(j_act_orb,j_act_orb,1)
|
|
print*,one_body_dm_mo_alpha(j_act_orb,j_act_orb,1),one_body_dm_mo_beta(j_act_orb,j_act_orb,1)
|
|
enddo
|
|
print*,'accu = ',accu
|
|
|
|
end
|
|
|