rvec f[],
const real chargeA[],
const real chargeB[],
- matrix box,
+ const matrix box,
const t_commrec *cr,
int natoms,
matrix lrvir,
return energy;
}
-real ewald_charge_correction(const t_commrec *cr, t_forcerec *fr, real lambda,
- matrix box,
- real *dvdlambda, tensor vir)
+real ewald_charge_correction(const t_commrec *cr,
+ const t_forcerec *fr,
+ const real lambda,
+ const matrix box,
+ real *dvdlambda,
+ tensor vir)
{
real vol, fac, qs2A, qs2B, vc, enercorr;
*
* Copyright (c) 1991-2000, University of Groningen, The Netherlands.
* Copyright (c) 2001-2004, The GROMACS development team.
- * Copyright (c) 2013,2014,2015,2017,2018, by the GROMACS development team, led by
+ * Copyright (c) 2013,2014,2015,2017,2018,2019, by the GROMACS development team, led by
* Mark Abraham, David van der Spoel, Berk Hess, and Erik Lindahl,
* and including many others, as listed in the AUTHORS file in the
* top-level source directory and at http://www.gromacs.org.
rvec f[],
const real chargeA[],
const real chargeB[],
- matrix box,
+ const matrix box,
const t_commrec *cr,
int natoms,
matrix lrvir,
*
* Should only be called on one thread. */
real
-ewald_charge_correction(const t_commrec *cr, t_forcerec *fr, real lambda, matrix box,
- real *dvdlambda, tensor vir);
+ewald_charge_correction(const t_commrec *cr,
+ const t_forcerec *fr,
+ real lambda,
+ const matrix box,
+ real *dvdlambda,
+ tensor vir);
#endif
const real *chargeA, const real *chargeB,
gmx_bool bHaveChargePerturbed,
const rvec x[],
- matrix box, rvec mu_tot[],
+ const matrix box,
+ const rvec mu_tot[],
int ewald_geometry, real epsilon_surface,
rvec *f,
real *Vcorr_q,
const real *chargeA, const real *chargeB,
gmx_bool bHaveChargePerturbed,
const rvec x[],
- matrix box, rvec mu_tot[],
+ const matrix box,
+ const rvec mu_tot[],
int ewald_geometry, real epsilon_surface,
rvec *f,
real *Vcorr_q,
real chargeA[], real chargeB[],
real c6A[], real c6B[],
real sigmaA[], real sigmaB[],
- matrix box, const t_commrec *cr,
+ const matrix box, const t_commrec *cr,
int maxshift_x, int maxshift_y,
t_nrnb *nrnb, gmx_wallcycle *wcycle,
matrix vir_q, matrix vir_lj,
real chargeA[], real chargeB[],
real c6A[], real c6B[],
real sigmaA[], real sigmaB[],
- matrix box, const t_commrec *cr,
+ const matrix box, const t_commrec *cr,
int maxshift_x, int maxshift_y,
t_nrnb *nrnb, gmx_wallcycle *wcycle,
matrix vir_q, matrix vir_lj,
void
do_force_listed(struct gmx_wallcycle *wcycle,
- matrix box,
+ const matrix box,
const t_lambda *fepvals,
const t_commrec *cr,
const gmx_multisim_t *ms,
* on the set of listed interactions */
void
do_force_listed(struct gmx_wallcycle *wcycle,
- matrix box,
+ const matrix box,
const t_lambda *fepvals,
const t_commrec *cr,
const gmx_multisim_t *ms,
t_nrnb *nrnb,
gmx_wallcycle_t wcycle,
const t_mdatoms *md,
- rvec x[],
+ rvec *x,
history_t *hist,
rvec *forceForUseWithShiftForces,
gmx::ForceWithVirial *forceWithVirial,
gmx_enerdata_t *enerd,
t_fcdata *fcd,
- matrix box,
- real *lambda,
+ const matrix box,
+ const real *lambda,
const t_graph *graph,
- rvec mu_tot[],
- int flags,
+ const rvec *mu_tot,
+ const int flags,
const DDBalanceRegionHandler &ddBalanceRegionHandler)
{
/* do QMMM first if requested */
*/
-void do_force_lowlevel(t_forcerec *fr,
- const t_inputrec *ir,
- const t_idef *idef,
- const t_commrec *cr,
- const gmx_multisim_t *ms,
- t_nrnb *nrnb,
- gmx_wallcycle *wcycle,
- const t_mdatoms *md,
- rvec x[],
- history_t *hist,
- rvec f_shortrange[],
- gmx::ForceWithVirial *forceWithVirial,
- gmx_enerdata_t *enerd,
- t_fcdata *fcd,
- matrix box,
- real *lambda,
- const t_graph *graph,
- rvec mu_tot[2],
- int flags,
- const DDBalanceRegionHandler &ddBalanceRegionHandler);
+void
+do_force_lowlevel(t_forcerec *fr,
+ const t_inputrec *ir,
+ const t_idef *idef,
+ const t_commrec *cr,
+ const gmx_multisim_t *ms,
+ t_nrnb *nrnb,
+ gmx_wallcycle *wcycle,
+ const t_mdatoms *md,
+ rvec *x,
+ history_t *hist,
+ rvec *f_shortrange,
+ gmx::ForceWithVirial *forceWithVirial,
+ gmx_enerdata_t *enerd,
+ t_fcdata *fcd,
+ const matrix box,
+ const real *lambda,
+ const t_graph *graph,
+ const rvec *mu_tot,
+ int flags,
+ const DDBalanceRegionHandler &ddBalanceRegionHandler);
/* Call all the force routines */
#endif