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