void print_time(FILE *out,
gmx_walltime_accounting_t walltime_accounting,
- gmx_int64_t step,
+ int64_t step,
t_inputrec *ir,
const t_commrec *cr)
{
static void print_large_forces(FILE *fp,
const t_mdatoms *md,
const t_commrec *cr,
- gmx_int64_t step,
+ int64_t step,
real forceTolerance,
const rvec *x,
const rvec *f)
bool nonFinite = !std::isfinite(force2);
if (force2 >= force2Tolerance || nonFinite)
{
- fprintf(fp, "step %" GMX_PRId64 " atom %6d x %8.3f %8.3f %8.3f force %12.5e\n",
+ fprintf(fp, "step %" PRId64 " atom %6d x %8.3f %8.3f %8.3f force %12.5e\n",
step,
ddglatnr(cr->dd, i), x[i][XX], x[i][YY], x[i][ZZ], std::sqrt(force2));
}
* the printing on other ranks. But we can only avoid that with
* an expensive MPI barrier that we would need at each step.
*/
- gmx_fatal(FARGS, "At step %" GMX_PRId64 " detected non-finite forces on %ju atoms", step, numNonFinite);
+ gmx_fatal(FARGS, "At step %" PRId64 " detected non-finite forces on %ju atoms", step, numNonFinite);
}
}
static void post_process_forces(const t_commrec *cr,
- gmx_int64_t step,
+ int64_t step,
t_nrnb *nrnb,
gmx_wallcycle_t wcycle,
const gmx_localtop_t *top,
gmx_enerdata_t *enerd,
int flags, int ilocality,
int clearF,
- gmx_int64_t step,
+ int64_t step,
t_nrnb *nrnb,
gmx_wallcycle_t wcycle)
{
* \param[in] enerd The energy data; the non-bonded group energies need to be added to enerd.term[F_EPOT] before calling this routine
* \param[in] inputrec The input record
*/
-static void checkPotentialEnergyValidity(gmx_int64_t step,
+static void checkPotentialEnergyValidity(int64_t step,
const gmx_enerdata_t &enerd,
const t_inputrec &inputrec)
{
if (energyIsNotFinite || (averageKineticEnergy > 0 &&
enerd.term[F_EPOT] > c_thresholdFactor*averageKineticEnergy))
{
- gmx_fatal(FARGS, "Step %" GMX_PRId64 ": The total potential energy is %g, which is %s. The LJ and electrostatic contributions to the energy are %g and %g, respectively. A %s potential energy can be caused by overlapping interactions in bonded interactions or very large%s coordinate values. Usually this is caused by a badly- or non-equilibrated initial configuration, incorrect interactions or parameters in the topology.",
+ gmx_fatal(FARGS, "Step %" PRId64 ": The total potential energy is %g, which is %s. The LJ and electrostatic contributions to the energy are %g and %g, respectively. A %s potential energy can be caused by overlapping interactions in bonded interactions or very large%s coordinate values. Usually this is caused by a badly- or non-equilibrated initial configuration, incorrect interactions or parameters in the topology.",
step,
enerd.term[F_EPOT],
energyIsNotFinite ? "not finite" : "extremely high",
const t_inputrec *inputrec,
gmx::Awh *awh,
gmx_enfrot *enforcedRotation,
- gmx_int64_t step,
+ int64_t step,
double t,
gmx_wallcycle_t wcycle,
ForceProviders *forceProviders,
static inline void launchGpuRollingPruning(const t_commrec *cr,
const nonbonded_verlet_t *nbv,
const t_inputrec *inputrec,
- const gmx_int64_t step)
+ const int64_t step)
{
/* We should not launch the rolling pruning kernel at a search
* step or just before search steps, since that's useless.
const t_inputrec *inputrec,
gmx::Awh *awh,
gmx_enfrot *enforcedRotation,
- gmx_int64_t step,
+ int64_t step,
t_nrnb *nrnb,
gmx_wallcycle_t wcycle,
const gmx_localtop_t *top,
const t_inputrec *inputrec,
gmx::Awh *awh,
gmx_enfrot *enforcedRotation,
- gmx_int64_t step,
+ int64_t step,
t_nrnb *nrnb,
gmx_wallcycle_t wcycle,
gmx_localtop_t *top,
const t_inputrec *inputrec,
gmx::Awh *awh,
gmx_enfrot *enforcedRotation,
- gmx_int64_t step,
+ int64_t step,
t_nrnb *nrnb,
gmx_wallcycle_t wcycle,
gmx_localtop_t *top,
t_state *state)
{
int i, m, start, end;
- gmx_int64_t step;
+ int64_t step;
real dt = ir->delta_t;
real dvdl_dum;
rvec *savex;