for (int ftype = 0; ftype < F_NRE; ftype++)
{
- if (((interaction_function[ftype].flags & IF_BOND) && ftype != F_CONNBONDS)
+ if (((interaction_function[ftype].flags & IF_BOND) && NRAL(ftype) == 2 && ftype != F_CONNBONDS)
|| (interaction_function[ftype].flags & IF_CONSTRAINT) || ftype == F_SETTLE)
{
const InteractionList& il = molt->ilist[ftype];