Prepare for virial calculation.
Buggy virial calculation.

The per-iteration values of the virial seem to be OK.  The summation is
somehow wrong.  Probably differences in neighborlist handling as
compared to the LAMMPS implementation.