LCPQ
/
quantum_package
mirror of https://github.com/LCPQ/quantum_package
10
0
Fork 0
Code Issues Releases Wiki Activity

Fixed bugs with correlation_energy_ratio

This commit is contained in:
Anthony Scemama 2017-05-02 16:46:31 +02:00
parent 963097e5d6
commit 592f978e62
2 changed files with 12 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
plugins/CAS_SD_ZMQ/cassd_zmq.irp.f
View File

@ -51,17 +51,18 @@ program fci_zmq
E_CI_before(1:N_states) = CI_energy(1:N_states)
double precision :: correlation_energy_ratio
correlation_energy_ratio = E_CI_before(1) - hf_energy_ref
correlation_energy_ratio = correlation_energy_ratio / (correlation_energy_ratio + pt2(1))
correlation_energy_ratio = 0.d0
do while ( &
(N_det < N_det_max) .and. &
(maxval(abs(pt2(1:N_states))) > pt2_max) .and. &
(correlation_energy_ratio < correlation_energy_ratio_max) &
(correlation_energy_ratio <= correlation_energy_ratio_max) &
)
correlation_energy_ratio = E_CI_before(1) - hf_energy_ref
correlation_energy_ratio = correlation_energy_ratio / (correlation_energy_ratio + pt2(1))
correlation_energy_ratio = (CI_energy(1) - hf_energy_ref) / &
(E_CI_before(1) + pt2(1) - hf_energy_ref)
correlation_energy_ratio = min(1.d0,correlation_energy_ratio)
print *, 'N_det = ', N_det
print *, 'N_states = ', N_states

12
plugins/Full_CI_ZMQ/fci_zmq.irp.f
View File

@ -19,7 +19,7 @@ program fci_zmq
hf_energy_ref = ref_bitmask_energy
endif
pt2 = 1.d0
pt2 = -huge(1.d0)
threshold_davidson_in = threshold_davidson
threshold_davidson = threshold_davidson_in * 100.d0
SOFT_TOUCH threshold_davidson
@ -51,17 +51,17 @@ program fci_zmq
n_det_before = 0
double precision :: correlation_energy_ratio
correlation_energy_ratio = E_CI_before(1) - hf_energy_ref
correlation_energy_ratio = correlation_energy_ratio / (correlation_energy_ratio + pt2(1))
correlation_energy_ratio = 0.d0
do while ( &
(N_det < N_det_max) .and. &
(maxval(abs(pt2(1:N_states))) > pt2_max) .and. &
(correlation_energy_ratio < correlation_energy_ratio_max) &
(correlation_energy_ratio <= correlation_energy_ratio_max) &
)
correlation_energy_ratio = E_CI_before(1) - hf_energy_ref
correlation_energy_ratio = correlation_energy_ratio / (correlation_energy_ratio + pt2(1))
correlation_energy_ratio = (CI_energy(1) - hf_energy_ref) / &
(E_CI_before(1) + pt2(1) - hf_energy_ref)
correlation_energy_ratio = min(1.d0,correlation_energy_ratio)
print *, 'N_det = ', N_det
print *, 'N_states = ', N_states