<div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial"><div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial">Dear QE developers,<div><br></div><div>I am using QE 4.3.2 with PLUMED 1.3 (PWscf with metadynamics), and encountered a problem. I think it may be a bug.</div><div>The problem occurs most of the time, and only occasionally the error did not occur.</div><div>The message was always the same:</div><div><br></div><div><div>     Molecular Dynamics Calculation</div><div> </div><div>      Starting temperature  =  9000.00 K</div><div> </div><div>      temperature is set once at start</div><div> </div><div>      mass Mg               =    24.30</div><div>      mass Si               =    28.09</div><div>      mass O                =    16.00</div><div>      Time step             =    20.00 a.u.,  0.9676 femto-seconds</div><div> </div><div>  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</div><div>      from randy : error #      1038</div><div>      j out of range</div><div>  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</div><div> </div><div>      stopping ...</div></div><div><div><br></div></div><div>I checked the relevant sources, and found that it occurred right in the md_init call.</div><div>The error was reported by randy function in the ramdom_numbers.f90 file when it was first called by set_random_seed.</div><div>In the random_numbers.f90:</div><div><br></div><div><div>     50       IF ( first ) THEN</div><div>     51         &nbsp
;!</div><div>     52          first = .false.</div><div>     53          idum = MOD( ic - idum, m )</div><div>     54          !</div></div><div><br></div><div>ic = 150889, if you pass a random seed that is larger than ic, then idum will be a negative number after line 53.</div><div>I wrote a little program to call set_random_seed, most of the time idum will be a negative number and the test in</div><div>randy</div><div>   63       IF( j > ntab .OR. j <  1 ) call errore('randy','j out of range',j)</div><div>fails. j will be a negative number when I test this.</div><div><br></div><div>But in actual run, j is always 1038 whenever the error occurs, which does not match what I see in my little test program. It is very strange.</div><div><br></div><div>If I remove the line 53, everything seems OK in my test program but I am not sure if this line should be corrected in QE.</div><div>So is it a bug in QE, or I just made some mistakes (the random seed can be negative, and the error is caused by other problems)?</div><div><br></div><div>I work on a Linux PC cluster with intel compilers and openmpi, and QE is linked with the mkl library.</div><div><br></div><div>I am not a subscriber of this mail list. Contact me by email to <a href="mailto:dingmaotu@126.com">dingmaotu@126.com</a>. Thank you!</div><div><br></div><div>Regards,</div><div>Li Ding</div><div>Institute of Geology and Geophysics, Chinese Academy of Science</div></div></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>