LCPQ/quantum_package
10
0
Fork 0
mirror of https://github.com/LCPQ/quantum_package synced 2024-12-25 05:43:47 +01:00
Code Issues Releases Wiki Activity

potential numerical precision bug

Browse Source
This commit is contained in:
Yann Garniron 2018-09-04 14:27:10 +02:00
parent 9a0f900d8c
commit 34d9fa0165
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
plugins/dress_zmq/dress_stoch_routines.irp.f
View File

@ -37,11 +37,11 @@ logical function testTeethBuilding(minF, N)
allocate(tilde_w(N_det_generators), tilde_cW(0:N_det_generators)) allocate(tilde_w(N_det_generators), tilde_cW(0:N_det_generators))
tilde_cW(0) = 0d0 tilde_cW(0) = 0d0
do i=1,N_det_generators do i=1,N_det_generators
tilde_w(i) = psi_coef_generators(i,dress_stoch_istate)**2 tilde_w(i) = psi_coef_generators(i,dress_stoch_istate)**2
tilde_cW(i) = tilde_cW(i-1) + tilde_w(i) tilde_cW(i) = tilde_cW(i-1) + tilde_w(i)
enddo enddo
tilde_cW(N_det_generators) = 1d0
n0 = 0 n0 = 0
do do
@ -130,7 +130,7 @@ END_PROVIDER
dress_M_m(i) = U * (((i*i)+i)/2) dress_M_m(i) = U * (((i*i)+i)/2)
end do end do
dress_M_m(1) = 1 dress_M_m(1) = min(dress_M_m(1), 2)
dress_M_m(dress_N_cp_max) = N_det_generators+1 dress_M_m(dress_N_cp_max) = N_det_generators+1
do i=1,N_j do i=1,N_j
@ -556,6 +556,7 @@ end function
tilde_w(i) = psi_coef_generators(i,dress_stoch_istate)**2 tilde_w(i) = psi_coef_generators(i,dress_stoch_istate)**2
tilde_cW(i) = tilde_cW(i-1) + tilde_w(i) tilde_cW(i) = tilde_cW(i-1) + tilde_w(i)
enddo enddo
tilde_cW(N_det_generators) = 1d0
pt2_n_0(1) = 0 pt2_n_0(1) = 0
do do