Use RNG correctly for v-rescale thermostat

[alexxy/gromacs.git] / src / gromacs / mdlib / coupling.c

diff --git a/src/gromacs/mdlib/coupling.c b/src/gromacs/mdlib/coupling.c
index a08e2dcfad986ed70ca3e3c125f3c0526c744ef1..0898371728aca74a1a835ed6d9f471378154813f 100644 (file)
--- a/src/gromacs/mdlib/coupling.c
+++ b/src/gromacs/mdlib/coupling.c
@@ -1520,7 +1520,7 @@ void vrescale_tcoupl(t_inputrec *ir, gmx_int64_t step,
 
             Ek_new  = vrescale_resamplekin(Ek, Ek_ref, opts->nrdf[i],
                                            opts->tau_t[i]/dt,
-                                           ir->ld_seed, step);
+                                           step, ir->ld_seed);
 
             /* Analytically Ek_new>=0, but we check for rounding errors */
             if (Ek_new <= 0)