#include "gromacs/math/functions.h"
#include "gromacs/math/utilities.h"
#include "gromacs/math/vec.h"
-#include "gromacs/mdlib/ppforceworkload.h"
#include "gromacs/mdtypes/md_enums.h"
+#include "gromacs/mdtypes/simulation_workload.h"
#include "gromacs/nbnxm/atomdata.h"
#include "gromacs/nbnxm/nbnxm.h"
#include "gromacs/nbnxm/pairlist.h"
const nbnxn_atomdata_t *nbat,
const interaction_const_t *iconst,
rvec *shift_vec,
- const gmx::ForceFlags &forceFlags,
+ const gmx::StepWorkload &stepWork,
int clearF,
gmx::ArrayRef<real> f,
real * fshift,
/* Reaction-field */
krsq = iconst->k_rf*rsq;
fscal = qq*(int_bit*rinv - 2*krsq)*rinvsq;
- if (forceFlags.computeEnergy)
+ if (stepWork.computeEnergy)
{
vcoul = qq*(int_bit*rinv + krsq - iconst->c_rf);
}
fscal = qq*(int_bit*rinvsq - fexcl)*rinv;
- if (forceFlags.computeEnergy)
+ if (stepWork.computeEnergy)
{
vcoul = qq*((int_bit - std::erf(iconst->ewaldcoeff_q*r))*rinv - int_bit*iconst->sh_ewald);
}
Vvdw_rep = c12*rinvsix*rinvsix;
fscal += (Vvdw_rep - Vvdw_disp)*rinvsq;
- if (forceFlags.computeEnergy)
+ if (stepWork.computeEnergy)
{
vctot += vcoul;
}
}
- if (forceFlags.computeEnergy)
+ if (stepWork.computeEnergy)
{
ggid = 0;
Vc[ggid] = Vc[ggid] + vctot;