int npargs;
npargs = asize(pa);
- if (!parse_common_args(&argc, argv, PCA_CAN_VIEW | PCA_CAN_BEGIN | PCA_CAN_END, NFILE, fnm,
- npargs, pa, asize(desc), desc, 0, nullptr, &oenv))
+ if (!parse_common_args(
+ &argc, argv, PCA_CAN_VIEW | PCA_CAN_BEGIN | PCA_CAN_END, NFILE, fnm, npargs, pa, asize(desc), desc, 0, nullptr, &oenv))
{
return 0;
}
}
if (bORT)
{
- fort = xvgropen(opt2fn("-ort", NFILE, fnm), "Calculated orientations",
- "Time (ps)", "", oenv);
+ fort = xvgropen(
+ opt2fn("-ort", NFILE, fnm), "Calculated orientations", "Time (ps)", "", oenv);
if (bOrinst && output_env_get_print_xvgr_codes(oenv))
{
fprintf(fort, "%s", orinst_sub);
}
if (bODT)
{
- fodt = xvgropen(opt2fn("-odt", NFILE, fnm), "Orientation restraint deviation",
- "Time (ps)", "", oenv);
+ fodt = xvgropen(opt2fn("-odt", NFILE, fnm),
+ "Orientation restraint deviation",
+ "Time (ps)",
+ "",
+ oenv);
if (bOrinst && output_env_get_print_xvgr_codes(oenv))
{
fprintf(fodt, "%s", orinst_sub);
xvgr_legend(out_disre, 2, drleg, oenv);
if (bDRAll)
{
- fp_pairs = xvgropen(opt2fn("-pairs", NFILE, fnm), "Pair Distances", "Time (ps)",
- "Distance (nm)", oenv);
+ fp_pairs = xvgropen(
+ opt2fn("-pairs", NFILE, fnm), "Pair Distances", "Time (ps)", "Distance (nm)", oenv);
if (output_env_get_print_xvgr_codes(oenv))
{
fprintf(fp_pairs, "@ subtitle \"averaged (tau=%g) and instantaneous\"\n", ir->dr_tau);
gmx_fatal(FARGS,
"Number of disre pairs in the energy file (%d) does not match the "
"number in the run input file (%d)\n",
- ndisre, ilist.size() / 3);
+ ndisre,
+ ilist.size() / 3);
}
snew(pairleg, ndisre);
int molb = 0;
k = fa[3 * i + 2];
mtopGetAtomAndResidueName(topInfo.mtop(), j, &molb, &anm_j, &resnr_j, &resnm_j, nullptr);
mtopGetAtomAndResidueName(topInfo.mtop(), k, &molb, &anm_k, &resnr_k, &resnm_k, nullptr);
- sprintf(pairleg[i], "%d %s %d %s (%d)", resnr_j, anm_j, resnr_k, anm_k,
+ sprintf(pairleg[i],
+ "%d %s %d %s (%d)",
+ resnr_j,
+ anm_j,
+ resnr_k,
+ anm_k,
ip[fa[3 * i]].disres.label);
}
set = select_it(ndisre, pairleg, &nset);
}
/* Subtract bounds from distances, to calculate violations */
- calc_violations(disre_rt, disre_rm3tav, nbounds, pair, bounds, violaver,
- &sumt, &sumaver);
+ calc_violations(
+ disre_rt, disre_rm3tav, nbounds, pair, bounds, violaver, &sumt, &sumaver);
fprintf(out_disre, " %8.4f %8.4f\n", sumaver, sumt);
if (bDRAll)
gmx_fatal(FARGS,
"Number of orientation restraints in energy file (%d) does "
"not match with the topology (%d)",
- blk->sub[0].nr, nor);
+ blk->sub[0].nr,
+ nor);
}
if (bORA || bODA)
{
gmx_fatal(FARGS,
"Number of orientation experiments in energy file (%d) does "
"not match with the topology (%d)",
- blk->sub[0].nr / 12, nex);
+ blk->sub[0].nr / 12,
+ nex);
}
fprintf(foten, " %10f", fr.t);
for (i = 0; i < nex; i++)
}
if (bORA)
{
- FILE* out = xvgropen(opt2fn("-ora", NFILE, fnm), "Average calculated orientations",
- "Restraint label", "", oenv);
+ FILE* out = xvgropen(
+ opt2fn("-ora", NFILE, fnm), "Average calculated orientations", "Restraint label", "", oenv);
if (bOrinst && output_env_get_print_xvgr_codes(oenv))
{
fprintf(out, "%s", orinst_sub);
}
if (bODA)
{
- FILE* out = xvgropen(opt2fn("-oda", NFILE, fnm), "Average restraint deviation",
- "Restraint label", "", oenv);
+ FILE* out = xvgropen(
+ opt2fn("-oda", NFILE, fnm), "Average restraint deviation", "Restraint label", "", oenv);
if (bOrinst && output_env_get_print_xvgr_codes(oenv))
{
fprintf(out, "%s", orinst_sub);
}
if (bODR)
{
- FILE* out = xvgropen(opt2fn("-odr", NFILE, fnm), "RMS orientation restraint deviations",
- "Restraint label", "", oenv);
+ FILE* out = xvgropen(opt2fn("-odr", NFILE, fnm),
+ "RMS orientation restraint deviations",
+ "Restraint label",
+ "",
+ oenv);
if (bOrinst && output_env_get_print_xvgr_codes(oenv))
{
fprintf(out, "%s", orinst_sub);
if (bDisRe)
{
- analyse_disre(opt2fn("-viol", NFILE, fnm), teller_disre, violaver, bounds, index, pair,
- nbounds, oenv);
+ analyse_disre(opt2fn("-viol", NFILE, fnm), teller_disre, violaver, bounds, index, pair, nbounds, oenv);
}
{
const char* nxy = "-nxy";