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

Bugs in SRMC corrected

Browse Source
This commit is contained in:
Anthony Scemama 2016-06-24 09:28:06 +02:00
parent 59d0186209
commit fa84ab6fa2
2 changed files with 2 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
src/PROPERTIES/properties_general.irp.f
View File

@ -47,18 +47,6 @@ BEGIN_PROVIDER [ double precision, wf_extension ]
SOFT_TOUCH wf_extension_min wf_extension_max SOFT_TOUCH wf_extension_min wf_extension_max
END_PROVIDER END_PROVIDER
BEGIN_PROVIDER [ double precision, srmc_pop_weight ]
implicit none
BEGIN_DOC
! Weight of the SRMC population
END_DOC
srmc_pop_weight = pop_weight_mult
srmc_pop_weight_min = min(srmc_pop_weight,srmc_pop_weight_min)
srmc_pop_weight_max = max(srmc_pop_weight,srmc_pop_weight_max)
SOFT_TOUCH srmc_pop_weight_min srmc_pop_weight_max
END_PROVIDER
BEGIN_PROVIDER [ double precision, drift_mod, (size_drift_mod) ] BEGIN_PROVIDER [ double precision, drift_mod, (size_drift_mod) ]
implicit none implicit none

4
src/SAMPLING/srmc_step.irp.f
View File

@ -188,7 +188,7 @@ END_SHELL
elec_coord(elec_num+1,1) += p*time_step elec_coord(elec_num+1,1) += p*time_step
elec_coord(elec_num+1,2) = E_loc elec_coord(elec_num+1,2) = E_loc
elec_coord(elec_num+1,3) = srmc_weight(i_walk) * pop_weight_mult elec_coord(elec_num+1,3) = srmc_weight(i_walk) * srmc_pop_weight_mult
do l=1,3 do l=1,3
do i=1,elec_num+1 do i=1,elec_num+1
@ -359,7 +359,7 @@ BEGIN_PROVIDER [ double precision, srmc_pop_weight_mult ]
BEGIN_DOC BEGIN_DOC
! Population weight of SRMC ! Population weight of SRMC
END_DOC END_DOC
pop_weight_mult = pop_weight(srmc_projection) srmc_pop_weight_mult = srmc_pop_weight(srmc_projection)
END_PROVIDER END_PROVIDER
BEGIN_PROVIDER [ integer, srmc_projection ] BEGIN_PROVIDER [ integer, srmc_projection ]