#include "gromacs/math/vec.h"
#include "gromacs/math/vecdump.h"
#include "gromacs/mdlib/forcerec_threading.h"
-#include "gromacs/mdlib/ppforceworkload.h"
#include "gromacs/mdlib/qmmm.h"
#include "gromacs/mdlib/rf_util.h"
#include "gromacs/mdlib/wall.h"
#include "gromacs/mdtypes/inputrec.h"
#include "gromacs/mdtypes/md_enums.h"
#include "gromacs/mdtypes/mdatom.h"
+#include "gromacs/mdtypes/simulation_workload.h"
#include "gromacs/pbcutil/ishift.h"
#include "gromacs/pbcutil/mshift.h"
#include "gromacs/pbcutil/pbc.h"
const real *lambda,
const t_graph *graph,
const rvec *mu_tot,
- const gmx::ForceFlags &forceFlags,
+ const gmx::StepWorkload &stepWork,
const DDBalanceRegionHandler &ddBalanceRegionHandler)
{
// TODO: Replace all uses of x by const coordinates
t_pbc pbc;
/* Check whether we need to take into account PBC in listed interactions. */
- const auto needPbcForListedForces = fr->bMolPBC && forceFlags.computeListedForces && haveCpuListedForces(*fr, *idef, *fcd);
+ const auto needPbcForListedForces = fr->bMolPBC && stepWork.computeListedForces && haveCpuListedForces(*fr, *idef, *fcd);
if (needPbcForListedForces)
{
/* Since all atoms are in the rectangular or triclinic unit-cell,
forceOutputs,
fr, &pbc, graph, enerd, nrnb, lambda, md, fcd,
DOMAINDECOMP(cr) ? cr->dd->globalAtomIndices.data() : nullptr,
- forceFlags);
+ stepWork);
}
const bool computePmeOnCpu =
{
/* Do reciprocal PME for Coulomb and/or LJ. */
assert(fr->n_tpi >= 0);
- if (fr->n_tpi == 0 || forceFlags.stateChanged)
+ if (fr->n_tpi == 0 || stepWork.stateChanged)
{
int pme_flags = GMX_PME_SPREAD | GMX_PME_SOLVE;
- if (forceFlags.computeForces)
+ if (stepWork.computeForces)
{
pme_flags |= GMX_PME_CALC_F;
}
- if (forceFlags.computeVirial)
+ if (stepWork.computeVirial)
{
pme_flags |= GMX_PME_CALC_ENER_VIR;
}