MC probability bug in PMTRAION flag
Cho S
chot... at gmail.com
Fri Dec 15 03:39:22 UTC 2017
After looking into the code, I noticed these in 'mc_ensembles.F' between
lines 700-708:
IF (nboxes .EQ. 2) THEN
IF (rand .LT. 0.75E0_dp) THEN
box_number = 1
ELSE
box_number = 2
ENDIF
ELSE
box_number = 1
ENDIF
It looks like the pmtraion move is indeed 75% biased toward the first box.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20171214/865af195/attachment.htm>
More information about the CP2K-user
mailing list