const char *ffdir)
{
#define MAXATOMSPERRESIDUE 16
- int k, m, i0, ni0, whatres, resind, add_shift, nvsite, nadd;
+ int k, m, i0, ni0, whatres, add_shift, nvsite, nadd;
int ai, aj, ak, al;
int nrfound = 0, needed, nrbonds, nrHatoms, Heavy, nrheavies, tpM, tpHeavy;
int Hatoms[4], heavies[4];
int *o2n, *newvsite_type, *newcgnr, ats[MAXATOMSPERRESIDUE];
t_atom *newatom;
char ***newatomname;
- char *resnm = nullptr;
int cmplength;
bool isN, planarN, bFound;
/* generate vsite constructions */
/* loop over all atoms */
- resind = -1;
+ int resind = -1;
for (int i = 0; (i < at->nr); i++)
{
if (at->atom[i].resind != resind)
{
resind = at->atom[i].resind;
- resnm = *(at->resinfo[resind].name);
}
+ const char *resnm = *(at->resinfo[resind].name);
/* first check for aromatics to virtualize */
/* don't waste our effort on DNA, water etc. */
/* Only do the vsite aromatic stuff when we reach the
}
if (bWARNING)
{
- fprintf(stderr,
- "Warning: cannot convert atom %d %s (bound to a heavy atom "
- "%s with \n"
- " %d bonds and %d bound hydrogens atoms) to virtual site\n",
- i+1, *(at->atomname[i]), tpname, nrbonds, nrHatoms);
+ gmx_fatal(FARGS, "Cannot convert atom %d %s (bound to a heavy atom "
+ "%s with \n"
+ " %d bonds and %d bound hydrogens atoms) to virtual site\n",
+ i+1, *(at->atomname[i]), tpname, nrbonds, nrHatoms);
}
if (bAddVsiteParam)
{