}
else
{
- pull_print_group_x(out, pull->dim, &pull->group[pcrd->group[0]]);
+ pull_print_group_x(out, pull->dim, &pull->group[pcrd->group[0]]);
}
}
pull_print_coord_dr(out, pull->dim, pcrd);
fprintf(out, "\n");
}
-void pull_print_output(t_pull *pull, gmx_large_int_t step, double time)
+void pull_print_output(t_pull *pull, gmx_int64_t step, double time)
{
if ((pull->nstxout != 0) && (step % pull->nstxout == 0))
{
const t_pbc *pbc, double t,
dvec dr, double *dev)
{
- static gmx_bool bWarned = FALSE; /* TODO: this should be fixed for thread-safety,
- but is fairly benign */
+ static gmx_bool bWarned = FALSE; /* TODO: this should be fixed for thread-safety,
+ but is fairly benign */
const t_pull_coord *pcrd;
int m;
double ref, drs, inpr;
double dt, double t)
{
- dvec *r_ij; /* x[i] com of i in prev. step. Obeys constr. -> r_ij[i] */
- dvec unc_ij; /* xp[i] com of i this step, before constr. -> unc_ij */
- dvec *rnew; /* current 'new' positions of the groups */
- double *dr_tot; /* the total update of the coords */
- double ref;
- dvec vec;
- double d0, inpr;
- double lambda, rm, mass, invdt = 0;
- gmx_bool bConverged_all, bConverged = FALSE;
- int niter = 0, g, c, ii, j, m, max_iter = 100;
- double a;
- dvec f; /* the pull force */
- dvec tmp, tmp3;
+ dvec *r_ij; /* x[i] com of i in prev. step. Obeys constr. -> r_ij[i] */
+ dvec unc_ij; /* xp[i] com of i this step, before constr. -> unc_ij */
+ dvec *rnew; /* current 'new' positions of the groups */
+ double *dr_tot; /* the total update of the coords */
+ double ref;
+ dvec vec;
+ double d0, inpr;
+ double lambda, rm, mass, invdt = 0;
+ gmx_bool bConverged_all, bConverged = FALSE;
+ int niter = 0, g, c, ii, j, m, max_iter = 100;
+ double a;
+ dvec f; /* the pull force */
+ dvec tmp, tmp3;
t_pull_group *pdyna, *pgrp0, *pgrp1;
t_pull_coord *pcrd;
bAbs = TRUE;
}
}
-
+
fprintf(fplog, "\nWill apply %s COM pulling in geometry '%s'\n",
EPULLTYPE(ir->ePull), EPULLGEOM(pull->eGeom));
fprintf(fplog, "with %d pull coordinate%s and %d group%s\n",