LCPQ/qmcchem
10
1
Fork 1
mirror of https://gitlab.com/scemama/qmcchem.git synced 2024-12-23 04:43:30 +01:00
Code Releases Activity

Merge lpqlx139:qmcchem

Browse Source
This commit is contained in:
Anthony Scemama 2021-05-31 14:00:39 +02:00
commit 56a1b47c61
1 changed files with 3 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
src/PROPERTIES/properties_energy.irp.f
View File

@ -252,12 +252,10 @@ BEGIN_PROVIDER [ double precision, E_loc ]
enddo enddo
! Avoid divergence of E_loc and population explosion ! Avoid divergence of E_loc and population explosion
if (do_pseudo) then if (do_pseudo .and. (qmc_method == t_DMC) ) then
double precision :: delta_e double precision :: delta_e
! delta_e = E_loc-E_ref delta_e = E_loc-E_ref
! E_loc = E_ref + erf(delta_e*time_step_sq)/time_step_sq E_loc = E_ref + delta_e * dexp(-dabs(delta_e)*time_step_sq)
E_loc = max(2.d0*E_ref, E_loc)
! continue
endif endif
E_loc_min = min(E_loc,E_loc_min) E_loc_min = min(E_loc,E_loc_min)
E_loc_max = max(E_loc,E_loc_max) E_loc_max = max(E_loc,E_loc_max)