<div dir="ltr">Hi again,<div><br>Problem solved. At least found.</div><div><br>I was not in the MD restart but in the wavefunction restart.</div><div><br></div><div>There is a difference between a wavefunction.restart created by a square number of MPI processes and created by other number.</div><div><br></div><div>One created by a normal amount of MPI processes can be "correctly" read by a square number of MPI processes (let say 128 to 144)</div><div>One created by a square amount of MPI processes can be "correcly" read by a  square number of MPI processes (let say 144 to 36)</div><div>One created by a square amount of MPI processes cannot be  "correctly" read by a square number of MPI processes  (wfn converged to a different one, hence a drift) (let say 64 to 128)</div><div><br></div><div>On our HPC, a non-square MPI processes don't lead to poor performance (scaling 128 to 144 is <i>as expected)</i> but lead to an incompatibility between restart.wfn (and also accuracy I guess)</div><div><br></div><div>Maybe it's how the decomposition of matrices is done. Or integrals.</div><div><br></div><div>Anyway, if i'll stick to a square number of MPI as "specified in the output" even if in term of pure performance it's the same compared to non-square (quid of accuracy?) ,  it is good.</div><div><br></div><div><br>On Wednesday, March 18, 2015 at 1:07:25 PM UTC+1, Rolf David wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div><div>Good morning everyone,<div><br></div><div>I have a QM/MM simulation which runs for 0.5ps (1000steps) > Everything looks 'normal' (well not so much good, but not exploding). Const Qty is fluctuating up and down but not in a direction or the other.</div><div><br></div><div>I want to restart the calculation for another 0.5ps to reach 1ps > I use my last.restart, put it in EXT_RESTART section. Everything is read, coordinate, kin energy, temp, pot energy, const qty. but Const Qty goes up. And only up.</div><div><br></div><div>In fact the picture attached shows it. You see at precisely 1000 something goes nuts in the Const Qty. So I'm inclined to think there is something wrong in my 'restart' job. </div><div><br></div><div>My restart input is the same as my first run except for renaming the project (0-1000 to 1000-2000)</div><div><br></div><div><span style="line-height:18px">&EXT_RESTART</span></div><div><span style="line-height:18px">  </span><span style="line-height:18px">RESTART_DEFAULT TRUE</span></div><div><span style="line-height:18px">  RESTART_FILE_NAME run_0-1000.restart</span><span style="line-height:18px"><br></span></div><div><span style="line-height:18px">&END EXT_RESTART</span></div></div></div><div><span style="line-height:18px"><br></span></div><div><span style="line-height:18px">If you have some idea of what i'm doing wrong (for the restart), maybe RESTART_DEFAULT FALSE and set up RESTART_POS, RESTART_VEL, RESTART_QMMM, RESTART_THERMOSTAT ?</span></div><div><span style="line-height:18px"><br></span></div><div><span style="line-height:18px">Thanks</span></div><div><span style="line-height:18px"><br></span></div></div></blockquote></div></div>