#include "gromacs/mdlib/lincs.h"
#include "gromacs/mdlib/shake.h"
#include "gromacs/mdtypes/inputrec.h"
-#include "gromacs/mdtypes/mdatom.h"
#include "gromacs/pbcutil/pbc.h"
#include "gromacs/topology/idef.h"
#include "gromacs/topology/ifunc.h"
t_inputrec ir_;
//! Local topology
std::unique_ptr<InteractionDefinitions> idef_;
- //! MD atoms
- t_mdatoms md_;
//! Computational time array (normally used to benchmark performance)
t_nrnb nrnb_;
tensor virialScaledRef_;
//! If the free energy is computed
bool compute_dHdLambda_;
+ //! If there are atoms with perturbed mass
+ bool hasMassPerturbed_ = false;
+ //! Lambda value
+ real lambda_ = 0.0;
//! For free energy computation
real dHdLambda_;
//! For free energy computation (reference value)