{
auto v = statePropagatorData_->velocitiesView();
- relax_shell_flexcon(
- fplog_, cr_, ms, isVerbose_, enforcedRotation_, step, inputrec_, imdSession_,
- pull_work_, step == nextNSStep_, static_cast<int>(flags), localTopology_, constr_,
- energyData_->enerdata(), statePropagatorData_->localNumAtoms(), x, v, box, lambda,
- hist, &forces, force_vir, mdAtoms_->mdatoms(), nrnb_, wcycle_, shellfc_, fr_,
- runScheduleWork_, time, energyData_->muTot(), vsite_, ddBalanceRegionHandler_);
+ relax_shell_flexcon(fplog_,
+ cr_,
+ ms,
+ isVerbose_,
+ enforcedRotation_,
+ step,
+ inputrec_,
+ imdSession_,
+ pull_work_,
+ step == nextNSStep_,
+ static_cast<int>(flags),
+ localTopology_,
+ constr_,
+ energyData_->enerdata(),
+ statePropagatorData_->localNumAtoms(),
+ x,
+ v,
+ box,
+ lambda,
+ hist,
+ &forces,
+ force_vir,
+ mdAtoms_->mdatoms(),
+ nrnb_,
+ wcycle_,
+ shellfc_,
+ fr_,
+ runScheduleWork_,
+ time,
+ energyData_->muTot(),
+ vsite_,
+ ddBalanceRegionHandler_);
nShellRelaxationSteps_++;
}
else
Awh* awh = nullptr;
gmx_edsam* ed = nullptr;
- do_force(fplog_, cr_, ms, inputrec_, awh, enforcedRotation_, imdSession_, pull_work_, step,
- nrnb_, wcycle_, localTopology_, box, x, hist, &forces, force_vir,
- mdAtoms_->mdatoms(), energyData_->enerdata(), lambda, fr_, runScheduleWork_, vsite_,
- energyData_->muTot(), time, ed, static_cast<int>(flags), ddBalanceRegionHandler_);
+ do_force(fplog_,
+ cr_,
+ ms,
+ inputrec_,
+ awh,
+ enforcedRotation_,
+ imdSession_,
+ pull_work_,
+ step,
+ nrnb_,
+ wcycle_,
+ localTopology_,
+ box,
+ x,
+ hist,
+ &forces,
+ force_vir,
+ mdAtoms_->mdatoms(),
+ energyData_->enerdata(),
+ lambda,
+ fr_,
+ runScheduleWork_,
+ vsite_,
+ energyData_->muTot(),
+ time,
+ ed,
+ static_cast<int>(flags),
+ ddBalanceRegionHandler_);
}
energyData_->addToForceVirial(force_vir, step);
}
{
const bool isVerbose = legacySimulatorData->mdrunOptions.verbose;
const bool isDynamicBox = inputrecDynamicBox(legacySimulatorData->inputrec);
- return builderHelper->storeElement(std::make_unique<ForceElement>(
- statePropagatorData, energyData, freeEnergyPerturbationData, isVerbose, isDynamicBox,
- legacySimulatorData->fplog, legacySimulatorData->cr, legacySimulatorData->inputrec,
- legacySimulatorData->mdAtoms, legacySimulatorData->nrnb, legacySimulatorData->fr,
- legacySimulatorData->wcycle, legacySimulatorData->runScheduleWork, legacySimulatorData->vsite,
- legacySimulatorData->imdSession, legacySimulatorData->pull_work, legacySimulatorData->constr,
- legacySimulatorData->top_global, legacySimulatorData->enforcedRotation));
+ return builderHelper->storeElement(
+ std::make_unique<ForceElement>(statePropagatorData,
+ energyData,
+ freeEnergyPerturbationData,
+ isVerbose,
+ isDynamicBox,
+ legacySimulatorData->fplog,
+ legacySimulatorData->cr,
+ legacySimulatorData->inputrec,
+ legacySimulatorData->mdAtoms,
+ legacySimulatorData->nrnb,
+ legacySimulatorData->fr,
+ legacySimulatorData->wcycle,
+ legacySimulatorData->runScheduleWork,
+ legacySimulatorData->vsite,
+ legacySimulatorData->imdSession,
+ legacySimulatorData->pull_work,
+ legacySimulatorData->constr,
+ legacySimulatorData->top_global,
+ legacySimulatorData->enforcedRotation));
}
} // namespace gmx