10
0
mirror of https://github.com/LCPQ/quantum_package synced 2024-12-23 04:43:50 +01:00

Reduced deterministic set in pt2 stoch

This commit is contained in:
Anthony Scemama 2017-05-05 15:24:04 +02:00
parent 49b413f486
commit 78fe5aeda6

View File

@ -543,12 +543,13 @@ end subroutine
comb_step = 1d0/dfloat(comb_teeth) comb_step = 1d0/dfloat(comb_teeth)
first_det_of_comb = 1 first_det_of_comb = 1
do i=1,N_det_generators do i=1,N_det_generators
if(pt2_weight(i)/norm_left < comb_step*.25d0) then if(pt2_weight(i)/norm_left < comb_step) then
first_det_of_comb = i first_det_of_comb = i
exit exit
end if end if
norm_left -= pt2_weight(i) norm_left -= pt2_weight(i)
end do end do
call write_int(6, first_det_of_comb-1, 'Size of deterministic set')
comb_step = (1d0 - pt2_cweight(first_det_of_comb-1)) * comb_step comb_step = (1d0 - pt2_cweight(first_det_of_comb-1)) * comb_step