int gmx_sorient(int argc, char* argv[])
{
t_topology top;
- int ePBC = -1;
+ PbcType pbcType = PbcType::Unset;
t_trxstatus* status;
int natoms;
real t;
bTPS = (opt2bSet("-s", NFILE, fnm) || !opt2bSet("-n", NFILE, fnm) || bCom);
if (bTPS)
{
- read_tps_conf(ftp2fn(efTPS, NFILE, fnm), &top, &ePBC, &xtop, nullptr, box, bCom);
+ read_tps_conf(ftp2fn(efTPS, NFILE, fnm), &top, &pbcType, &xtop, nullptr, box, bCom);
}
/* get index groups */
rmin2 = gmx::square(rmin);
rmax2 = gmx::square(rmax);
- rcut = 0.99 * std::sqrt(max_cutoff2(guess_ePBC(box), box));
+ rcut = 0.99 * std::sqrt(max_cutoff2(guessPbcType(box), box));
if (rcut == 0)
{
rcut = 10 * rmax;
if (bTPS)
{
/* make molecules whole again */
- gpbc = gmx_rmpbc_init(&top.idef, ePBC, natoms);
+ gpbc = gmx_rmpbc_init(&top.idef, pbcType, natoms);
}
/* start analysis of trajectory */
do
gmx_rmpbc(gpbc, natoms, box, x);
}
- set_pbc(&pbc, ePBC, box);
+ set_pbc(&pbc, pbcType, box);
n = 0;
inp = 0;
for (p = 0; (p < nrefgrp); p++)