wallcycle_start(wcycle,ewcSTEP);
- GMX_MPE_LOG(ev_timestep1);
-
if (bRerunMD) {
if (rerun_fr.bStep) {
step = rerun_fr.step;
}
}*/
- GMX_MPE_LOG(ev_timestep2);
-
/* We write a checkpoint at this MD step when:
* either at an NS step when we signalled through gs,
* or at the last step (but not when we do not want confout),
(bNS ? GMX_FORCE_NS : 0) | force_flags);
}
- GMX_BARRIER(cr->mpi_comm_mygroup);
-
- /* if (bTCR)
+/* if (bTCR)
{
mu_aver = calc_mu_aver(cr,state->x,mdatoms->chargeA,
mu_tot,&top_global->mols,mdatoms,gnx,grpindex);
fprintf(fplog,sepdvdlformat,"Constraint",0.0,dvdl);
}
enerd->term[F_DHDL_CON] += dvdl;
-
- GMX_MPE_LOG(ev_timestep1);
-
}
/* MRS -- now done iterating -- compute the conserved quantity */
* the update.
* for RerunMD t is read from input trajectory
*/
- GMX_MPE_LOG(ev_output_start);
-
mdof_flags = 0;
if (do_per_step(step,ir->nstxout)) { mdof_flags |= MDOF_X; }
if (do_per_step(step,ir->nstvout)) { mdof_flags |= MDOF_V; }
}
wallcycle_stop(wcycle,ewcTRAJ);
}
- GMX_MPE_LOG(ev_output_finish);
/* kludge -- virial is lost with restart for NPT control. Must restart */
if (bStartingFromCpt && bVV)
so scroll down for that logic */
/* ######### START SECOND UPDATE STEP ################# */
- GMX_MPE_LOG(ev_update_start);
bOK = TRUE;
if (!bRerunMD || rerun_fr.bV || bForceUpdate)
{
unshift_self(graph,state->box,state->x);
}
- GMX_BARRIER(cr->mpi_comm_mygroup);
- GMX_MPE_LOG(ev_update_finish);
-
if (vsite != NULL)
{
wallcycle_start(wcycle,ewcVSITECONSTR);