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