init_enerdata(mtop.groups.groups[SimulationAtomGroupType::EnergyOutput].nr, inputrec->fepvals->n_lambda, enerd);
/* Set up interactive MD (IMD) */
- t_gmx_IMD *imdSession =
- init_IMD(inputrec, cr, ms, &mtop, mdlog,
- MASTER(cr) ? globalState->x.rvec_array() : nullptr,
- filenames.size(), filenames.data(), oenv, mdrunOptions);
+ auto imdSession = makeImdSession(inputrec, cr, wcycle, enerd, ms, &mtop, mdlog,
+ MASTER(cr) ? globalState->x.rvec_array() : nullptr,
+ filenames.size(), filenames.data(), oenv, mdrunOptions);
if (DOMAINDECOMP(cr))
{
enforcedRotation ? enforcedRotation->getLegacyEnfrot() : nullptr,
deform.get(),
mdModules_->outputProvider(),
- inputrec, imdSession, &mtop,
+ inputrec, imdSession.get(), &mtop,
fcd,
globalState.get(),
&observablesHistory,
// clean up cycle counter
wallcycle_destroy(wcycle);
- // Free PME data
+// Free PME data
if (pmedata)
{
gmx_pme_destroy(pmedata);