- v = do_nonbonded_listed(ftype, nbn, iatoms+nb0, idef->iparams, x, f, fshift,
- pbc, g, lambda, dvdl, md, fr, grpp, global_atom_index);
+ v = do_pairs(ftype, nbn, iatoms+nb0, idef->iparams, x, f, fshift,
+ pbc, g, lambda, dvdl, md, fr, grpp, global_atom_index);