as you know, potentials like s-w, eam and tersoff are added as "pair"
potentials. I've added another three-body potential, also in the same
way, and I'd like to calculate E_pot of every atom in the system. I
can't see an easy way to do it, because in pair potentials E_pot of
atoms is calculated using the method single(), which takes as an
argument only distance two between atoms, so it won't work for
Any ideas how it should be done?