*/
#include "gmxpre.h"
+#include "gromacs/legacyheaders/force.h"
+
#include "config.h"
#include <assert.h>
#include <math.h>
#include <string.h>
-#include "gromacs/legacyheaders/typedefs.h"
+#include "gromacs/bonded/bonded.h"
+#include "gromacs/legacyheaders/coulomb.h"
+#include "gromacs/legacyheaders/domdec.h"
+#include "gromacs/legacyheaders/gmx_omp_nthreads.h"
#include "gromacs/legacyheaders/macros.h"
-#include "gromacs/legacyheaders/force.h"
-#include "gromacs/legacyheaders/nonbonded.h"
+#include "gromacs/legacyheaders/mdrun.h"
#include "gromacs/legacyheaders/names.h"
#include "gromacs/legacyheaders/network.h"
-#include "gromacs/legacyheaders/ns.h"
+#include "gromacs/legacyheaders/nonbonded.h"
#include "gromacs/legacyheaders/nrnb.h"
-#include "gromacs/legacyheaders/bondf.h"
-#include "gromacs/legacyheaders/txtdump.h"
-#include "gromacs/legacyheaders/coulomb.h"
+#include "gromacs/legacyheaders/ns.h"
#include "gromacs/legacyheaders/pme.h"
-#include "gromacs/legacyheaders/mdrun.h"
-#include "gromacs/legacyheaders/domdec.h"
#include "gromacs/legacyheaders/qmmm.h"
-#include "gromacs/legacyheaders/gmx_omp_nthreads.h"
-
+#include "gromacs/legacyheaders/txtdump.h"
+#include "gromacs/legacyheaders/typedefs.h"
#include "gromacs/legacyheaders/types/commrec.h"
#include "gromacs/math/vec.h"
#include "gromacs/pbcutil/ishift.h"
t_fcdata *fcd,
gmx_localtop_t *top,
gmx_genborn_t *born,
- t_atomtypes *atype,
gmx_bool bBornRadii,
matrix box,
t_lambda *fepvals,
{
wallcycle_sub_start(wcycle, ewcsBONDED);
calc_bonds(cr->ms,
- idef, x, hist, f, fr, &pbc, graph, enerd, nrnb, lambda, md, fcd,
- DOMAINDECOMP(cr) ? cr->dd->gatindex : NULL, atype, born,
+ idef, (const rvec *) x, hist, f, fr, &pbc, graph, enerd, nrnb, lambda, md, fcd,
+ DOMAINDECOMP(cr) ? cr->dd->gatindex : NULL,
flags);
/* Check if we have to determine energy differences
{
lam_i[j] = (i == 0 ? lambda[j] : fepvals->all_lambda[j][i-1]);
}
- calc_bonds_lambda(idef, x, fr, &pbc, graph, &(enerd->foreign_grpp), enerd->foreign_term, nrnb, lam_i, md,
+ calc_bonds_lambda(idef, (const rvec *) x, fr, &pbc, graph, &(enerd->foreign_grpp), enerd->foreign_term, nrnb, lam_i, md,
fcd, DOMAINDECOMP(cr) ? cr->dd->gatindex : NULL);
sum_epot(&(enerd->foreign_grpp), enerd->foreign_term);
enerd->enerpart_lambda[i] += enerd->foreign_term[F_EPOT];